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CHAPTER I 


INTRODUCTION 

A tangential wall jet is the flow of high velocity 
fluid emanating from a narrow slot and blowing tangentially 
over a rigid wall. The flow of a wall jet mixing with 
an external moving stream has drawn considerable basic 
and applied interest in the past because it incorporates 
the characteristics of both a boundary layer and a free 
jet. Tangential wall jets are generally used in practice: 
(a) to control turbulent boundary layer separation on high 
lift aerofoils and thereby achieve a large increase in 
lift (Williams and Alexander, 1958), (b) to prevent 

separation and improve the pressure recovery in wide 
angle diffusers (Ramaprian, 1969; Nicoll and Ramaprian, 
1970) , (c) to cool a surface exposed to a stream of hot 

gas as in the case of combustion chambers and the exhaust 
nozzles of rocket motors (Samuel and Joubert, 1964; Seban, 
1960; Papell and Trout, 1959), (d) to heat a surface 

exposed to cold temperatures (Wieghardt, 1946) . 

This thesis is concerned with the behavior of 
two-dimensional incompressible turbulent wall jets sub- 
merged in a boundary layer when they are used to prevent 
boundary layer separation on plane surfaces. The main 
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motivation for studying this topic comes from the large 
increase in lift that can be obtained from high lift aero- 
foils if the flow can be kept attached to the aerofoil 
surface by the use of wall jets (Williams and Alexander, 
1958) . The effect of the jet from the blowing slot is to 
increase the kinetic energy of the flow in the boundary 
layer near the surface, thereby enabling it to advance into 
a high adverse pressure gradient region without separating. 

1 . 1 Previous Work 

Wall jets have been investigated quite extensively 
in the past. However, most work is concerned with the case 
where the ratio of jet velocity to free-stream velocity 
is large in the range of 2 to infinity and with a negligible 
upstream boundary layer at the slot. The reason for the 
concentration of effort on the higher velocity ratios is 
mainly because of the fact that the velocity profiles in 
such cases contain only a velocity maximum instead of also 
a minimum and the flow can be analyzed approximately by 
methods of velocity profile similarity. However, the use- 
ful range of velocity ratios lies between 1 and 2, since 
it is impractical to maintain higher velocity ratios, 
especially in supersonic flows. Even in the cases where 
the study of wall jets at low velocity ratios was attempted, 
the momentum deficit of the upstream boundary layer at the 
slot was small. In practicial applications, however, the 
jet usually mixes with a thick upstream boundary layer 
that is approaching separation, giving rise to a velocity 
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profile shown in Fig. 1.1(b) rather than that of a "simple" 
wall jet in Fig. 1.1(a). Therefore, the study of wall 
jets submerged in a thick upstream boundary layer xs very 
useful. This class of flows will be refer-rtt to as "blown 
boundar” layers" since the wall jet is submerged in a 
boundary layer. 

1.1(a) Previous Work on Wall Jets 
With Upstream Boundary Layers 

In describing the previous work, attention will 
be restricted to those authors who have considered a 
wall jet submerged in an upstream boundary layer, which 
was referred to earlier as "blown boundary layex . " Irwin 
(1974) gives a fairly comprehensive description of the 
work in the literature on blown boundary layers. 

Carriere, Eichelbrenner and Poisson-Ouinton 
(1959) appear to have been the first to measure detailed 
mean velocity profiles downstream of the slot. They 
attempted to use an empirical integral method to predict 
the development of the flow. Thomas (1962, 1965) gave 
a crude empirical method determining the blowing momen- 
tum required to prevent separation. Even though Thomas' 
empirical method is simple, it does not seem very sound 
as it is based on very little experimental evidence. 

Bradshaw and Gee (1962) presented mean flow 
measurements in a blown boundary j ayer along with measure- 
ments in simple wall jets on curved and plane surfaces. 

They identified two essentially different modes of 
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"separation, " one involving reversed flow at the wall, 
the other having reversed flow in the wake at the velocity 
minimum. McGahan (1965) carried out mean flew measurements 
in blown boundary layers on a cylinder with its axis 
aligned with the flow direction. The pressure gradient 
was adjusted by sticking paper strips on the outside of 
an outer concentric porous cylinder and by placing an 
end plate to block the flow, at the downstream end. 

McGahan proposed an integral calculation method for pre- 
dicting the flow development when the upstream boundary 
layer was thick. It agreed with his data fairly well 
except near to the point of separation at the wall. 

Gartshore and Newman (1969) described an integral 
calculation method for wall jets in arbitrary pressure 
gradients primarily for the use with small upstream 
boundary layers . However they did attempt to t ccount for 
a large upstream boundary layer by changing the starting 
conditions properly. They also presented measurements 
of mean quantities for wall jets under zero and adverse 
pressure gradients with small upstream boundary layers . 

Kind (1971) extended Gartshore and Newman's 
calculation method to deal with strongly curved wall 
jets developing in the presence of moderately thick up- 
stream boundary layers. The upstream boundary layer was 
accounted for by using a revised method of starting and 
assuming that the stagnation pressure remains constant 
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along streamlines outside the wall jet. He also presented 
mean flow measurements on wall jets on the cylindrical 
trailing edge of a circulation control airfoil. 

Most of the above mentioned previous work on 
blown boundary layers was mainly concerned with high jet 
velocity ratios (greater than 2.0) and measurements of 
mean quantities under small upstream boundary layers . 

Goradia and Colwell (1971) measured mean velocity 
profiles of several wall jets under adverse pressure 
gradients in a two-dimensional diffuser with low jet 
volicity ratios in the neighborhood of 1 to approximately 
2. The velocity data were used to formulate empirical 
relationships between parameters such as the form factor 
H and energy form factor and to derive empirical expres- 
sions for the velocity profiles. The measurements were 
also utilized for the calculations of wall shear and shear 
distribution by numerical methods. English (1970) 
considered the flow over a slotted flap w’hich is equiva- 
lent to considering a wall jet with the total head at the 
slot equal to that of the freestream. The measurements 
of mean velocity and shear stress were made under zero 
and adverse pressure gradients. The slot width was some- 
what larger than that usually used for wall jets, so that 
a region of potential flow often exists in the flow from 
the slot over a large portion of the flap surface. 

Kacker and Whitelaw? (1968, 1971) investigated wall 
jets under a zero pressure gradient with the jet velocity 
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ratios in the range of 0.75 and 2.74 and with a small 
upstream boundary layer. They made measurements of mean 
velocity, turbulence intensities, turbulent shear stress, 
and spectra. They derived the results of tidy viscosity, 
mixing length and Prandtl-Kolmogorov length scale from 
the measurements. However, the restriction to zero 
pressure gradient conditions limits the usefulness of 
the data. The upstream boundary layer was too small to 
have a very significant influence on the flow. 

Ramaprian (1973, 1975) reporte measurements on 
the wall of a conical diffuser with annular injection at 
the entrance and with the jet velocity ratios in the 
range of 1.27 to 2.38. Measurements of mean velocity, 
turbulence intensities, shear stress, and spectra were 
reported. He used Spalding and Patankar's (1967) method 
with a mixing length model for the prediction of his 
flows and the method yielded satisfactory predictions 
of the wall jet development. The empirical constants 

were adjusted to give best agreement with his data. 

Irwin (1974) studied a number of blown boundary 
layers with adverse pressure gradients. Measurements of 
mean velocity, turbulence intensities, shear stress and 
spectra were made for the case of a self-preserving strong 
wall jet under an adverse pressure gradient flow with 
negligible upstream boundary layer. Another case he stud- 
ied was that of a strong wall jet (jet velocity ratio = 3.3) 
under an adverse pressure gradient with a small upstream 
boundary layer. The effect of the upstream boundary layer 
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in this case persisted quite far downstream of the slot, 
but was eventually absorbed by the wall jet. Measurements 
of mean velocity, turbulence intensities and shear stress 
were reported in this case. Irwin also reported measure- 
ments of mean velocity on two strong wall jets (jet velocity 
ratio -3.3 and 2.5) with highly adverse pressure gradients 
and with large upstream boundary layers. Only one case of 
a wall jet under adverse pressure gradient conditions with 
a relatively smaller jet excess velocity (jet velocity 
ratio =1,65) and large upstream boundary layer was reported 
by Irwin. In this case, the wall jet was just sufficient 
to prevent separation, but no detailed measurements of 
turbulence were made. Irwin developed a theoretical 
prediction method applicable to blown boundary layers . 

It essentially uses the computing method of Spalding 
and Patankar (1967, 1969) and the turbulence model pro- 
posed by Launder, Reece and Rodi (1973) along with the 
modifications of the turbulence model to account for the 
effect of the wall on the turbulence and the streamline 
curvature. He reported good predictions of his measure- 
ments in blown boundary layers and also other different 
types of flows, which included isolated wall jets, normal 
boundary layers and curved wall jets. 

The above mentioned previous work on blown 
boundary layers reveals that previous studies of wall 
jets under thick upstream boundary layers under zero and 
adverse pressure gradients are very limited. In some 
cases where attempts have been made, the data were 
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generally limited to the mean quantities rather than 
detailed data on turbulence. Hence, there is a need for 
the study of wall jets under thick upstream boundary 
layers and with low jet velocity ratios, i.e., less than 

2 . 

1.2 Asymmetric Jet Velocity Profile 

Almost all of the previous investigators of wall 
jets have used a uniform jet velocity profile shown in 
Fig. 1.2(a) with negligible upper and lower wall boundary 
layers of the jet nozzle. In cases where the jet bound- 
ary layers were considerable, the velocity profile was 
symmetric about the centerline of the slot as shown in 
Fig. 1.2(b). However, it is interesting to see how the 
flow development is affected if the jet velocity profile 
is made asymmetric as shown in Fig. 1.2(c) instead of 
uniform or symmetric for a given value of the total 
jet momentum. For the same total jet momentum, an 
asymmetric profile will have higher maximum jet velocity 
than a symmetric profile. The idea of using an asymmetric 
velocity profile stems from the following arguments. In 
a practical situation, the upstream boundary layer at 
the slot is the one corresponding to a flow approaching 
separation and has a large deficit of momentum. There- 
fore, it is reasonable to state that a greater part of 
the jet momentum should be made available nearer to the 
slot lip than nearer to the wall to reduce the momentum 
deficit of the upstream boundary layer. However, one 
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might argue that momentum is also needed near the wall 
to prevent separation there. But relatively more momentum 
is required near the slot lip than near the wall and the 
momentum of the jet in the lower half of the jet should 
be sufficient to prevent wall flow separation. In com- 
parison, the symmetric jet velocity profile might be 
able to prevent separation at the bottom wall, but a large 
momentum deficit region might develop in the outer layer 
far downstream which might eventually lead to separation 
at the lower wall. Also, the asymmetric profile should 
result in less frictional losses at the wall as the 
velocity gradients are relatively smaller at the wall. 

Thus, the jet momentum is more effectively utilized 
in the case of an asymmetric jet velocity profile 
in reducing the momentum deficit of the upstream boundary 
layer instead of wasting it on wall friction. 

1.3 Objectives of the Thesis 

The objectives of the present work are as follows: 

1. To obtain experimental data on low jet 
velocity ratio wall jets in thick upstream boundary lay- 
ers advancing into zero and strong adverse streamwise 
pressure gradients. 

2. To obtain detailed turbulence data on the 
type of wall jets considered here in order to aid in the 
future development of turbulence models and prediction 
methods and also to supplement the few existing sets of 
turbulence data for this type of flow. 
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3. To investigate the effectiveness of an 
asymmetric jet velocity profile, as compared to a uniform 
profile, in the control of separation and its influence 
on the development of the flow downstream of the slot. 

4. To predict the present class of wall jets 
using an existing theoretical method for such flows 
and to investigate the effectiveness of an asymmetric 
jet velocity profile from a theoretical basis. 

1.4 Outline of the Thesis 

Chapter 2 describes the experimental aparatus . 

The constructional and performance details of the wind 
tunnel and the wall jet flow system are given as well as 
how the asymmetric jet velocity profile was produced 
at the slot. Chapter 3 is devoted to the description of 
the instrumentation used in obtaining the experimental 
data . 

Chapter 4 gives the experimental results for the 
zero pressure gradient flow. A brief description of the 
flow conditions and the process of setting the zero 
pressure gradient are given. The effect of the asymmetric 
jet velocity profile on the flow development as indicated 
by different measured quantities is described. The 
turbulence data for the zero pressure gradient case were 
studied in more detail in order to be useful as an aid 
in developing turbulence models and prediction methods 
in the future. The spectral data were obtained only 
for the zero pressure gradient flow. The measured 
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quantities included mean velocity, turbulence intensities, 
shear stress, spectra and skin friction. The derived 
quantities included, turbulent kinetic energy, eddy 
viscosity, Prandtl mixing length, Prandtl-Kolmogorov 
length scale, correlation coefficients, dissipation and 
production rates of turbulent kinetic energy and bursting 
periods . 

Chapter 5 is devoted to the experimental results 
for the adverse pressure gradient flow. The flow condi- 
tions and the process of setting the adverse pressure 
gradient are given. Fewer turbulence data are presented 
as compared to the zero pressure gradient flow. The 
pressure gradient was adjusted to be representative of 
the practical situation on aerofoils and jet flaps. The 
effect of the asymmetric jet velocity profile on the flow 
development is described for the case of adverse pressure 
gradient flow. 

Chapter 6 is concerned with the description of the 
prediction method used to predict the present flows . The 
prediction method developed by Irwin (1974) for blown 
boundary layers was used. Only the salient features of 
the method are given. 

Chapter 7 deals with the computed results. The 
prediction results for the present flows are compared with 
the experimental data. A main distinction is made between 
the predictions using the experimental starting conditions 
and the predictions using the "automatic starting pro- 
cedure" of Irwin. The former represents the experimental 
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asymmetric jet velocity profile and the latter represents 
the uniform jet velocity profile. The advantages of an 
asymmetric profile over a uniform profile are discussed 
on a theoretical basis. Finally, the case of a linear 
jet velocity profile was taken to represent an ideal 
asymmetric profile, which had the greatest momentum 
near the slot top. Computations were made with the 
linear velocity profile and compared with the predic- 
tions using a uniform profile for the adverse pressure 
gradi ;nt flow. 

Chapter 8 lists the conclusions of the thesis. 


CHAPTER II 


EXPERIMENTAL APPARATUS 

This chapter deals with the construction and per- 
formance details of the wind tunnel and the wall jet-flow 
system used in the present work. The function of the 
wind tunnel is to supply the test section with low tur- 
bulence air at a given velocity and temperature. The 
function of the wall jet is to supply the secondary 

necessary for tangential flow injection at the wall. 

Figure 2.1 shows the schematic of the wind tunnel 
including the wall jet system. The wind tunnel is an 
open circuit blower type tunnel and uses the ambient 
air. The important components of the wind tunnel include: 
the heat exchanger for temperature control, the blower, 
the plenum chamber, the contraction and the test section. 

The wall jet flow system consists of the wall jet assem- 
bly and the air supply system. 

The different components of the wind tunnel are 
described below in the order that the flow passes through 
them, followed by a description of the wall jet flow system. 

2.1 Description of the Wind Tunnel Components 

2.1.1 Heat Exchanger 

The first component in the path of the airflow 
is the heat exchanger. A passenger car radiator used for 
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this purpose is a finned tube type with a frontal area 
of 60X48 cm. The heating or cooling of the ambient air 
was accomplished by allowing hot or cold water through 
the radiator. By suitably controlling the flow rate of 
water it was possible to control the temperature of the 
air in the test section within ±0.1 6 C at 25°C. The heat 
exchanger was mainly used to cool the air in the present 
experiments. The room air conditioning system, supported 
by two portable heaters each of 1000W capacity was used 
to heat the ambient air. By properly controlling the 
heat input to the heaters, it was possible to control the 
temperature of the air in the test section within ±0.1°C 
at 25 °C . 

2.1.2 Blower System 

The next component in the air circuit is a FARR 
HP-2A class 2 type rear access air filter of 61X61X31 
cm size enclosed in a box covered on all four sides. 

This filter is capable of removing dust and foreign 
particles down to 5 microns size and larger with a 95% 
efficiency. With less efficiency, it filters particles 
down to 2 microns size. A damper was placed between the 
heat exchanger and the filter to control the amount of 
air flow and hence the velocity of the air in the test 
section. The damper was made of a plexiglas sheet that 
slides in an aluminum frame mounted ip the access space 
between the filter frame and the heat exchanger. 
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The air then enters the suction side cf an Aero- 
vent, BIA Aerofoil type, belt driven, centrifugal, gen- 
eral purpose blower. The blower is driven by a Reliance 
1.49 KW , 1970 RPM, open, drip-proof motor. The blower 
can deliver 40m^/min of air at 6.6 .m static pressure. 

The blower wheel is 31 cm in diameter and has aerofoil 
type blades. 

2.1.3 Plenum Chamber 

The air delivered by the blowc. : enters the plenum 
chamber. The primary function of the plenum chamber is 
to suppress any large scale fluctuations produced by the 
blower and to settle the air. The plenum chamber is made 
of 0.7 mm thick sheet metal and is 99 cm wide, 155 cm 
long, and 114 cm high. The plenum chamber is connected 
to the blower on the upstream side and to the ductwork 
on the downstream side through rubber sheets to avoid 
the transmission of vibrations from the blower to the 
downstream ductwork and finally to the test section. 

The inlet and outlet of the plenum chamber, dimensioned 
55X55 cm and 53X53 cm respectively, were facing each 
other to start with. In this case, the air from the 
blower was entering directly into the downstream ductwork, 
unaffected by the plenum chamber. To eliminate this prob- 
lem, a wooden baffle of 114X112X1.25 cm in size was placed 
inside the chamber between the inlet and the outlet. The 
baffle divides the plenum chamber into two equal compart- 
ments joined by about a 114X41 cm gap at each end of 
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the baffle. The function of the baffle is to circumvent 
the air from the blower to the sides, thereby preventing 
the air from entering directly into the downstream duct 
work. The baffle plate fully covers the height of the 
plenum chamber. 

2.1.4 Transition Section 

The ductwork between the plenum chamber and the 
test section is classified as the transition section. 

It includes two contractions and a rectangular duct 
containing the screens and the honeycomb material. The 
purpose of the first contraction is two-fold: 

1. It serves as a transition piece between the 
plenum chamber and the downstream sections, and 

2. It increases the flow velocity and reduces 
the turbulence level. 

The first contraction has the inlet and outlet dimensions 
of 53X53 cm and 39X23 cm respectively and is 46 cm in 
length. It is made of 0.7 mm thick sheet metal. The air 
after passing through the first contraction enters a 
rectangular duct also made of 0.7 mm thick sheet metal 
with dimensions, 46X28X39 cm. This rectangular duct 
houses the honeycomb material and three screens whose 
primary function is to reduce the turbulence level in the 
air stream. The screens and honeycomb were selected 
according to guidelines given by Bradshaw and Pankhurst 
(1964) . The aluminum honeycomb is 3.2 mm in cell size 
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and 4.3 cm in thickness. The honeycomb eliminates the 
large scale turbulence and swirl in the air stream. The 
scale of turbulence is further reduced by three screens 
in series. The first two aluminum screens are 16 (hor- 
izontal) X 18 (vertical) mesh screens with a wire diameter 
of 0.25 mm and an open area ratio 0.63. The third 
aluminum screen is a 32X32 square mesh screen with a wire 
diameter of 0.18 mm and an open area ratio 0.602. The 
screens are placed 14 cm apart to allow the turbulence 
in the wake of each screen to decay before the next screen 
is reached. 

The next component in the transition section is 
a second contraction made of 0.7 mm thick sheet metal 
with inlet and outlet dimensions of 40X24 cm and 10.2X24 
cm, respectively, giving a contraction ratio of nearly 
4 to 1 . This contraction was designed by Simpson and 
Wyatt (1972) according to the design method of Jordinson 
(1961) . It further reduces the turbulence intensity of 
the air stream. Measurements made before the test sec- 
tion was installed indicated a very flat velocity profile 
at the outlet of the second contraction. The velocity 
was virtually constant in the spanwise direction. The 
turbulence intensity u t /U 00 was about 0.2% in the free- 
stream of the test section. 
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2.1 , _ 5_ Test Section 

The test section is 24 cm wide and has a total 
length of 196 cm. The side walls arc 6.4 mm thick 
float type plate glass 22 cm in height and 196 cm in 
length. The bottom wall is 1.9 cm thick hard "fin-form" 
plywood with a very smooth surface finish and is made of 
two parts. The first part is 58.5 cm in length and the 
second part is 128.3 cm. The 9.2 cm gap between them is 
filled by the wall jet, the constructional details of 
which are given later. The leveled bottom wall rests on 
a steel platform bolted to the concrete floor. 

A piece of sandpaper 16.5 cm (length) X 24 cm 
(width) is glued to the bottom wall at the beginning of 
the test section immediately after the contraction outlet. 
This sandpaper is a "NORTON CLOSEKOTE" silicon carbide 
floor sanding paper with 24 grit size (mesh number) and 
with an average particle size of the abrasive grain 
equal to 1.04 mm. Several other types of roughness 
elements such as, (a) 6.4 mm square rod, (b) 2.4 mm 
round rod, (c) a rectangular strip 12.7X3.2 mm and com- 
binations of these were also tried before selecting the 
sandpaper. The purpose of the roughness element is to 
produce a thick turbulent boundary layer at the end of 
the first 58.5 cm of the test section. The velocity 
profile at 45 cm away from the beginning of the test 
section and at the center of the test section was measured 
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with a rectangular mouth boundary layer total head pitot 
tube (Fig. 3.1) , and with each of the above mentioned 
roughness elements in position. The present sandpaper 
roughness element was selected on t’.e basis of maximum 
momentum thickness Reynolds number and maximum strength 
of the wake component (Coles, 1962) obtained from the 
measurements at 45 cm away from the beginning of the 
test section. The reason for having a large momentum 
thickness Reynolds number is because of the fact that 
the present wall jet experiments were planned to be 
carried out with large momentum deficit upstream boundary 
layer meeting the wall jet. Three brass rods of 1.6 mm 
diameter and of proper length are glued to the side and 
top edges of the contraction outlet. These tripping 
devices fix the point of transition on the top and side- 
wall boundary layers. 

2.1.5(a) Top Wall for the Zero Pressure Gradient Studies 
The top wall used for zero pressure gradient 
studies is a 9.5 mm thick plexiglas sheet with access 
holes at several stations to insert the measuring probes 
and can be adjusted to various positions. The edges of 
the top wall are sealed against the glass side walls by 
squeezing foam weather stripping in between them. The 
nominal height of the test section between the contraction 
outlet and the wall jet is 9.6 cm. The height of the 
remaining part of the test section can be adjusted by 
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moving the upper wall to create a zero px-ersure gradient 
flow. The end section of the tunnel was kept open for 
zero pressure gradient studies. 

2.1.5(b) Top Mall for the Adverse Pressure Gradient Studies 

A considerable amount of effort has been spent 
in selecting the proper method for producing an adverse 
pressure gradient without inducing strong three-dimensional 
flow effects. Two of the possible methods were: 

1. A test section with a solid top wall and an 
increasing cross-sectional area with an increasing static 
pressure, and 

2. A test section with a perforated top wall 
which allows bleeding off of a portion of the flow in the 
test section thereby increasing the static pressure 
downstream. The second method was selected in view of 
the limitations on the size of the test section and the 
severity of the required pressure gradient. The pressure 
distribution in this method can be varied by suitably 
covering some of the perforations uniformly across the 
tunnel. The static pressure difference between the 
inside. of the test section and the ambient atmosphere 
necessary to bleed the flow was created by a perforated 
plate attached to the end section of the tunnel. By 
properly covering portions of the end plate, the level of 
static pressure inside the tunnel can be varied. 
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A number of commercially available perforated 
sheets of different materials were considered for the 
top wall of the present application. A 3 . 4 mm thick 
"masonite" perforated sheet with an open area ratio 0.33 
was finally selected for the present application. 

Figure 2.2 shows the perforated sheet in position. 

Figure 2.3 shows one single cell of the perforated plate 
with a scale in the inset. The cells are square in shape 
and are 19 rranin size (center to center) . The perforated 
plate was cut to size in such a way that the test sec- 
tion width is spanned by 12 cells. Two aluminum angles 
(25X13X3 mm) , one on each side were attached by screws 
to two rectangular (13 a 3 mm) aluminum strips lying 
above the aluminum angles with the perforated top wall 
sandwiched in between them. The aluminum angle is placed 
inside the flow in such a way that its longer side is 
perpendicular to the tunnel floor and runs all the way from 
the beginning to the end of the tunnel. This aluminum 
angle serves two purposes : 

1. It acts as a reinforcement for the structurally 
fragile "masonite” top wall, and 

2. It assists in side wall boundary layer control 
by acting as a baffle between the main flow and the side 
wall boundary layer. 

The first 4Jj cells on each side of the center 
of the tunnel at any given streamwise section were open 



Fig. 2.2 Perforated Top Wall in Position for the Adverse 
Pressure Gradient Flow 
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for flow bleeding. The next half cell was covered with 
the aluminum angle. After the aluminum ?• -le, one cell 
was available for sidewall boundary layer bleeding. 

Small plastic mouldings of a channel cross section with 
air tight rubber beading pressed into them were slipped 
on to the edges of the top wall. The top wall was then 
pushed into position with a slight pressure against the 
glass side walls to ensure proper sealing at the edges. 
The perforations in the top wall were also used to in- 
sert the various probes without further disturbing the 
flow . 

The end plate of proper size to span the end 
section of the tunnel was cut from the same perforated 
sheet from which the top wall was made. The end plate 
contained six rows of cells. The static pressure level 
inside the tunnel was controlled by properly covering 
these cells with a commercially available duct tape. 

Fig. 2.4 shows the end plate. 

2,2 Description and Performance of 
The Wall Jet Flow System 

The wall jet flow system can be divided into two 
parts: the air supply system and the wall jet assembly, 

each of which is described below in detail followed by 
the details on the performance of the wall jet flow 
system. 
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2.2.1 Air Supply System 

The air supply to the wall jet is obtained from 

the compressed air supply available in the laboratory. 

The compressed air supply can provide filtered air at a 

6 2 

maximum pressure of 2.76X10 N/m and at the flow rates 
suitable to the present requirements. This compressed 
air is fed to the wall jet through a 13 mm I.D. copper 
tube. A pressure regulator in the line regulates the 
input pressure to the wall jet so as to maintain a constant 
mass flow rate and it is located 178 cm upstream of the 
wall jet inlet along the copper tubing. The last 41 cm 
of the copper tubing has 25 mm I.D. to suit the inlet 
to the wall jet. 

The temperature of the air entering the wall 
jet assembly was always found to be lower than the 
required level. Two electrical heating tapes were 
wrapped around the copper tubing between the pressure 
regulator and the wall jet inlet to heat the compressed 
air. These heating tapes, supplied by the Fischer 
Scientific Co., are 305X1.3 cm in size with a maximum 
attainable temperature of 249 °C. The heating tapes 
operate on 115V supply and have a power output of 340 
Watts. The electrical power input to the tapes is con- 
trolled by a 115V, 10A variable auto transformer. By 
properly adjusting the variable auto transformer, the 
temperature of the air coming out of the wall jet can 
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be controlled to the required value within t 0.3°C at 
25° C. The temperature of the wall jet air was measured 
by a thermometer fitted into the copper tubing near to 
the wall jet assembly inlet with its bulb located in the 
flow. Because of the low residence time for the air, the 
heat losses between the point where the thermometer is 
located and the exit of the wall jet were found to be very 
small. This was indicated by the air temperature differ- 
ence between the inlet and the exit of the wall jet assembly 
which was found to be of the order of 0.3°C. Because of 
the high mass flow rate of the mainstream as compared to 
the wall jet mass flow rate and both being at the same 
temperature (25°C), the above mentioned temperature dif- 
ference should not cause any appreciable error in measure- 
ments . 

2.2.2 Wall Jet Assembly 

Fig. 2.5 shows the constructional details of 
the wall jet assembly. It is basically a 2-D nozzle 
intended to convert high pressure low velocity fluid to 
a low pressure high velocity fluid. The high velocity 
fluid is then injected into the main flow along the bottom 
wall of the tunnel. 

A specially cut, varnished, poplar wood section 
forms one wall of the nozzle and a 12.7 mm thick aluminum 
divider plate forms the other wall. The flow passage 
between the aluminum divider plate and another poplar wood 
section similar in shape to the first one can be used as wall 
suction system to bleed the upstream boundary layer. 



12. Cylindrical Baffle made of steel sheet 

13. Bottom surface of the tunnel 
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However, this flow passage has been blocked at its top by 
a 13X3X1.6 mm aluminum angle in place as shown in Fig. 

2.5. The exact dimentions of the poplar wood sections are 
shown in Figs. 2.6 and 2.7. The poplar wood was chosen 
as it is less susceptible to moisture absorption and 
warpage . 

A specially extruded aluminum end piece attached 
to the aluminum divider plate by screws from the top 
completes the nozzle shape. The dimensions of the ex- 
truded aluminum end piece are shown in Fig. 2.8. The 
primary purpose of this extruded piece is to turn the 
fluid leaving the bottom portion of the nozzle sc that 
the fluid is injected into the main flow, parallel to 
the tunnel wall. This extrusion also serves as an ex- 
tension of the nozzle as there is some reduction in the 
flow area when the flow passes underneath the extrusion. 
Shims are placed between the aluminum divider plate and 
the extruded end piece wherever necessary to insure a 
uniform opening at the wall jet exit along the length 
of the wall jet. Table 2.1 shows the jet exit gap or the 
"slot height" along the length of the wall jet which has 
a nominal value of 2.72 mm. 

The wall jet assembly is held together by two 
3 mm thick aluminum end plates on each side. The entire 
wall jet assembly is mounted inside a 2.6 mm thick alu- 
minum box section with 114X45 mm outside dimensions and 
with the top face being removed. A threaded hole was 
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TABLE 2.1 

VARIATION OF THE SLOT HEIGHT ALONG 
THE LENGTH OF THE WALL JET 


z 

(nun) 

c 

(nun) 

0.0 

2.72 

12.7 

2.72 

25.4 

2.72 

38.1 

2.72 

50.8 

2 .72 

63.5 

2.72 

76.2 

2.72 

88.9 

2.72 

101.6 

2.72 

114.3 

2.72 

127.0 

2.72 

139 . 7 

2.72 

152.4 

2 . 72 

165 . 1 

2.72 

177.8 

2.72 

190.5 

2.72 

203 . 2 

2 .72 

215.9 

2.71 

228.6 

2.69 

241.3 

2.69 

Z = Distance from the right 

v = slot height 

end (looking downstream) 
of the wall jet 

c 
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cut at the bottom of the aluminum box section on the wall 
jet nozzle side to receive the 25 mm I.D. copper tube 
supplying the air. 

The air from the compressed air supply enters 
the wall jet assembly at its bottom. A baffle plate 
(Fig. 2.9) attached to the aluminum box section diverts 
the flow from the center to the ends , preventing the air 
from being blown directly into the nozzle. The air then 
passes through a 13 mm thick honeycomb with 3 mm cell 
size. After that, the air passes over a 30 square mesh 
steel screen with a wire diameter of 0.17 mm and an open 
area ratio 0.65. The honeycomb and screen act as resist- 
ance to the flow and help to distribute the flow uniformly 
over the entire cross-sectional area of the nozzle 
entrance . 

2.2.3 Performance of the Wall Jet 

The performance of the wall jet is described 
below from the view points of (a) the significant features 
of the present wall jet, and (b) the two-dimensional 
behavior of the flow coming out of the slot. 

2.2.3(a) Significant Features of the Wall Jet 

Fig. 2.10 shows the flow path inside the wall 
jet nozzle in the assembled position with a scale in the 
inset. The present wall jet design is unique in view of 
the few important considerations given below. 
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This wall jet design produces an asymmetric mean 
velocity profile at the exit of the nozzle with a rela- 
tively larger amount of momentum concentrated in the top 
portion of the slot than in the bottom portion. Pig. 2.11 
shows the velocity and turbulence profiles measured with 
a hot wire at x/y c = 0.292. It can be seen here that 
the u t profile is typical of a two-dimensional turbulent 
channel flow. The advantages of having an asymmetric 
mean velocity profile at the nozzle exit are spelled out 
in several places in this work. It is one of the ob- 
jectives of this work to investigate the effect of an 
asymmetric jet velocity profile on the downstream dev- 
elopment of the wall jet flow. 

The asymmetric mean velocity profile at the exit 
of the nozzle is produced in the following way. The flow 
takes a steep 90° turn as it passes from the vertical 
section of the nozzle to the horizontal section. In this 
process the flow experiences a strong pressure gradient 
normal to the streamlines as the streamlines are curved 
around the corner with the static pressure decreasing 
as one moves away from the bottom surface of the nozzle. 

The fluid near the bottom surface of the nozzle exper- 
iences more deceleration than the fluid near the top 
surface of the nozzle because the static pressure is 
greater near the bottom surface than at the top. Hence 
the bottom wall boundary layer is much thicker than the 
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top wall boundary layer, thereby resulting in an asymmet- 
ric velocity profile at the exit. 

Here there is a possiblity that the flow might 
separate near the small radius round corner of the 
bottom wall because of the strong flow curvature. This was 
prevented by gradually decreasing the cross-sectional 
area of the nozzle around the corner by a ratio of 
approximately 2.5:1. The reduction in the flow area 
keeps the flow from separating. To support this fact, 
flow visualization studies have been made around the 
bottom round corner of the nozzle. The bottom wooden 
surface was first coated with a thin layer of wax. 

Water was sprayed in the form of a fine mist around the 
corner and along the length of the nozzle. The jet was 
turned on with operating velocities typical of our 
experiments. If separation occurs around the corner, 
water droplets should have remained there without being 
carried away by the flow. However, all the water drop- 
lets were carried away by the flow, indicating that the 
flow was not separating around the corner. 

Here it should be mentioned that Bowles (1977) 
took measurements of the exit mean velocity and turbulence 
profiles on a wall jet of a similar design but with a 
larger slot exit gap. The asymmetric velocity profile 
typical of the present wall jet design was observed in 
his case also. However the flow cross-sectional area 
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was kept constant around the corner in his wall jet, 
which increases the possiblity of separation near the 
corner with reattachment of the flow downstream. More 
evidence of this behavior is given by the shape of the 
turbulence profile since it has a maximum near the wall 
but well above the bottom wall sublayer, as for a 
reattaching turbulent boundary layer. The details of 
Bowles 'work are included in the Appendix A. 

The disadvantage of the present wall jet design 
is that the frictional losses should be higher as com- 
pared to those for the jet design with a uniform velocity 
profile at the exit because of the steep flow curvature 
necessary for producing the asymmetric velocity profile . 

An attempt has been made to estimate the amount of fric- 
tional losses in the present design. An energy balance 
was made on the control volume surrounding the corner 
of the nozzle, using the velocity profile measurements 
at the slot and the static pressure measurements inside 
the nozzle. 

The estimated frictional losses came to approxi- 
mately 15% of the total energy input. Thus, there is 
a trade off between the increased frictional losses in 
producing an asymmetric velocity profile and the advan- 
tages of having an asymmetric velocity profile. Some 
more practical details have to be worked out before 
exploring the present jet design commercially. Analytical 
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details of the estimation of frictional losses are given 
in Appendix B. 

The thickness of the extruded end piece shown in 
Fig. 2.8 at its downstream lip was kept to a minimum 
value of 0.8 mm within the limits of the commercially 
available manufacturing capabilities. This thin lip 
feature prevents a large region of separation and base 
flow downstream of the lip, which is typical of wall 
jet flows with a thick lip. 

Another one of the useful features of the 
extruded end piece shown in Fig. 2.8 is that the down- 
stream portion of the top surface of this extrusion is 
slanting downwards. This gives rise to an increasing 
flow cross-sectional area as the distance between the 
jet body and the top wall of the tunnel increases as 
one goes downstream. In effect, the flow on this 
slanting surface experiences deceleration, resulting in 
a larger velocity gradient between the jet and the up- 
stream boundary layer flow at the lip, resulting in 
better mixing. There is a possibility of flow separation 
around the crest of the extruded end piece because of the 
sharp edge there. This was checked by making flow visual- 
ization studies around the crest. No observable separa- 


tion of the flow was found there. 
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2.2.3 (b) Two-Dimensionality of the Wall Jet Flow 

The uniformity 01 the dynamic pressure along the 
wall jet exit was used as a check of the two-dimensional 
behavior. The exit dynamic pressure was measured with 
a square brass tube, the dimensions of which are shown 
in Fig. 2.12. The brass tube was kept with its cross 
section perpendicular to the flow and flush with the 
bottom wall of the tunnel. It is connected by a flex- 
ible tubing to an inclined tube manometer. Thus the 
measured dynamic pressure is the average dynamic pressure 
over the opening of the square brass tube. 

The top wall of the wind tunnel was removed 
during the wall jet flow measurements, thereby ensuring 
that the static pressure at the exit is atmospheric. The 
wall jet flow measurements were made with no external 
stream. The dynamic pressure at the exit of the wall 
jet was measured at 25 mm intervals along the length of 
the jet starting from the center. The average exit 
velocity as measured by the square brass tube was set 
approximately at 42.7 m/sec. A second set of measurements 
was also taken at 36.6 m/sec. To start with, there was 
a 5% variation in the average exit velocity at various 
stations along the length of the wall jet. Particularly, 
there was a deficiency of flow at about 38 mm on each side 
from the center line of the tunnel indicated by the low 
value of the average velocity. To correct this, three 
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holes (Fig. 2.9) of 4 mm diameter were drilled into the 
baffle plate at the entrance of the wall jet plenum. 

These holes allow more air to be supplied to the fluid 
deficient area, thus equalizing the average velocity at 
several locations. However, this step alone did not solve 
the problem completely. There was still a significant 
variation in the average dynamic head, although it was 
less than what it was without the above mentioned holes. 

Another modification was made in the wall jet 
plenum to improve flow uniformity. A 0.1 mm thick steel 
sheet was rolled into a cylinder and inserted into the 
copper feed tube with a small portion of the steel sheet pro- 
jecting out as shown in Fig, 2.9. The height of the 
projecting portion of the steel sheet war adjusted to 
produce the least variation in the average value of the 
dynamic head measured by the square brass tube. The steel 
sneet was then glued to the _opper tube along its edges . 

The steel sheet acts as a baffle to divert the flow to 
the fluid deficient area. The projecting portion of 
the steel sheet on the left half (looking upstream) is 
approximately 1.6 mm and that on the right half (looking 
upstream) is approximately 0.8 mm. 

With the above modifications, the wall jet main- 
tained a uniform flow indicated by the average velocity. 

The average velocity as measured by the square brass tube 
varied within ±0.9% of its value at the center for the 
entire length of wall jet. There was a small region of 
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high velocity at 59 mm from the center to the left (look- 
ing upstream) where the average velocity was higher than 
the central value by 1.4 55-.. 

2.2,3 (c) Calibration of the Wall Jet Flow 

The static pressure inside the wall jet nozzle is 
related to the amount of the flow at the exit of the wall 
jet. Therefore the wall jet exit flow rate can be cali- 
brated against the static pressure measurements inside 
the nozzle. A static pressure tap was made in the end 
plate on the left side (looking downstream) of the 
wall jet nozzle. The position of the static pressure 
tap is given in Fig. 2.5. The static pressure at this 
location was measured by an inclined tube manometer. 

The wall jet was calibrated for different static 
pressure readings by measuring the velocity profile at 
the center of the wall jet. The velocity profile measure- 
ments were made with a rectangular mouth boundary layer 
total head probe shown in Fig. 3.1. The measurements were 
made after removing the top wall of the wind tunnel and 
without the external stream. The static pressure inside 
the wall jet was altered by adjusting the pressure regulator 
in the air input line. The nominal values of the maximum 
velocities in the velocity profile at the center of the jet 
were 30.5, 36.6, and 42.7 m/sec, correpsonding to the three 
different static pressure settings for which the measurements 
were made. Thus for a given static pressure inside the 
wall jet, the flow rate can be calculated from the measured 
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velocity profiles at the jet exit. By means of this method, 
the wall jet nozzle can also be used as a flow meter. 

As a further check on the two-dimensional behavior 
of the wall jet, velocity profiles were also taken at 
2.5 cm and 7.6 cm on each side of the center and for the 
above mentioned static pressure settings. The maximum 
variation in velocity at any given height from the floor 
of the tunnel and at any given static pressure setting 
was within 3% . In all the above mentioned velocity 
profile measurements the asymmetric nature of the jet exit 
velocity profile was observed. 


CHAPTER III 


INSTRUMENTATION 

This chapter describes the different instruments that 
were used for measuring the mean and turbulent components 
of velocity, the Reynolds shear stress, and the skin 
friction. The instrumentation includes the pitot probes, 
Preston tubes, hot-wire probes, and the hot-wire signal 
processing equipment. A brief description of the calibra- 
tor is also included. The estimated uncertainties in 
the measurements using the above instrumentation are pre- 
sented . 

3.1 Pitot Probes and Preston Tubes 

Pitot tube traverses were made only to check the 
two-dimensionality of the flow. The pitot tube used for 
this purpose is shown in Fig. 3.1. It is a total head 
rectangular-mouth boundary layer probe. The same pitot 
tube was also used 

1. For measuring the jet exit velocity profile 
during the calibration of the wall jet, and 

2. For the upstream boundary layer velocity 
profile measurements necessary for selection of the 
roughness element, as described in Chapter 2. 
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All other static and total pressure measurements were made 
with a United Sensor and Control Corporation model 
PAA-8-KL pitot tube. It had proper connections for mea- 
suring both the total head and the static head. 

The skin friction was measured by means of a 
Preston tube using Patel's (1965) calibration curve. 

The Preston tube used in the zero pressure gradient 
studies has the mouth dimensions of 0.5 mm I.D. and 
0.6 mm O.D. The Preston tube used ir e reverse pressure 
gradient studies has the mouth dimer > xw: of 0.78 mm I.D. 

and 0.88 mm O.D. 

3.2 Hot-wire Probes 

A TSI model 1274-10 normal hot-film probe was 
used to measure the mean velocity U and the turbulence 
intensity u./U • A TSI model 1273-T 1.5 slant hot- 
wire was used to measure the Reynolds shear stress -uv 
and the turbulence intensities v./U and w, /U . The 
45° slant wire probe is rotatable about its axis through 
a rotating mechanism described below. The shear stress 
and the normal stresses were deduced from the slant-wire 
measurements at seven angular positions about its axis. 

The choice of a single rotatable slant-wire over X wires 
was made for three reasons: 

1. Any uncertainty about possible interaction 
between two wires and four prongs is eliminated. 
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2. The need for two anemometers and two linearizers 
together with a critical matching of the two hot-wire 
calibrations is avoided 

3. w can also be measured with the rotatable 
single slant wire. 

A TSI model 1210-20 rake hot wire probe was used 
to measure the velocity in the freestream at several 
locations which was required to adjust the pressure grad- 
ient in the case of zero pressure gradient flow. 

3.2.1 Slant-wire Rotating Mechanism 

Fig. 3.2 shows the slant-wire rotating mechanism 
with the probe mounted in it. The probe is free to rotate 
inside a cylindrical steel guide with its face turned 
smooth on a lathe to reduce the flow blockage. The steel 
guide is firmly attached to a rectangular rod which in 
turn is attached to another round rod in the traversing 
mechanism through a swivel joint. This swivel joint 
helps in rotating the slant wire probe about a horizontal 
axis through the joint. The probe can also be rotated 
about a vertical axis by turning the round rod in the 
traversing mechanism. These two movements are needed for 
the alignment of the probe with the free-streamlines . The 
probe stem itself passes through another slotted spindle 
with 12 slots cut on its rim at 30° intervals and can be 
attached to the probe stem by a set screw. The probe 
stem then passes through a rectangular steel block with 
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a passage hole for the probe stem. This steel block is 
held in position by means of the rectangular rod mentioned 
above. A small steel strip that fits exactly into the 
slots of the slotted soindle is hinged to the steel block 
and can be raised or lowered into the slot by means of 
a fishing line attached to it and brought outside the 
tunnel through one of the probe holes. By raising the 
steel strip, the probe is free to rotate about its axis 
and by lowering the steel strip into a slot, the probe 
position can be locked at any one of the twelve orientations . 
Two small steel stops, one on each side of the steel 
block, prevent any horizontal movement of the probe. The 
first stop is attached to the rectangular rod by a 
spring and the second stop is attached to the probe stem 
by a set screw. The probe is rotated from the end of 
the tunnel by the cable carrying the hot-wire signal. 

With the above mechanism for rotating the slant wire, 
it was possible to rotate the probe with only an eccentri- 
city of 0.25 mm between the probe stem axis and the axis 
of the mechanism. 

3.2.2 Traversing Mechanism 

The traversing mechanism used for the traversal 
of different probes consists of a brass screw having a 
movement of 0.635 mm per turn along with proper mountings. 

The probe can be set to an accuracy of .0254 mm by means 
of this traversing mechanism. The traversing mechanism 
is shown in Fig. 3.3. 
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a passage hole for the probe stem. This steel block is 
held in position by means of the rectangular rod mentioned 
above. A small steel strip that fits exactly into the 
slots of the slotted spindle is hinged to the steel block 
and can be raised or lowered into the slot by means of 
a fishing line attached to it and brought outside the 
tunnel through one of the probe holes. By raising the 
steel strip, the probe is free to rotate about its axis 
and by lowering the steel strip into a slot, the probe 
position can be locked at any one of the twelve orientations. 
Two small steel stops, one on each side of the steel 
block, prevent any horizontal movement of the probe. The 
first stop is attached to the rectangular rod by a 
spring and the second stop is attached to the probe stem 
by a set screw. The probe is rotated from the end of 
the tunnel by the cable carrying the hot-wire signal. 

With the above mechanism for rotating the slant wire, 
it was possible to rotate the probe with only an eccentri- 
city of 0.25 mm between the probe stem axis and the axis 
of the mechanism. 

3.2.2 Traversing Mechanism 

The traversing mechanism used for the traversal 
of different probes consists of a brass screw having a 
movement of 0.635 mm per turn along with proper mountings. 

The probe can be set to an accuracy of .0254 mm by means 
of this traversing mechanism. The traversing mechanism 
is shown in Pig. 3.3. 
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3.3 Electronic Signal Processing 

A TSI model constant temperature anemometer control 
unit and a TSI 1055 linearizer were used to produce the 
hot-wire signals. The root-mean-sguare voltages were 
measured using a HP model 400E voltmeter connected in 
series to a true integrating voltmeter, consisting of 
a Tektronix model DC503 counter. The signal was averaged 
over a lu second interval. Two readings were taken for 
each data point and then an average was taken giving an 
effective averaging time interval of 20 seconds. The mean 
voltages were measured with a true integrating voltmeter, 
consisting of a Wavetek model 131 function generator 
with a voltage controlled frequency circuit, connected 
in series to an Anadex model CF600 counter. The mean 
voltages were also averaged over 10 second intervals . Two 
readings were taken for each data point and then an average 
was taken giving an effective averaging time interval of 
20 seconds . The spectra were measured with a Princeton 
Applied Research Inc. model 4512 FFT real time spectrum 
analyzer. This spectrum analyzer was used over a frequency 
analysis range of 10Hz to 40KHz and a sensitivity range of 0.1 
to 10 volts. The output of the spectrum analyzer was 
recorded on a Honeywell 320 XY recorder. 

3 . 4 Calibrations 

The calibration of different probes was done in 
a TSI model 1125 calibrator which can supply nearly 
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turbulence free air at a given temperature. The calibrator 
is connected to a filtered compressed air supply available 
in the laboratory. The temperature of the air coming 
out of the calibrator can be electrically controlled by 
a heat exchanger in the flow path. The calibrator has 
different sets of nozzles that can be fitted in position 
for different velocity ranges and also it has two inner 
chambers for low velocity calibration of the probes. 

The linearizer was adjusted to give a linear 
calibration within ±1% between the velocity and the output 
voltage of the linearizer. The linear calibration curve 
was obtained by a least square fit of the calibration 
points. During most of the experimental runs, the anemo- 
meter remained drift free. This was indicated by the 
calibrations done before and after each experimental 
run. The anemometer circuit had a flat frequency response 
in the range of frequencies (0-10 KHz) encountered in 
the present flows. 

3.5 Uncertainties in the Measurements 

The method of Kline and McClintock (1953) was 
used to estimate the uncertainties in the calculated data 
obtained from the primary measurements . The uncertainty 
of a particular variable R is denoted by AR. The uncertainty 
figures given in table 3.1 are the maximum possible values 
calculated for the cases of maximum possible error. 
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TABLE 3.1 

UNCERTAINTIES IN THE DATA 



Measurement 

Uncertainty 

Device Used 


U 

AU = ±0.4 m/sec 

Pitot tube 

U 

AU = ±0.15m/sec 

Normal hot-film 
and hot-wire 

C f 

AC f = ±.00033 

Preston tube 

~2 

u 

A(i?)/u 2 = ±3% 

Normal hot-film 

-uv 

A(-uv)/-uv = ±10% 

Slant hot-wire 

T 

V 

A(v^/v^" = ±12% 

Slant hot-wire 

w 

A(w 2 )/w 2 = ±12% 

Slant hot-wire 


The uncertainties in the mean velocity measure- 


ments with pitot tube and the skin friction measurements 

with Preston tube are mainly due to the uncertainties 

in the pressure measurements. The uncertainty in the probe 

calibration is the main source of uncertainty in the 

~2 

measurements of U and u with the normal hot-film. The 
~ ~7 

data of -uv, v , and w are mainly affected by the un- 
certainty in the probe calibration, uncertainty in the 

determination of the constant (Appendix C) and the 

~~2 

uncertainty in the measurement of u . 


CHAPTER IV 


EXPERIMENTAL RESULTS FOR THE ZERO 
PRESSURE GRADIENT FLOW 

The experimental results for the zero pressure 
gradient flow are presented in this chapter. The objec- 
tives of the experiments are: 

1 . To obtain experimental data with thick 
upstream boundary layers and with low jet velocity ratios, 

2 . To obtain detailed turbulence data which would 
aid in the future developemnt of turbulence models and 

3. To observe the effect of an asymmetric jet 
velocity profile on the flow development downstream of 
the slot. 

The measured quantities presented here are the 
mean velocity, turbulence intensities, shear stress, 
skin friction and spectra. The derived quantities pre- 
sented include the integral and profile parameters, eddy 
viscosity, mixing length, Prandtl-Kolmorgor gth 

scale, turbulent kinetic energy, correlation iicients, 

production and dissipation rates of turbulent tic 

energy and bursting periods. A brief description of the 
flow conditions is given first followed by the presenta- 
tion of the measured and derived experimental results . 
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4 . 1 Flow Conditions 

The method of setting the zero pressure gradient 
flow, the two-dimensional nature of the flow and the 
qualities of the upstream boundary layer are described 
below in sequence. 

The flow was set for a zero pressure gradient 
condition by measuring the free-stream velocity at several 
stations beginning from the slot by means of a rake hot- 
wire probe. The top wall of the tunnel was adjusted to 
give the same free-stream velocity at different stations 
along the tunnel. The initial adjustments of the pressure 
gradient were done by measuring the free-stream velocity 
with a pitot tube and the final setting was done with a 
rake hot-wire probe. The variation of the free-stream 
velocity measured with the hot-wire probe is shown in 
Table 4.1. The free-stream velocity was constant from the 
slot to the exit of the tunnel with an average variation 
of ±1%. The measured free-stream turbulence intensity 
u./U^ was found to be about 0.2%. The height of the top 
wall above the bottom wall at different locations is shown 
in Fig. 4.1. The wall jet was adjusted to have a maximum 
velocity in the slot velocity profile approximately 
equal to 37.8 m/sec for the measurements under - zero 
pressure gradient. The flow rate through the jet slot 
was kept cons L.ant during the measurements by maintaining 
a constant static pressure difference between the inside 
of the jet nozzle and the free-stream above the jet exit. 
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1.000 


35.05 

1.001 

53.74 

0.998 

74.48 

0.997 

92.29 

0.996 

108.36 

0.996 

146.18 

1.001 

182.83 

1.000 
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1.000 

260.80 

0.998 

289.72 
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The important consideration which affects the 
quality of the flow in the present situa ion is the two- 
dimensional nature of the mean flow. The two-dimensional 
behavior of the flow was checked at several points : 

1. Upstream of the slot at approximately 45 cm 
from the contraction exit 

2. Downstream of the slot at approximately 
5 cm from the slot exit, and 

3 . Downstream of the slot at approximately 61 
cm from the slot exit. 

These checks were made with the pitot tube traversals after 
setting the pressure gradient and the wall jet velocity 
to proper values. The velocity profiles were measured 
at the center of the tunnel and at 7.6 cm from the center 
on each side at each of the locations mentioned above. 

The velocity at any given y location on either side of 
the center was found to be within a maximum of ±4% and 
within an average of ±1.5% of the velocity at the center 
of the tunnel. A further check on the two-dimensional 
nature of the flow was made in section (4.2.3) by examin- 
ing the two-dimensional integral momentum equation with 
the present measurements. 

One of the objectives of the present measurements 
is to obtain data with a large upstream boundary layer at 
the slot. Therefore, the nature of the upstream boundary 
layer was studied by making velocity profile measurements 
with a normal hot-film at 21 cm upstream of the slot. 

Fig. 4.2 shows the velocity profile at 21 cm upstream of the 
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slot plotted on the conventional semilogarithmic co- 
ordinates . The skin friction was determined from the 
Preston tube measurements. The velocity profile measure- 
ments are compared with the logarithmic law of the wall 
represented by the solid line in Fig. 4.2. The constants 
recommended by Patel (1965) were used in the logarithmic 
law of the wall given below: 

= 5.5 Log 10 -J- + 5.45 (4.1) 

T 

The velocity profile data agree well with the logarithmic 
law of the wall at the upstream station mentioned above. 

The momentum thickness Reynolds number is equal to 3142 
at this station. The strength of the wake component A(U/U ) 
(Coles, 1962) is equal to 2.2. This value of A(U/U ) agrees 
within 10% of its value given by Coles, for the same momentum 
thickness Reynolds number and for a "normal" turbulent bound- 
ary layer under a zero pressure gradient. Therefore, the 
boundary layer upstream of the slot is a developed turbulent 
boundary layer. 

4.2 Mean Flow Data 

The mean flow data presented here include the 
skin friction, the mean velocity profiles and the integral 
and profile parameters. 

4.2.1 Skin Friction C^ 

The variation of skin friction coefficient C^ 
obtained from the Preston tube measurements is shown in 
Fig. 4.3 (a). Fig. 4.3 (b) shows the variation of C f 
with x/y c on a logarithmic scale. It can be seen that 




66 


there is a power law variation between skin friction 
coefficient C f and x/y c . 

It has been found in the literature that for self 
preserving wall jet flows with and without pressure 
gradients, the skin friction varies according to the 
following relation: 


w 


^pu; 


= A 


max 

where A and B are constants. 


U Y 
max ma x 

v 


B 


(4.2) 


This relation holds only for self-preserving wall 
jet flows with only a velocity maximum and with a negli- 
gible upstream boundary layer. Irwin (1974) suggested 
on the basis of the recommendations made by Guitton (1970) , 
Patel (1962) , Bradshaw and Gee (1962) , Kruka and Eskinazi 
(1964) , and McGahan (1965) that A = 0.026 and B = -0.18 in 

the equation (4.2). Fig. 4.3 (c) shows the variation of 
2 

C- = I,/ (%pU ) in the present case in comparison with 
rm w max 

the above equation for C- . The present data also show a 
power law variation of . However, the constants A and 
B are different than those suggested by Irwin. For the 
present data, A - 0.102 and B = -0.33. The difference 
in the constants A and B may be due to the non-self- 
preserving nature of the present flow. 

4.2.2 Mean Velocity U 

Figs. 4.4 (a-e) show the normal hot-film measurements 
of mean velocity profiles in the U/U vs y/6 coordinates. 


The negative sign on the x/y c values in Fig. 4.4 (a) 
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indicates that they are the distances measured upstream 
of the slot. The value x/y c = -16.4 corresponds to the 
crest of the extruded aluminum end piece (section 2.2.2) 
and x/y = -77.1 corresponds to a position 21 cm upstream 
of the slot. 

Fig. 4.4(b) shows the velocity profile at x/y c = 

0. 292. Here it can be seen that the velocity profile 
in the jet is asymmetric with a relatively greater con- 
centration of momentum in the upper half of the wall jet. 
This asymmetric velocity profile in the jet is typical 

of the present wall jet design as discussed in Chapter 
2. By comparing the upstream boundary layer velocity 
profiles at x/y = -77.1 and x/y = 0.292, one can see 
that the wall jet body did not introduce any major changes 
in the mean velocity profile of the upstream boundary layer. 
The upstream boundary layer above the slot shown in Fig. 
4.2(b) is very thick and it has a large deficit of momentum, 
satisfying one of the conditions under which the present 
measurements were intended to be made, as given in Chapter 

1. Figs. 4.4 (c) / (d) and (e) show the mean velocity 

profiles from x/y c = 7.45 and onwards. The upstream 
boundary layer has been completely absorbed by the jet 
as the flow proceeds in the downstream direction. The 
velocity maxima and minima cannot be identified after 
x/y c = 146.2. At x/y c = 435.2, the velocity profile 
looks similar to a normal turbulent boundary layer. 
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Figs. 4.5(a) and (b) show the velocity profiles 

plotted on a conventional semi logarithmic plot with 

U/U T and U^y/y as the co-ordinates . The velocity profiles 

were shown for stations x/y = 7.45 and onwards. The fric- 

- 1 c 

tional velocity U^. was obtained from the Preston tube 
measurements . The velocity profile measurements were 
compared with the logarithmic law of the wall (Equation 
4.1), with the constant recommended by Patel (1965) . It 
can be seen from Figs. 4.5 (a) and (b) that the experi- 
mental data agree well with the logarithmic law of the 
wall for all the stations, except at x/y c = 7.45 where a 
defined logarithmic region has not been formed yet because 
of its proximity to the slot. 

The good agreement between the logarithmic law 
of the wall and the mean velocity data indicates that the 
skin friction measurements are accurate. This also 
indicates that the uncertainty in the skin friction 
measurements given in Chapter 3 may have been over 
estimated. The momentum thickness Reynolds number and 
the wake component A (U/U ) (Coles, 1962) are equal to 
6312 and 3.0, respectively, at x/y c = 435.2. The wake 
component at x/y c = 435.2 is 10% higher than its value 
given by Coles for a normal zero pressure gradient tur- 
bulent boundary layer at the same or higher momentum 
thickness Reynolds numbers. The mean velocity data 
are tabulated in Appendix D along with the data of u^. 
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Fig. 4.5 (c) shows the streamline pattern around 
the slot lip. It shows that there is no appreciable 
variation of the static pressure in the direction normal 
to the bottom wall of the tunnel as the streamlines are 
not steeply curved in the region of the slot lip. 

4.2.3 Integral Parameters 

The integral parameters evaluated are the boundary 
layer thickness (6), displacement thickness (6-^) , momentum 
thickness (6 2 ), s ^ a P e factor (H) , and the momentum thick- 
ness Reynolds number (Re 2 ) . Figs. 4.6 (a-e) show the 
development of the integral parameters 6, 6^, 6 2 < H, and 
Re 2 respectively. The rates of increase of the momentum 
thickness, displacement thickness, and the boundary 
layer thickness seem to be very slow in the present ex- 
periments. The shape factor H tends to become a constant 
approximately equal to 1.35 at far downstream stations, 
which closely agrees with the value of H = 1.4 indicated 
by Schubauer and Klebanoff (1955, 1956) for a flat 
plate turbulent boundary layer. 

The two-dimensional nature of the flow was exam- 
ined by applying the two-dimensional integral momentum 
equation (4.3) to the present data. 

dRe 2 
dRe x 
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dx 
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Fig, 4.6(b) Growth of the Displacement Thickness 
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Fig. 4.6(c) Growth of the Momentum Thickness 
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Fig. 4.6(d) Variation of the Shape Factor 
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Fig. 4.6(e) Variation of the Momentum Thicknes 
Number Re^ 


H5D.0D 


B 


S5Q.0D 


Reynolds 


Re 2 obtained using the two-dimensional 
integral momentum equation (4.3) 
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or, 


x/y. 


x/y. 


Re 2 


— dRe 


X 


x/y c =16.5 


+ (Re 2 ) 


x/y =16.5 

w 


"T T 

u -v 


o D - 


dy 


x/v^=16 . 5 


where 


U fi, U x 

Re 0 = and Re = 

v x v 

f ~7 ~7 2 

The contribution of the normal stresses term J[ (u - v )/U ro ]dy 
to the integral momentum equation was found to be smaller 
than 5% and hence it was neglected. The measured skin 
friction coefficients were used in the equation (4.3). 

Pig. 4. 6 ® shows the values of Re 2 obtained from the two- 
dimensional integral momentum equation (4.3) along with 
the experimental values of Re 2 - The experimental values 
of Re 2 agree with the values of Re 2 obtained from equation 
(4.3) within ±10% on the average, establishing the two- 
dimensionality of the flow. 

4,2.4 Effect of the Asymmetric Jet Velocity 
Profile on the Mean Flow 

The effect of the asymmetric jet velocity profile 
on the mean flow development is discussed below. A crude 
comparison of the development of the mean velocity profile 
for the present flow with that of similar flows in the 
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literature (Kacker and Whitelaw, 1968) shows that the 
development of the jet requires a greater x/y distance 
in the present case. The flow takes a longer distance in 
the present case to achieve the state of a normal turbulent 
boundary layer. The effect of the jet momentum can be 
seen even as far as x/y c = 435.2 due to the low value of 
Re 2 there. This can be attributed only to the asymmetric 
nature of the jet velocity profile. The asymmetric 
velocity profile spreads the jet momentum more uniformly 
in the layer by supplying more momentum to the momentum 
deficient upstream boundary layer than wasting the same 
jet momentum in the form of friction at the wall. In 
comparison, the uniform jet velocity profile has a 
relatively large concentration of momentum near the wall 
which results in greater frictional losses. Thus uhe 
momentum of the jet is carried away to a much longer 
distance in the case of an asymmetric jet velocity 
profile . 

The rates of increase of momentum thickness, 
displacement thickness and boundary layer thickness seem 
to be very slow for the present case where a highly 
momentum deficient upstream boundary layer is meeting a 
wall jet of moderate momentum. In a case like this, 
one would expect the integral thicknesses to increase 
rapidly, indicating deficiency of momentum if it was a 
uniform profile in the jet. The uniform jet velocity 
profile cannot meet the momentum requirements of the 


ar v 


o ,< 


upstream boundary layer as efficiently as an asymmetric 
jet velocity profile. This results in the development of 
momentum deficient regions in the case of a uniform jet 
velocity profile. Hence the integral thicknesses can be 
expected to increase rapidly in the case of uniform jet 
velocity profile. Therefore, the present slow growth rate 
of the integral thicknesses can only be attribututed to the 
asymmetric nature of the jet velocity profile. 

4.2.5 Profile Parameters 

The profile parameters are the quantities related 
to the mean velocity profile. The profile parameters pre- 
sented here include Y . . Y, tr „ and U_ ,• _ . With 

max' mm half max min 

reference to Fig. 1.1 (b) , Y is the position of the max- 

ITlaX 

imum velocity U max , Y m ^ n is the position of minimum velocity 


U . , and Y 


half 


is the position .where U = (U + U._- )/2. 


The velocity maximum and the velocity minimum could not be 
identified after x/y c = 146.2. Hence, all the profile para- 
meters were plotted only up to x/y c = 146.2. 

4.2.5 (a) Development of Y max , Y half , and Y^ 

Fig. 4.7 (a) shows the development of Y m _ v and 

lUci X 

Yhaif* Fig. 4.7 (b) shows the development of Y m ^ n . The 
developments of Y max and Y^ al ^ are of particular interest 
here. It has been found in previous research that Y max 
can be expressed as a universal linear function of the 
distance x/y c for self-preserving wall jet flows under 


zero and adverse pressure gradients. The survey of the 
existing data on self-preserving wall jets by Narayan 
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and Narasimha (1975) confirms the universal behavior of 
the Y distribution. The data ri in (19 73) also 

Ula.X 

reveal the universal behavior of the Y distribution. 

max 

The universal Y max distribution for the self-preserving 

wall jet flows is shown in Fig. 4.7 (a). However, it 

is not clear whether the non-self-preserving wall jets 

exhibit similar universal behavior for the Y _ develop- 

ment. The data of Ramaprian (1973) show that even for 

non-self-preserving wall jets under adverse pressure 

gradients, the Y max development is universal. The data 

of Irwin (Newman and Irwin, 1975) do not indicate a 

universal behavior for the Y „ development in the case 

max 

of non-self-preser /ing adverse pressure gradient flows. 

The present flow is a non-self-preserving flow. 

The Y distribution for the present flow does not seem 
to follow the universal distribution as revealed by the 
Fig. 4.7 (a). The growth rate of Y max in the present 
case is higher than that of the universal Y max distribution. 
This can be partly attributed to the asymmetric jet 
velocity profile as explained below. The experimental 
asymmetric jet velocity profile has its maximum velocity 
located nearer to the slot lip to start with. Hence, 
the asymmetric velocity profile has relatively more 
momentum near the slot lip. The inner layer correspond- 
ing to the region of maximum velocity spreads out more 
rapidly in an effort to meet the momentum requirements 
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of the upstream boundary layer. Therefore, the position 
of maximum velocity (Y max ) moves more rapidly into the 
outer layer. In comparison, a wall jet with uniform jet 
velocity profile has its momentum distributed evenly 
across the slot. Hence, the growth rate of Y m „„ is 
relatively slow for a uniform jet velocity profile. 

Also, a wall jet with uniform jet velocity profile retains 
a considerable part of injected momentum nearer to the 
wall, thereby losing a relatively larger amount of momen- 
tum as surface frictional losses. 

The data of Irwin (1973), Irwin (Newman and 
Irwin, 1975), Ramaprian (1973), Gartshore and Newman 
(1969) , and the data of various authors presented in 
Narayan ' s (1973) work indicated that the growth rate of 
Y half a ^ wa Y s higher than that of Y max for wall jets 
with uniform injection at the slot. The above data 
also indicate that the development of Y^aif is not 
universal d it depends strongly on the conditions 
upstream of the slot and the ratio of the jet velocity 
to the free-stream velocity. However, the growth rate of 
Y half present experiments is almost the same as that 

of Y max shown in Fig. 4.7 (a). The slower growth rate of 
Y half ; '' n t * ie present experiments can be attributed to the 
asymmetric jet velocity profile. The experimental 
asymmetric profile has relatively higher concentration of 
momentum in the upper half of the jet and nearer to the 
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momentum deficient region in the outer layer. Therefore, 

the outer layer tries to extract more momentum from the 

inner layer rather than extracting momentum from the free- 

stream. This results in slower growth of the outer layer 

and hence the slower growth rate of 

4.2.5 (b) Variation of If and U . 

max min 

Fig. 4.8 (a) shows the variat ' n of u max / u j ave and 
Umi n / U jave with X /Y c * Pi 9‘ 4 - 8 shows the variation 

of (D ma* - D min )/u j aTC " ith x/ >V U Jave is the ■ niforI " 
jet velocity for an equivalent jet with uniform profile 

and having the same momentum as the experimental asymmetric 
jet. Wall jet data in the literature show that the varia- 
tion of U max and U m ^ n depends on several parameters like 
the ratio of jet velocity to the free-stream velocity, 
the pressure gradient and the conditions upstream of the 
slot. Hence an attempt has not been made to compare the 
present variation of U and U • with the data in the 
literature. It can be seen from Fig. 4.8 (a) that the 
rate of decay of U max is greater than the growth rate of 
U min* Ramaprian’s (1973) data for different non-self- 
preserving wall jets under adverse pressure gradients 
show that for a given ratio of jet velocity to free- 
stream velocity, the decay of maximum velocity is the 
same as that for plane self-preserving wall jets under 
zero pressure gradient. However, this is not true in 
general . 


(Umax“Umi 
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4 . 3 Turbulence Data 

The measured turbulence data include the results 
of the measured turbulence intensities (u^/U^, v t /U ro and 
w. /U ) , the measured turbulent shear stress ~uv and the 

U 00 

T 

measured spectra of u . The quantities derived from the 
measured turbulence data include the eddy viscosity, 
mixing length, Prandtl-Kolmogorov length scale, turbulent 
kinetic energy, correlation coefficients, production and 
dissipation rates of turbulent kinetic energy and bursting 
periods. The results of the measured and derived turbulent 
quantities are given below in detail. 

4.3.1 Turbulence Intensity u t /U ro 

The u t data were obtained from normal hot- film 
traversals at several stations starting from x/y = 0.292. 
Figs. 4.9 (a-c) show the variation of u t /U ro vs y/^ for 
different x locations. The negative sign on the value of 
x/y e indicates that those stations are located upstream 
of the slot. The positions of the stations x/y = -16.4 
and x/y = -77.1 are as given in section 4.2.2. The pro- 
file of u./U at x/y = 0.292 shows that the flow inside 
the jet is not an inviscid one, instead it is a fully 
turbulent flow, because of the high turbulence intensities 
exist!.. c thei-e. The level of turbulence intensity grad- 
ually decreases as one proceeds downstream starting from 
the j et - 

Fig. 4.10 (a) and (b) show the variation of u t 
in wall coordinates. Most of the data are limited to 
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Fig. 4.9(c) Distribution of the Turbulence Intensity u /U at 
x/y c = 146.2 to 435.2 t “ 

Uj./U^ data by Klebanoff (1955) for a turbulent 

boundary layer on a flat plate (reproduced from 
Rotta (1962)) 
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the outer edge of the viscous sublayer with few stations 
as exception. The point of first maximum in u fc was 
always found to be near the wall and at the edge of the 
viscous sublayer. The point of minimum u fc in the inner 
layer generally corresponded with the region of maximum 
velocity even though the exact position was generally 
nearer to the wall than the point of maximum velocity. 
Sufficient slant-wire data were not available in the inner 
layer to describe the location of mimimum u t with respect 
to the point of zero shear stress. The point of second 
maximum in u t was found to be in between the points of 
velocity maximum and the velocity minimum and nearer to 
the maximum velocity gradient. The u fc was found to be 
either minimum or constant around the point of velocity 
minimum. The steep drop in u fc far away in the outer 
layer was found to begin from the point of maximum velocity 
gradient between the velocity minima and the freestream. 

The data for u^. are given in Appendix D. 

4.3.2 Turbulence Intensities v^/U^ , w t /U 00 

And Shear Stress -uv 

A rotatable slant-wire (Chapter 3) was used to 



obtain v^ , w^, and the shear stress -uv. The u data 

obtained from the normal hot-film were used in the solution 

~2 ~ 

of the simultaneous equations for -uv, v and w . The 
details of obtaining v^, w^, and -uv from the slant wire 
data including the relevant mathematical details are given 
in Appendix C. The data of v^_ , w t , and -uv for all the 
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slant-wire data stations are tabulated in Appendix E. 
Appendix E also shows the interpolated data of u t and 
3 (U/U^) /3 (y/6) obtained from the normal hot-film data. 
Because of the size limitations set by the slant-wire 
probe supports, the slant-wire data at almost all the 
stations were taken only above the point of velocity 
maxima with a few exceptions. The results of the tur-r 
bulence intensities v.^/U^, w t /U ro , and the shear stress 
-uv are described below in detail. 

4.3.2 (a) Turbulence Intensities v./U and w. /U 

t a> C oo 

Pigs. 4.11 (a), (b) and 4.12 (a), (b) show the 
variation of v.j./U^ and w^/U^ respectively at several 
stations, starting from x/y = 16.5. The slant wire 

W 

data at x/y = 0.292 and x/y = -16.4 were not shown on 

G C 

the plots . 

The point of maximum in v^ and w^ between the 
points of velocity maximum and velocity minimum has been 
found to be in the region of maximum velocity gradient. 

The level of magnitude of v t and w^ dropped very rapidly 
between x/y = 16.5 and x/y = 74.5 and then the level 
of magnitude remains the same further downstream. For 
stations up to x/y c = 108.4, the magnitudes of both v t 
and w t remained nearly equal at a given station. However, 
for stations beyond x/y c = 108.4, w fc was found to be 
higher than v. indicating the tendency to become a normal 
boundary layer (Klebanoff , 1955) . 
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Fig, 4.11(a) Distribution of the Turbuli 

v /U at x/y =16.5 to 74 
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Fig. 4.11(b) Distribution of the Turbulence Intensity 
v t /U co at x /y c = 108.4 to 435.2 

v t /U oo data by Kiebanoff (1955) for a 

turbulent boundary layer on a flat plate 
(reproduced from Rotta (1962)) 
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Fig. 4.12(b) Distribution of the Turbulence Intensity 
w./U at x/y = 108.4 to 435.2 

t CO ■> c 

w /U data by Kiebanoff (1955) for a 

L 00 

turbulent boundary layer on a flat plate 
(reproduced from Rotta (1962)) 


4.3.2 (b) Shear Stress -uv 


Figs. 4.13 (a) and (b) show the distribution of 

2 

-uv/U at different stations. The slant-wire data at 
00 

x/y_ = 0.292 and x/y_ = -16.4 were not shown on the plots. 

W W 

The points of zero shear can be seen distinctly in the 

-uv distributions shown in Figs. 4.13 (a) and (b) . At 

a given station, there should be two points of zero 

shear corresponding to the two points of zero velocity 

gradients. The two points of zero velocity gradients 

refer to the points of maximum and minimum velocities 

in the velocity profile. The first point of zero shear 

encountered as one goes away from the bottom wall is called 

the inner point of zero shear and the second zero shear 

point is called the outer point of zero shear. 

In the present -uv data, both the inner ana outer 

points of zero shear can be observed only at stations 

x/y = 108.4 and x/y = 146.2. Only the outer point of 
c c 

zero shear can be observed for the stations upstream of 
x/y c = 108.4, since it was not possible to make slant 
wire measurements closer to the wall at those stations. 

This was because of the size limitations set by the slant 
wire probe supports. The points of zero shear do not 
coincide with the points of zero velocity gradients. Both 
the inner and outer points of zero shear were found to be 
closer to the bottom wall than the corresponding points 
of zeio velocity gradient. This result is in agreement 
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Fig. 4.13(b) Tistribution of the Shear Stress at 

x/y = 74.48 to 435.2 
J c 

Shear stress data by Klebanoff (1955) for 

a turbulent boundary layer on a flat 
plate (reproduced from Rotta (1962)) 





with a similar fact reported in the literature on wall 
jet data. 
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The velocity maxima and minima cannot be identified 
beyond x/y = 146.2. Therefore, the -uv profiles beyond 
x/y = 146.2 do not show points of zero shear. The points 
of maximum shear stress between the inner point of zero 
shear and the freestream were found to be in the region 
of maximum velocity gradients. 

4.3.3 Effect of the Asymmetric Jet Velocity Profi. 1 e on 
The Turbulent Quantities 

The profiles of u t , v fc , w fc , and -uv, at x/y p = 435.2 
are compared with those of a turbulent boundary layer on 
a flat plate by Klebanoff (1955) in Figs. 4.9 (c) , 4.11 
(b) , 4.12 (b) , and 4.13 (b) respectively. The present 
profile of u^ at x/y c = 435.2 is significantly different 
from that of a flat plate boundary layer. A region of 
constant u^ can be found in the present data at x/y c = 

435.2. The profiles of v^ and w^_ at x/y c = 435.2 also 
differ considerably from those for a flat plate boundary 
layer. The influence of the jet seems to persist as far 
as x/y c = 435.2. This was indicated by the departure of 
the turbulence profiles from that of a flat plate boundary 
layer. Kacker and Whitelaw's (1968, 1971) experiments 
show that the distribution of the turbulence intensities 
approach normal flat plate boundary layer patterns in a 
much shorter distance for flows similar to the present 
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flow and with a uniform jet velocity profile. This supports 
the fact that an asymmetric profile carries the influence 
of the jet to a much longer distance than a uniform jet. 

The shape of the -uv distribution at x/y = 435.2 
shown in Fig. 4.13 (b) is not comparable to that of the 
flat plate boundary layer by Klebanoff. The -uv drops 
to a -very small value even at y/5 = 0.5, indicating that 
there are no appreciable velocity gradients beyond y/6- 
0.5. It means that the momentum of the jet is distributed 
evenly across the layer thereby reducing the momentum 
deficient region and the large velocity gradients associ- 
ated with it. This is caused by the asymmetric velocity 
profile in the jet which mainly distributes the momentum 
evenly across the layer instead of concentratina it in the 
region very near the wall as in the case of a uniform jet 
velocity profile. 

Kacker and Whitelaw's (1968, 1971) experiments 
on flows similar to the present flow show that the -uv 
distribution approaches the present distribution for 
x/y c = 435.2 at a much shorter distance. This indicates 
that the flow with a uniform jet velocity profile has a 
tendency to return to a normal boundary layer pattern at 
a much shorter distance than an asymmetric profile. Hence, 
the asymmetric profile carries the jet momentum to a much 
longer distance than a uniform profile, since one can see 
the effect of the jet as far as x/y c = 
of asymmetric profile. 


435.2 in the case 
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4.3.4 Turbulent Kinetic Energy k 

Figs. 4.14 (a) and (b) show the distribution of 
turbulent kinetic energy. The point of maximum kinetic 
energy coincides with the region of maximum u^, v t , and 
W.J. and also the region of maximum velocity gradients. A 
region of constant turbulent kinetic energy can be seen 
in the profiles starting from x/y c = 53.74. The turbulent 
kinetic energy decreases rapidly up to x/y c = 53.74; 
later on the decrease is more gradual. 

4.3.5 Correlation Functions R uv and a^ 

The correlation function R uv is defined as 

R uv = -uv/(u t v t ) . The correlation function a^, also 

called the Bradshaw's turbulence intensity parameter, 

2 2 2 

is defined as a 1 = -uv/(u t + v fc + These two coeffi- 

cients show the correlation between the different fluctuating 
components of velocity at the same point. The above 
correlation functions are very useful in understanding the 
structure of turbulence. The distributions of R uv and a^ 
for the present flow are given below. 

4.3.5 (a) Correlation Function R uv 

Figs. 4.15 (a) and (b) show the distribution of 
the correlation function R uy at different stations. This 
function is either negative or positive depending on the 
sign of -uv and it is zero at the point of zero shear stress. 
The value of R uv reaches a maximum value lying between 



Fig. 4.14 Distributions of the Turbulent Kinetic Energy 
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+0.4 and +0.5 for points beyond the outer point of zero 
shear. The minimum value of R^ v between the inner and 
outer points of zero shear is around -0.5 to start with, 
gradually decreasing in absolute value as the value of 
x/y c increases. R uv is always positive at x/y c = 220.7 
and beyond. At x/y c = 435.2, R uv is nearly a constant 
varying between +.35 and +.45 up to y/ 6 = 1.' 0 and then 
gradually decreasing to zero in the freestream. 

The distribution of R uv at x/y c = 435.2 is compared 
in Fig. 4.15 (b) with the distribution of R uv for a tur- 
bulent boundary layer on a flat plate (Klebanoff , 1955) . 

The distribution of R „ at x/y„ = 435.2 is similar in 

uv J c 

pattern to the turbulent boundary layer on a flat plate. 

The magnitudes of R uv at x/y c = 435.2 also compare well 
with those of a flat plate boundary layer. 

4.3.5 (b) Correlation Function a^ 

Figs. 4.16 (a) and (b) show the distribution of 
the correlation function a-^. The distribution of a-^ ex- 
hibits similar features as in the case of R uv - The value 
of a^ becomes zero when the shear stress is zero. The value 
of a 1 reaches a maximum value ranging between +.15 and 
+.18 for points beyond the outer point of zero shear. The 
minimum value of a^ between the inner and outer points of 
zero shear is around -0.2 to start with, gradually decreas- 
ing in absolute value as x/y c increases. At x/y c = 435.2, 
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is nearly a constant, varying between +.11 and +.15 up 
to y/S = 1.0 and then gradually decreasing to zero in the 
freestream. The distribution of a-^ at x/y c = 4 35.2 is 
compared in Fig. 4.16 (b) with the distribution of a^ for a 
turbulent boundary layer on a flat plate obtained by Hinze 
(1959) from the data of Klebanoff (1954). The distribu- 
tion of a^ at x/y c = 435.2 is similar in pattern to the 
turbulent boundary layer on a flat plate. The magnitudes 
of a^ at x/y c = 435.2 also compare well with those of 
a flat plate bounuary layer. 

4.3.6 Eddy Viscosity 

The study of the eddy viscosity distributions is 
very useful in the development of turbulence models and 
prediction procedures. Fig. 4.17 (a-c) shows the distri- 
butions of eddy viscosity v e £f = -uv/(3U/3y) 
non-dimens ionalised with (U 6.) . The trends of the ex- 

co i 

perimental results are indicated by the dashed lines. It 
can be easily seen that there is a region of singularity 
exhibited at each station up to x/y c = 146.2. The eddy 
viscosity becomes either negative or undefined in these 
singular regions because of the following reasons. 

As mentioned earlier (section 4.3.2 (b) ) , the 
points of zero shear do not coincide with the points of 
maximum or minimum velocity. Therefore, the value of -uv 
is different from zero and 3U/3y is equal to zero at the 
point of maximum or minimum velocity. Hence the eddy 



Eddy viscosity distribution for a turbulent boundary layer on a flat plate obtained by Hinze (1959) from 
the data by Klebanoff (1954) 
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viscosity becomes undefined at the point of maximum or 
minimum velocity. Also,, -uv and 3U/9y are of different 
sign in the region between the point of zero shear and the 
point of velocity maximum or minimum. Hence the eddy 
viscosity becomes negative there. 

The eddy viscosity finally drops to zero in the 
freestream. The eddy viscosity distributions at x/y c = 

108.4 and x/y = 146.2 exhibit two regions of singularity, 
as they are the only stations where the slant wire measure- 
ments cover both the points of maximum and minimum velocity. 
The above mentioned features of the eddy viscosity dis- 
tributions were also observed by Kacker and Whitelaw 
(1968) in their wall jet experiments. 

At x/y c = 435.2, the eddy viscosity is nearly 
constant over most of the layer up to y/<5 = 0.6. The eddy 
viscosity starts increasing after y/6 = 0.6 and it attains 
a maximum value near y/6 = 0.8. Afterwards, it gradually 
goes to zero in the freestream. The eddy viscosity dis- 
tribution at x/y c = 435.2 is compared in Fig. 4.17(c) 
with the eddy viscosity distribution for a flat plate tur- 
bulent boundary layer obtained by Hinze (1959) from the 
data of Klebanoff (1954). The residual effect of the jet 
may be responsible for the different patterns of the eddy 
viscosity distributions for the flat plate boundary layer 
and for the present wall jet at x/y c = 435.2. On the whole, 
there does not seem to be any similarity behavior in 
the present viscosity distributions. 
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4.3.7 Length Scales L . and L. 

mix K 

Two length scales were calculated from the present 
measurements. They are the Prandtl mixing Length L m j_ x > 
defined as L m ^ x = [v Q f f / (| |^| ) i 2 and the Prandtl-Kolmogorov 
turbulent length scale L^ defined as L k = v e ff // ( C y' / ^ * 

is a constant and k is the turbulent kinetic energy 
2 2 2 

equal to %(u t + v + w fc ) . Both these length scales are 
very useful in the prediction methods . 

4.3.7 (a) Prandtl Mixing Length L m ^ x 

Figs. 4.18 (a-c) show the distribution of the 
Prandtl mixing length non-dimensionalised with 6 . Here 
also, the singularity can be observed around the points 
of zero shear as in the case of the eddy viscosity. The 
trends of the experimental results are indicated by the 
dashed lines. 

At x/y c = 16.5, the mixing length is essentially 
constant in the outer layer beyond the outer point of 
zero shear which is due to the effect of upstream boundary 
layer. Two distinct regions, i.e., the outer and inner 
regions, can be found in the distributions of L , . . The 

Till Jv 

length scale in the outer region is in general approximately 
twice that of in the inner region, and this ratio is much 
higher for x/y c = 146.2 and beyond. At x/y c = 435.2, 
the value of L m ^ x / 6 in the inner region is nearly a con- 
stant equal to around 0.04 and in the outer region it is 
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around 0.14 with a smooth transition between them. In 
general the tendency of L m ^ x /6 in the outer region is to 
approach a uniform value of around 0.08 at higher values 
of y/<5 . However, the value of L m ^ x /6 approaches a uniform 
value of around 0.12 at x/y c = 257.9 and 435.2 for large 
values of y/6 . The mixing length distributions do not 
show any similarity behavior among themselves. Kacker and 
Whitelaw {1968, 1971) observed similar mixing length 
distirbutions in their experiments on wall jets and the 
order of magnitudes of L m ^ x /6 were found to be the same 
as the present data. 

4.3.7 (b) Prandtl-Kolmogorov Length Scale 

Figs. 4.19 (a-c) show the distribution of Prandtl- 
Kolmogorov length scale L^, non-dimens ionalised with 6. 

A value of = 0.2 (Kacker and Whitelaw, 1968, 1971) was 
used in the evaluation of L^.. The length scale ex- 
hibits similar features as the mixing length, except that 
that magnitude of is different. At x/y c = 146.2 and 
beyond, the magnitude of in the outer region is about 
four times that of the inner region. At x/y c = 435.2, the 
value of L^/ 5 is around 0.13 in the inner region and around 
0.35 in the outer region with a smooth transition in 
between them. The Prandtl-Kolmogorov length scale dis- 
tributions do not show any similarity behavior among 
themselves. Kacker and Whitelaw (1968, 1971) observed 



Fig. 4; 19 Distributions of the Prandtl-Kolmogorov Length Scale 
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similar Prandtl-Kolmogorov length scale distributions in 
their experiments on wall jets, and the order of magni- 
tudes of L^/6 were found to be the same as the present 
data . 

Both the Prandtl mixing length model and the 
Prandtl-Kolmogorov model fail near the point of zero 
shear. The length scales L . and as defined above, 
become either negative or imaginary around the point of 
zero shear. This is because v becomes negative around 
the point of zero shear due to the different signs of the 
velocity gradient and the shear stress. 

Therefore, there is a necessity for a turbulence 

model that is suitable for wall jets where the point of 

zero shear and the point of velocity maxima or minima do 

not coincide. Launder (1969) has proposed an expression 

for the effective viscosity based on the Prandtl-Kolmogorov 

2 

model which includes the second order terms like 3 U/3y“ 
to take care of the region around the point of zero shear. 
However, no attempt has been made here to test the applic- 
ability of the Launder ' s expression for the effective 
viscosity to the present flow. 

4.3.8 Production of Turbulent Kinetic Energy 

The equation for the consex'vation of turbulent 
kinetic energy (Rotta, 1962) involves the advection, 
production, diffusion and the dissipation terms. A full 
scale energy balance has not been attempted here. However, 
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an attempt has been made to evaluate the production and 

dissipation terms in the equation for the turbulent kinetic 

energy from the present measurements. The dissipation rate 

is evaluated from spectral measurements and the dissipation 

results are presented in the next section. The production 

of turbulent kinetic energy can be divided into normal 

stress production and the shear stress production. The 

shear stress production is given by -uv ■*— and the normal 

__ 6 Y 

2 2 3U 

stress production is given by (u - v ) Both the shear 

stress and the normal stress production terms were evalu- 
ated from the present measurements and the results are 
shown below. 

Figs. 4.20 (a-c) show the distribution of the 
shear stress production of turbulent kinetic energy 

QTJ 3 

given by -uv - 5 — ( 6 /U ) on a semi-logarithmic scale. It 
dy 00 

can be seen here that there are two distinct levels of 
turbulent energy production corresponding to the inner 
and outer regions . These two production levels differ by 
an order of magnitude to start with and they gradually 
become equal. At x/y c = 108.4, the production level in 
the region between the two points of zero shear is of the 
same order of magnitude as in the region beyond the outer 
point of zero shear. As mentioned earlier (section 4.3.6), 
there is a small region around the points of zero shear 

0-Q 

where -uv and - 5 — are of different sign. In this region, 
dy 

the shear stress production becomes negative and these 


X/YC= 7V*8 X/rC= ??□.» 



Fig. 4.20 Distributions of the Shear Stress Production of Turbulent Kinetic Energy 
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negative production points are not shown in the plots . 

The negative production terms were found to be in the order 
of -uv —• 0$/U^) = 0.3X10 At x/y c = 146.2 and there- 
after, the slant-wire measurements extend into the layer 
near the wall and the production is maximum near the wall 
at those stations. Also, there is a minimum production 
region representing small velocity gradients for stations 
x/y = 146.2, 220.7 and 257.9. For these stations, the 
production attains a maximum also in the outer region. 

At x/y c = 435.2, the distribution of shear stress produc- 
tion is almost linear in the semi-logarithmic co-ordinates 
with a maximum near the wall and the shear stress production 
gradually decreases to zero in the freestream. 

The normal stress production term obtained from 
the present measurements was generally found to be two 
orders of magnitude lower than the shear stress production 
term. The normal stress terms were found to be relatively 
higher for stations very close to the slot. 

4.3.9 Spectra and Dissipation Measurements 
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where n is the frequency in Hz. The spectrum function 
F(n) was obtained at several points across the layer at 
each streamwise location. 

Figs. 4.21, 4.22, and 4.23 show various spectra of 

— 

u representing the present flow at various stations. 

The spectra shown are smooth lines drawn through the ex- 
perimental curves. Only representative spectra are 
presented here. F(n) is presented here in the non- 
dimensional form of Nq F(n). Nq is a constant of the 
spectrum analyser equal to 20 Hz for the 10 KHz range. 

The magnitudes of Nq F(n) are not given on the Nq F(n) 
axis. The position of each spectrum with respect to the 
Nq F(n) ordinate can be fixed with the aid of Table (4.2) , 
showing the magnitude of N Q F(n) at 1 KHz for each spectrum 

that is presented in Figs. 4.21, 4.22 and 4.23. 

T 

The u spectra were taken in the present case with 
two points in view: 

1. To find out the range of frequencies where 

n F(n) is constant, (i.e., F(n) « n which correspond to 
the most energetic frequencies and 

2. To locate the range of frequencies where F(n) « 

-5/3 

n • 

The two points given above lead to the determination of the 
dissipation rate e as discussed later in this section. 

Figs. 4.21 (a) and (b) show the spectra at x/y c = 
0.292 for different values of y/y . It can be seen from 















TABLE 4.2 MAGNITUDES OF Nq F(n) at 1 KH 

FOR THE SPECTRA PRESENTED IN 
FIGS. 4.21, 4.22, 4.23 
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Pig. 4.21 (a) that inside the jet a considerable amount 
of energy was contained in the high frequency range be- 
yond 3 KHz. A series of high frequency peaks were 
observed very near to the slot lip in the range of 3 to 
10 KHz. These peaks diminished away from the lip. One 
of the spectra taken near the lip has a peak at 756 Hz. 

This spectrum was not presented here. However, the energy 
contained in that peak was relatively small compared to 
the total energy under che spectral curve. As shown in 
Fig. 4.21 (a), the spectra inside the jet exhibit two 
ranges of n ^ slope. The first range corresponds to the 
low frequencies below 3 KHz and the second range corre- 
sponds to the high frequencies beyond 3 KHz. Beyond the 
lip, the spectra at x/y c = 0.292 (Fig. 4.21 (b) ) are similar 
to that of a normal turbulent boundary layer representing 
the upstream boundary layer. The region of n ^ slope was 
in the high frequency range for spectra immediately above 
the lip corresponding to the wall region of the upstream 
boundary layer. The energy content of the higher frequen- 
cies decreased away from the lip. The -5/3 slope was 
observed over a greater region of the spectrum for spectra 
far away in the outer layer . 

Fig. 4.22 (a) and (b) show the spectra representing 
the inner layer at various stations starting from x/y c = 
7.45. A peak was observed at 4.3 KHz for x/y c = 7.45 
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and at the y/> indicated. The range of -1 slope slowly 
shifts from high to low frequencies as one goes downstream 
in the inner layer. It should be noted that this shift 
is only relative and from an absolute viewpoint, the range 
of -1 slope at x/y c = 182. 8 in the inner layer still lies 
on the high frequency side around 1 KHz. The -5/3 slope 
was found over a greater region in the inner layer spectra 
at stations farther from the jet than at stations nearer 
to the jet. 

Fig. 4.23 (a) and (b) show the spectra represent- 
ing the outer layer at various stations starting from 
x/y c =7.45. It can be clearly seen that at all stations 
a considerable region of -5/3 slope exists. The region 
of =1 slope can be found only in the low frequency range 
between 200 to 700 Hz. 

It has been found in general that more energy is 
contained in the higher frequencies for the spectra near 
to the wall. As the value of y increased towards •che 
freestream, the energy content of the higher frequencies 
decreased. This result is in accord with findings of 
Klebanoff (1954) , for the developed zero pressure gradient 
turbulent boundary layer. The spectra also indicated that 
the -5/3 slope is found over a greater region of the spec- 
trum for points away from the wall than close to the wall. 
This result is also in accord with the findings of 
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Klebanoff (1954) for the developed zero pressure gradient 
turbulent boundary layer flow. 

4.3.9 (b) Dissipation Rate e 

The dissipation rate of turbulent kinetic energy 
was estimated from two different methods, the -5/3 spec- 
tral law and the Tchen's high mean velocity spectral model 
(Hinze, 1975) . The -5/3 law of the inertial subrange 
states that 


u 2 P (n) = 


K e 


2/3 


,2t L -2/3 
^ U ‘ 


n 


-5/3 


(4.4) 


where K is a constant equal to 0.49 (Corrsin, 1964; 
Bradshaw, 1967a) and U is the local mean velocity. The 
-5/3 law was used to obtain the dissipation rate in places 
where a clear -5/3 region was found. The same method was 
used even when clear -5/3 region was not found, by draw- 
ing a tangent of appropriate slope to the spectrum. How- 
ever, the region of -5/3 slope did not exist at all for 
some spectra near the jet and in the inner layer. No 
attempt has been made to evaluate e using -5/3 law in those 

cases. Bradshaw (1967b) suggested that the turbulence 

2 2 }- 

Reynolds number Re = ( u ) X/v where \= (15vu /e) must 

A 

be greater than 100 for an inertial subrange to exist. 

In the present case, 25 < Re^ < 98 for most spectra where 
-5/3 region has been found. 

Tchen's high mean vorticity model (Hinze, 1975) 
relates F(n) to the turoulent dissipation e by 
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00 


where f = U/U and n = 
00 


- a [n F(n) ] 
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T 

u 

u 2 


3f 

an 


(4.5) 
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and a is a constant . 


In the present results a = 0.8 was used, which is in accord 
with the value of a for IClebanoff's (19 55) flat plate 
boundary layer. Here the value of nF(n) to be used in 
the equation (4.5) is the value corresponding to -1 slope. 

The dissipation rate e was also evaluated using the 
equilibrium relationship 


(4.6) 

which naturally holds good only in the logarithmic velocity 
profile region where turbulence dissipation equals produc- 
tion (Rotta, 1962) . The -uv measured from slant-wire 
data was used in this method of evaluating e. 

Figs. 4.24 (a), (b) and (c) show the distribution 

4 

of ev/U^ at various stations starting from x/y c = 0.292. 
The three different values of e obtained from the three 
different methods mentioned above were plotted in these 

4 

figures. It should be noted that ev/U^ is plotted on 
a full logarithmic scale against y/6 . At points where 
two ranges of n ^ slope have been found, the values of e 
obtained from Tchen's theory for the second range corre- 
sponding to high frequencies were shown as solid symbols . 

In general e has a maximum value near the wall and then 
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it decreases rapidly by about three orders of magnitude 
as one goes towards the freestream. Tchen's theory gives 
higher values of e as compared to the -5/3 law at points 
very near to the wall. Beyond y/6 = 0.01, Tchen's theory 
agrees very well with the -5/3 law although there is some 
scatter in the far region of the outer layer. This region 
where Tchen's theory agrees with the -5/3 law generally 
corresponded with the logarithmic region in the velocity 
profile . 

The equilibrium assumption between production and 
dissipation was generally found to give lower value of e 
than either Tchen's theory or the -5/3 law. However, 
most of the slant wire data used in the evaluation of e 
using equilibrium assumption were taken beyond y/6 = 0.1 
where no logarithmic velocity profile exists and the valid- 
ity of the equilibrium assumption there is doubtful. 

, 3/2 

The dissipation length scale L £ = k ' /e was 

calculated using the value of e obtained from the -5/3 law. 
Here k is the turbulent kinetic energy. No definite 
pattern has been observed in the distribution of . 
However, it starts with a low value in the inner region 
and gradually increases as one goes towards the freestream. 
4.3.9 (c) Bursting Frequency 

Strickland and Simpson (1973, 1975) have shown that 
there is a one to one correspondence between the bursting 
frequency of wall shear stress spectra in the turbulent 
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boundary layer and the peak of the nF(n) spectral dis- 
tribution. This principle is used in the present case to 
obtain the bursting frequencies across the layer. However, 
most of the spectra taken here have shown a range of fre- 
quencies where -1 slope is valid instead of a defined 
peak. Therefore, it is difficult to select a precise 
single frequency at which nF(n) is a maximum for a given 
spectral distribution. Hence the bursting frequency was 
chosen as the frequency corresponding to the center of the 
range of frequencies over which -1 slope is valid. 

The bursting frequency was normalized with U ro 
and 5 giving rise to a non-dimensional bursting period 
U^/n^S , where n^ is the bursting frequency. Figs. 

4.25 (a) , (b) and (c) show the distributions of the 

bursting period in semi-logarithmic co-ordinates. For 
spectra where two ranges of -1 slope have been observed, 
the bursting periods for the second range corresponding 
to the high frequencies were shown as solid symbols. It 
can be seen that the bursting period starts with a value 
lying between 0.5 and 2.0 for points very close to the 
wall (y/6 < .01) and then decreases rapidly to about half 
of its starting value at a point somewhere in between 
y/6 = 0.01 and y/6 = 0.1. After y/6 = 0 . 1, the bursting 

period increases rapidly to values up to 4.0. This means 
that in the inner layer most of the energy is contained in 
the high frequency range. In comparison, the outer layer 



Fig. 4.25(a) Distribution of the Bursting Period at x /y c - 0.292 to 35.05 
The solid symbols correspond to the second n' 1 range (high frequency) 





Fig. 4.25(c) Distribution of the Bursting Period at x/y = 146.2 to 257. 


has the most energetic frequencies in the low frequency 
range. However, within the inner layer two ranges of 
bursting frequencies can be found, one being for points 
very close to the wall (y/6 < .01) and the other for 
points lying between 0.01 < y/6 < 0.1. 

It was mentioned above that the bursting period 

in the region .01 < y/6 <0.1 was found to be half of 

its value for points very close to the wall (y/6 < .01). 

The region 0.01 < y/6 < 0.1 corresponds to approximately 
+ + U x y 

20 < y < 300 where y = — — . The first spectral data 
point for most of the stations corresponds with a y+ value 
approximately equal to 9 and was taken at the closest 
possible distance from the wall. Therefore, the bursting 
period for points lying between 20 < y + < 300 is half 
of that at y + = 9 . This is in agreement with a similar 
result for plane turbulent boundary layer by Ueda and Hinzc 
(1975) . 

At x/y c = -77.1, the bursting periods for y + = 10. 
and y + = 95.0 were found to have values around 4.8 and 2.3 
respectively. These bursting periods for the upstream 
boundary layer at x/y c = -77.1 also agree with those 
reported by Ueda and Hinze (1975) for a plane turbulent 
boundary layer. 



CHAPTER V 


EXPERIMENTAL RESULTS FOR THE ADVERSE 
PRESSURE GRADIENT FLOW 

The experimental results for the adverse pressure 
gradient flow are presented in this chapter. The objec- 
tives of the experiments are: 

1. To obtain experimental data on a wall jet 
under an adverse pressure gradient with a thick upstream 
boundary layer and with a low jet velocity ratio; 

2. To obtain turbulence data on such a wall 
jet flow for future use in the development of turbulence 
models ; and 

3. To further observe the effect of an asymmetric 
jet velocity profile on the flow development downstream 

of the slot under an adverse pressure gradient. 

The measured quantities presented here are the 
mean velocity, Reynolds stresses, and skin friction. 

The derived quantities presented include the integral 
and profile parameters, eddy viscosity, mixing length, 
Prandtl-Kolmogorov length Scale, turbulent kinetic energy, 
correlation coefficients and the rate of production of turbu- 
lent kinetic energy. A brief description of the flow conditions 
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is first given, followed by the presentation of the meas- 
ured and derived experimental results . 

5 . 1 Flow Conditions 

The selection of the required adverse pressure 
gradient and the method of setting the pressure gradient 
are described below in detail. 

One of the important considerations in selecting 
the proper pressure gradient or the external velocity 
distribution was that it should represent a typical 
free-stream velocity distribution that occurs in practice 
such as on aerofoils and jet flaps. An external velocity 
distribution which has a steep decrease in velocity ini- 
tially followed by a more gradual decrease is generally 
typical of the external velocity distributions on aero- 
foils and jet flaps (Irwin, 1974) . The external velocity 
falls between 0.6 and 0.3 of its value at the slot in 
normal external velocity distributions on jet flaps 
(Irwin, 1974) . The other considerations in selecting the 
proper pressure gradient were: 

1. To keep the ratio of the jet velocity to the 
free-stream velocity around 1.5 because it falls in the 
practical working range, especially for high speed flows; 

2. To have the most severe pressure gradient that 
the flow can withstand without flow reversal at the above 
mentioned ratio of the jet velocity to the free-stream 
velocity . 
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The strength of the pressure gradient can be 
increased by increasing the level of the static pressure 
inside the tunnel. One can easily see that the higher 
the static pressure level inside the tunnel, the greater 
is the amount of flow through the perforations in the top 
wall and hence the stronger is the perssure gradient. 
However, the capacity of the blower limits the static 
pressure level that can be maintained inside the tunnel 
without drastically reducing the free-stream velocity 
at the slot. 

In order to satisfy the above requirements, a 
trial and error method had to be adopted to set up the 
proper pressure gradient. The flow through the perforated 
top wall was tailored by adding or removing duct tape 
of proper width uniformly across the upper wall of the 
tunnel at several streamwise locations. The perforated 
top wall was completely covered with the duct tape from 
the contraction exit to the slot. Also, the last 51 cm 
of the top wall was completely covered with duct tape. 

The height of the top wall (Fig. 5.1) above the bottom 
wall was approximately the same as that for the zero 
pressure gradient case except that it slightly increases 
with distance in the downstream direction. This gradual 
increase in the height of the top wall contii s -the 
process of reducing the free-stream velocity even after 
the flow bleeding stops in the last 51 cm of he tunnel 
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length. The static pressure level inside the tunnel v T as 
tailored by blocking the perforations in the end plate 
with a duct tape. 

Flow visualization studies were first performed 
using dried tea leaves and tufts to observe flow reversal 
for a given distribution of the duct tape on the top 
wall and the end plate. These studies have shown that the 
sidewall boundary layers were separating and prevention 
of it became the main concern. The side wall boundary 
layers were then sucked off by using false plexiglas 
side walls with a sharp leading edge. These false side 
walls are 4.8 mm in thickness and are attached to the 
aluminum angle that reinforces the top wall. The bottom 
edges of the false side walls were sealed against the bot- 
tom wall with foam weather strip. The distance between 
the false side walls was 19.4 cm compared to the distance 
of 24 cm between the glass side walls. The false side 
walls start at 21 cm from the slot and extend all the way 
to the end of the tunnel. The flow between the false 
side walls and the actual glass side walls was bled from 
the top wall through 13 cm long openings located at 95 cm 
from the slot and on either side of the center of the 
tunnel . The remaining portions of the top wall between 
the false and actual side walls was completely covered 
with duct tape. Upstream of the false side wall, the 
duct tape strips extended across the tunnel from one 
glass side wall to the other. 
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At this stage, several trials were made with the 
distribution of the duct tape on that portion of the 
top wall where the main flow was being bled and on the 
end plate. These trials were made in view of the three 
important considerations mentioned at the beginning of 
the chapter. Flow visualization studies were simultan- 
eously made with tufts to observe any flow reversal. Figs. 
2.2 and 2.4 show the final positions of the duct tape on 
the top wall and the end plate respectively. After fix- 
ing all the different parameters involved, the wall jet 
velocity was increased slightly to prevent any flow sep- 
aration due to the flow blackage by the measuring probes. 

The final free-stream velocity at the slot was 
approximately 18.3 m/sec and the maximum velocity in 
the wall jet velocity profile was 30.4 m/sec. The wall 
jet velocity was kept constant during the experiments by 
keeping a constant static pressure difference between 
the inside of the jet nozzle and the frees tream at the 
slot. The free-stream velocity at different x locations 
was calculated from the static pressure measurements in 
the freestream, assuming that the total pressure was 
constant along a streamline in the freestream. The final 
free-stream velocity distribution along the tunnel is 
given in Fig. 5.2 There is a steep decrease in the ex- 
ternal velocity initially followed by a more gradual 
decrease. The external velocity fell to 0.4 times its 


Fig. 5.2 Free-Stream Velocity Distribution 
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Fig. 5.3(a) Variation of Skin Friction 
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value at the slot in about 360 slot heights. The flow 
upstream of the slot was maintained under zero pressure 
gradient. It should be mentioned here that the flow was 
separating from the bottom wall at a point very close 
to the slot when no fluid was introduced through the jet 
and the entire downstream flow on the lower wall was 
separated. 

5 . 2 Mean Flow Data 

The mean flow data presented here include the skin 
friction, mean velocity profiles and the integral and pro- 
file parameters. 

5.2.1 Skin Friction 

The variation of skin friction coefficient 

obtained from the Preston tube measurements is shown in 

Fig. 5.3 (a) . Fig. 5.3 (b) shows the variation of C f with 

x/y c on a logarithmic scale. It can be seen that there 

is a power law variation between C- and x/y up to x/y = 

X c c 

143.0 . 

It has been found in the literature (section 4.2.1) 

that for self-preserving wall jets under pressure gradients, 

2 

the skin friction C^ m = t w / ( U ) varies according to 

the equation (4.2) . Fig. 5.3 (c) shows the variation of 

with U__, Y /v for the present case. No definite 
fm max max c 

pattern is found in the distribution of C fm except that 
it follows a power law in the initial region. The failure 
of Cf m to follow a power law all the way unlike the case 


CFm*10 




of zero pressure gradient flow can be explained in the 
following way. Because of the severity of the adverse 
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pressure gradient, the value of Y max increases very rapid- 
ly at far downstream stations. Therefore, at farther 
downstream stations, the Y max point cannot be considered 
as part of the inner layer and hence a defined correlation 

between C_ and U m , Y , /v cannot be found, 
fm max max 

5.2.2 Mean Velocity U 

Figs. 5.4 (a-c) show the normal hot-film meas- 
urements of mean velocity profiles in the U/U^vs y (cm) 
co-ordinates. The negative sign on the x/y c value in 
Fig. 5.4 (a) indicates that it is the distance measured 

upstream of the slot. The value x/y c = -16.4 corresponds 
to the crest of the extruded aluminum end piece (section 
2 . 2 . 2 ). 

Fig. 5.4 (a) also shows the velocity profile at 
x/y c = 0.292. Here it can be seen again that the velocity 
profile in the jet is asymmetric with a relatively greater 
concentration of momentum in the upper half of the wall 
jet. We will see in the later sections of this chapter 
how this asymmetric jet velocity profile helps in better 
control of separation compared to a uniform velocity 
profile in the jet. The upstream boundary layer at the 
jet (Fig. 5.4 (a)) is very thick and has a large deficit 
of momentum satisfying one of the proposed conditions 
under which the present measurements were intended to be 
made as given in Chapter 1. 
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Predictions using experimental initial profiles 
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Figs. 5.4 (b) and (c) show the mean velocity 

profiles from x/y c = 15.92 and onwards. It should be 

mentioned here that the entire height of the tunnel was 

almost completely filled with the shear layer at x/y = 

c 

177.7 and onwards. Therefore, the experimental data at 
x/y c = 177.7 and beyond have to be treated with caution 
from every point of view. The data beyond x/y c = 177.7 
are especially useful in finding out the effectiveness 
of an asymmetric jet velocity profile in the control of 
separation. The boundary layer thickness 6 could not be 
determined beyond x/y c = 143.0 since no frees tream was 
observed. Hence, the velocity profiles were plotted 
against y in cm. The free-stream velocity used in the 
data reduction at stations beyond x/y c = 143.0 was only 
an imaginary one that would produce the measured static 
pressure in an irrotational frees tream. 

As mentioned in section 5.1, the present pressure 
gradient was the most severe one that the flow could with- 
stand with a ratio of jet velocity to free-stream velocity 
at the slot around 1.5. In other words, any stronger 
pressure gradient beyond the present case, keeping the 
jet velocity to free-stream velocity ratio constant, 
would lead to flow separation downstream. Because of 
the severe pressure gradient imposed on the flow, the 
shear layer rapidly became thicker and eventually filled 
the whole tunnel as shown in Figs. 5.4 (b) and (c) . The 
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upstream boundary layer was not absorbed completely by 
the jet, unlike in the case of the zero pressure gradient 
flow. The velocity maxima was observed throughout the flow 
and the velocity minima could be observed until x/y = 
285.9. 

Fig. 5.5 (a) and (b) show the velocity profiles 
plotted on a conventional semi-logarithmic plot with 
U/U T and U y/v as the co-ordinates. The velocity pro- 
files are shown for stations x/y = 15.92 and onwards. 

The shear velocity U T was obtained from the Preston tube 
measurements . The velocity profile measurements were 
compared with the logarithmic law of the wall, 

u u T y 

0 - = 5.5 Log.Q -i- + 5.45 (5.1) 

with the constants recommended by Patel (1965) . It can 
be seen from Figs. 5.5 (a) and (b) that the experimental 
data agree well with the logarithmic law of the wall up 
to x/y c = 143.0. The exception is at station x/y c = 

15.92, where a defined logarithmic region is not present 
because of its proximity to the slot. A defined logarith- 
mic region can be observed even beyond x/y c = 143.0. 
However, the constant in the logarithmic law of the wall 
has to be different from 5.45 to fit the data. At far 
downstream stations, the data extended considerably into 
the viscous sublayer. The mean velocity data are tabulated 
in Appendix F along with the data of u. . 
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Fig. 5.5(b) Mean Velocity Profiles in the Wall Coorc 
at x/y^ = 177.7 to 360.8 

Logarithmic law of the Wall (Equation 5.1) 



5.2.3 Integral Parameters 


The integral parameters evaluated are the boundary 
layer thickness ( 6 ) , displacement thickness ( 6 ^ ) , momentum 
thickness ( 62 ) / shape factor (H) and the momentum thick- 
ness Reynolds number (Re 2 ) . 

Figs. 5.6 (a-e) show the development of the in- 
tegral parameters 6 /y c ., 5 1 /y c , 6 2 /y c , H and Re 2 
respectively. The value of 6 beyond x/y = 143.0 was 
taken equal to the value of y at the last data point. 
Because of the severe adverse pressure gradient imposed 
on the flow, 6 -^ and 6 2 increased very rapidly in the ini- 
tial region up to x/y - 177.7. The displacement thickness 
increased more rapidly than the momentum thickness up to 
x/y c = 177.7. The growth of 6 ^ and i $ 2 decreased beyond 
x/y c = 177.7. This was partly due to the fact that the 
outer edge of the outer layer was already touching the 
top wall of the tunnel at x/y c = 177.7 and hence the full 
velocity profile was not accounted for in the evaluation 
of 6 -^ and 6 2 at x/y c “ 177.7 and beyond. The reduction 
in the growth rates of 6 ^ and o 2 at far downstream stations 
was also partly due to dU^/dx being small there compared 
to the region near the jet. 

The two-dimensional nature of the flow was exam- 
ined by applying the two-dimensional integral momentum 
equation (5.2) to the present data. 




Fig. 5.6(b) Growth of the Displacement Thickness 
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(U 2 6„) = (U 2 6 0 ) 

co 2 co 2 


2 ' f /yC P~l 

” 2 x/y -15.92 “ 2 / 1 t dx J 

*x/y c =15 .92 x/y c =15.92 


- I 2 2 

2 u -v 


b U co 


.x/y c 


c=15 . 92 


7: ~1 ~I 2 

The contribution of the normal stress terms /[ (u -v ) /U 1 dy 

•AD 

to the integral momentum equation was found to be smaller 
than 5% and hence it was neglected. The measured distribu- 
tions of the skin friction coefficient, external velocity 
and displacement thickness were used in the equation 5.2. 
Fig. 5.6 (e) shows the values of Re 2 = U (X> &2 //v °btained 
from the two-dimensional integral momentum equation (5.2) 
along with the experimental values of Re^. The values of 
Re 2 obtained from the equation (5.2) agree with the ex- 
perimental values of Re 2 within 15% on the average up 
to x/y c = 177.7. The failure of the experimental data of 
Re 2 to agree with the values of Re 2 obtained from the 
integral momentum equation beyond x/y = 177,7 was mainly 
due to the fact that the entire tunnel height was filled 
by the shear layer beyond x/y c = 177.7. Because the shear 
layer was spread over the entire tunnel height, the full 



velocity profile was not accounted for in the evaluation 
of Re 2 from the experimental data. 

5.2.4 Effect of the Assymetric Jet Velocity 
Profile on the Mean Flow 

The effect of the asymmetric jet velocity profile 
on the mean flow development is discussed below. A pre- 
liminary comparison of the present flow with similar flows 
in the literature (Irwin, 1974; Goradia and Colwell, 1971; 
and Ramaprian, 1973, 1975) has shown that the present flow 
was able to withstand much more severe adverse pressure 
gradient without separation even with a modest jet momen- 
tum. Even though a direct comparison cannot be made, 
there was an indication from the references mentioned 
above that with the present pressure gradient and the 
injection velocities, the flow will separate downstream 
if a uniform velocity profile was used in the jet. 

The asymmetric jet velocity profile used in the 
present case was mainly responsible for the prevention of 
separation even under the severe adverse pressure gradient 
employed and the low injection velocities used. As 
mentioned earlier (section 4.2.4), the asymmetric jet 
velocity profile helps in distributing the jet momentum 
evenly across the layer by supplying more momentum to the 
momentum deficient upstream boundary layer than a uniform 
profile. The asymmetric profile also reduces the momentum 
losses due to friction at the wall by keeping a higher 
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concentration of jet momentum in the upper half of the 
jet than the lower half. In comparison, the uniform jet 
velocity profile tries to retain more momentum near the 
wall instead of supplying sufficient momentum to the mo- 
mentum deficient upstream boundary layer. This might 
result in the development of a large momentum deficient 
region downstream near the point of minimum velocity and 
eventual flow reversal there. The relatively large con- 
centration of momentum near the wall results in greater 
frictional losses at the wall in the case of a uniform 
profile. The above conditions might lead to eventual flow 
separation at the bottom wall in the case of the uniform 
jet velocity profile. An additional proof of the advan- 
tages of the asymmetric jet velocity profile over the 
uniform jet velocity profile was obtained from the theo- 
retical prediction of the present flows using uniform 
and asymmetric jet velocity profiles (Chapter 7 ) . 

The development of 6^ and &2 f° r the present 
flow was compared with that of a similar flow by Irwin 
( 1974 ) . Irwin's flow has a ratio of jet velocity to 
free-stream velocity of the same order and a relatively 
less severe pressure gradient. An exact comparison cannot 
be made because the exact, flow conditions are different 
in each case. The increase in 6-^ and 6 9 for a given 
x/y c distance was found to be lower for the present flow. 
The asymmetric jet velocity profile is mainly responsible 
for the slow growth of 6^ and 62 i- n the present flow. 
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The asymmetric profile distributes the momentum more 
evenly than a uniform jet and prevents build-up of a 
large momentum deficient region downstream near the min- 
imum velocity point, thereby resulting in slow growth 
of 6^ and S 2 . In comparison, the uniform jet velocity 
profile might resu .t in a large momentum deficient region 
downstream near the minimum velocity point as mentioned 
above. This will lead to large values of 6^ and 6 2 in 
the case of a uniform jet velocity profile, which is an 
unwanted effect in the efforts to prevent separation. 
5.2.5 Profile Parameters 


The mean velocity profile parameters presented 

The sig- 
nificance of these parameters was given in section 4.2.5. 


here include Y Y. n Y . , U and U . 

max, half, mm max mm 


No minima in velocity was found at x/y c = 360.84. There- 
fore, the last data point near the top wall was taken as 
the point of minimum velocity. 

5.2.5 (a) Development of Y , Y, . _ and Y . 

max half mm 


Fig. 5.7 (a) shows the development of Y max and 


Y half* 5. 7 (b) shows the development of Y 


mm 


Fig. 


5.7 (a) also shows the development of Y max for the zero 

pressure gradient flow. It can fc seen that the rate of 

growth of Y_„ is almost the same for both zero and ad- 

verse pressure gradients and the pressure gradient does 

not seem to influence the Y development. The universal 

max c 

Y development found in the literature (section 4.2.5 (a) 
max 
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Fig. 5.7(a) Development of and 
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for self-preserving flows with and without pressure grad- 
ients was also shown in Fig. 5.7 (a). As mentioned in 
section 4.2.5 (a), it is not clear whether the 
non-self-preserving wall jets exhibit universal behavior 
for the Y max development. However, the increased Y 
growth rate in the present case as compared to the univer- 
sal line for self-preserving flows can be partly attributed 
to the asymmetric jet velocity profile. The physical ex- 
planation for the increased Y m growth rate for an 
asymmetric jet velocity profile was given in detail in 
Chapter 4 under the section 4.2.5 (a) . The development 
of Y^alf is affected considerably by the pressure gradient. 
The growth rate of Y half is higher than that for the zero 
pressure gradient case and also it is higher than the 

gr-owth rate of Y for the adverse pressure gradient case. 
3 max r 

A comparison of the ratio of the growth rates 
of Y max and for the adverse pressure gradient flow 

with that of similar flows in the literature (Irwin, 1974; 
Ramaprian, 1973; Gartshore and Newman, 1969) shows that 
the present growth rate of is n ot very large compared 

to the Y max growth rate. The relatively slow growth rate 
of Y half can be a t tr ibuted to the asymmetric jet velocity 
profile. The physical explanation for the relatively slow 
growth rate of Y^ a ^ for an asymmetric profile was given 
in detail in Chapter 4 under the section 4.2.5 (a) . 


5.2.5 (b) Variation of U „ and U . 

max mm 

Figs. 5.8 (a-c) show the variation of U /U, 

max' Jave, 

U min /U Jave' and (U max~ U min ) /U Jave respectively with x/y c . 
U Jave ■’" S tlie jet velocity for an equivalent jet 

with an uniform profile and having the same momentum as 
the experimental asymmetric jet. Wall ;jet data in the 
literature show that the variation of U and U . depends 
on several parameters, such as the ratio of the jet veloc- 
ity to the free-stream velocity, the pressure gradient, 
and the conditions upstream of the slot. Hence an attempt 
has not been made to compare the present variation of U max 
and U . with the data in the literature. However, the U 
decay seems to follow a power law variation (Equation 5.3) 
as shown in Fig. 5.8 (a) . 


U 

max 

U, 

Jave 


3.513 (x/y c )“- 474 


(5.3) 


5.3 Turbulence Data 

The measured turbulence data include the turbulence 
intensities (u^/U^, v t /U 00 , and w i _/U 00 ) and the turbulent shear 
stress -uv. The quantities derived from the measured tur- 
bulence data include, the eddy viscosity, mixing length, 
Prandtl-Kolmogorov length scale, turbulent kinetic energy, 
correlation coefficients and the rate of production of 
turbulent kinetic energy. The results of the measured and 
derived turbulent quantities are given below in detail. 
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5*3.1 Turbulence Intensity u./U 
t, 00 

The u fc data were obtained from the normal hot-film 

traversals at several stations starting from x/y = 0.29 2. 

c 

Figs. 5.9 (a-c) show the variation of u.j./U^ vs. y(cms) . 

The negative sign on the value of x/y c indicates that 
those stations are located upstream of the slot. The pos- 
ition of x/y c = -16.4 is the same as given in section 
5.2.2. 

Figs. 5.10 (a) and (b) show the variation of u,. 

t 

in wall coordinates. For the stations near the slot, 
all of the data were out of the viscous sublayer, whereas 
for far downstream stations the data extend considerably 
into the sublayer. The point of first maximum in u^_ is 
always near the wall and at the outer edge of the sub- 
layer. The point of minimum u t in the inner layer gen- 
erally corresponds with the region of maximum velocity 
even though the exact position was generally nearer to the 
wall than the point of maximum velocity. The point of 
second maximum in u^ was found to be between the points 
of velocity maximum and the velocity minimum and nearer 
to the point of maximum velocity gradient. The point of 
minimum u fc in the outer layer was found to be around the 
point of minimum velocity. The steep drop in u fc far away 
in the outer layer was found to begin from the point of 
maximum velocity gradient between the velocity minima and 
the freestream. Some of the above features are not applicable 
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to u fc profiles beyond x/y c = 177.7, since the outer edge 
of outer layer was already touching the top wall at x/y c = 
177.7. Most of the above mentioned features of the u t 
profiles have been observed in the case of zero pressure 
gradient flow also. The u t data are given in Appendix F. 
5.3.2 Turbulence Intensities v^/U^,, w t /U oo 

And the Shear Stress -uv 

A rotatable slant-wire was used to obtain v^, 
w^ and the shear stress -uv, as in the case of zero pres- 
sure gradient flow (section 4.3.2). The details of 
obtaining the turbulence intensities and the shear stress 
from the slant wire data are given in Appendix C. The 
data of v^, w^, and -uv are tabulated in Appendix G. 
Appendix G also shows the interpolated data of u t and 
9 (U/U^) / 3 (y/S) obtained from the normal hot-film data. 
Because of the size limitations set by the slant wire 
probe supports, the slant wire data at some of the stations 
were taken only above the point of the velocity maximum. 
Slant-wire data were obtained only at four stations, 
compared to a greater number of slant-wire data stations 
for the case of zero pressure gradient flow. As mentioned 
in the Appendix C, the angle between the flow vector and 
the axis of the probe in the vertical plane was higher in 
the case of the adverse pressure gradient flow than for 
the zero pressure gradient flow. Details on the accuracy 
of the slant-wire data with respect to the flow inclina- 
tion are given in the Appendix C. The results of the 









turbulence intensities v t /U eo , , and the shear stress 

-uv are described below. 

5.3.2 (a) Turbulence Intensities v^/U^ and w t /U 0o 

Figs. 5.11 and 5.12 show the variation of v t /U 00 

and w t /U 00 at different stations. Generally, the positions 

of the maxima and minima in v. and w. were not well defined 

t t 

with respect to other variables . The maxima in v t and 
w^ were found to be nearer to the region of maximum ve- 
locity gradients for the stations x/y c = 15.92 and 44.39. 

5.3.2 (b) Shear Stress -uv 

Figs. 5.13 (a) and (b) show the distribution of 
-uv at different stations. The significance of the inner 
and the outer points of zero shear was explained in Chap- 
ter 4 under section 4.3.2 (b) . In the present -uv data, 
the inner point of zero shear was observed only at x/y c = 

100.3 and 177.7. It was not possible to make slant wire 
measurements up to the inner point of zero shear for sta- 
tions x/y = 15.92 and 44.39. This was mainly because of 
size limitations set by the slant-wire probe supports. 

The outer point of zero shear is shown in Figs. 5.13 (a) 
and (b) only for stations x/y = 15.92, 44.39, and 100.3. 
The outer point of zero shear for x/y = 177.7 can be 

V 

observed nearer to the last data point in the tabulated 
data (Appendix G) and it was not shown in the plots. 

The inner point of zero shear was observed to be closer 
to the bottom wall than the point of zero velocity 
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gradient corresponding to the velocity maximum. The outer 
point of zero shear was also found to be closer to the 
bottom wall than the point of zero velocity gradient cor- 
responding to the velocity minimum for stations x/y = 
15.92 and 44.39. The above results concerning the posi- 
tion of the points of zero shear are in agreement with 
similar results reported in the literature on the wall 
jet data. However, the outer point of zero shear at 
x/y c = 100.3 and 177.7 was found to be farther from the 
bottom wall than the point of zero velocity gradient 
corresponding to the minimum velocity. The slant-wire 
data were not corrected for the flow inclination (Appendix 
C) . The flow inclination angle near the outer point of 
zero shear was found to be higher at x/y c = 100.3 and 
177.7. Therefore, the increased flow inclination angle 
may be responsible for the observed discrepancy in the 

position of the outer point of zero shear at x/y = 100.3 

c 

and 177.7. The points of maximum shear stress were found 
to be in the region of maximum velocity gradients between 
inner and outer points of zero shear and also between the 
outer point of zero shear and the freestream. 

5.3.3 Turbulent Kinetic Energy k 

Fig. 5.14 shows the distribution of turbulent 
kinetic energy. Generally, the positions of the maxima 
and minima in the turbulent kinetic energy were not well 
defined with respect to other variables. The point of 
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Fig. 5.14 Variation of the Turbulent Kinetic Energy k 
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maximum kinetic energy falls in the range of maximum 
velocity gradients and the maximum shear stress. 
5.3.4 Correlation Functions R and a. 


The significance and the definition of the correla 
tion functions R uv and a^ are given in Chapter 4 under 
section 4.3.5. Figs. 5.15 and 5.16 show the distribu- 
tions of the correlation functions R uv and a^ respective.’ y. 
Both the functions become zero at the points of zero 
shear and they change sign according to the sign of ~uv. 
The maximum value of R uv between the outer point of zero 
shear and the freestream is around +0.57 to start with, 
gradually decreasing to about +0.27 at x/y = 100.3. 

At x/y c = 177.7, no positive values of R uV can be ob- 
served since the slant-wire data did not extend far beyond 
the outer point of zero shear. The minimum value of R 
between the inner and outer points of zero shear lies be- 
tween -0.7 and -0.55. 

Similarly, the maximum value of a^ between the 
outer point of zero shear and the freestream is around 
±0.15 "to start with, gradually decreasing to about +0.09 
at x/y c = 100.3. The minimum value of a^ between the 
inner and outer points of zero shear lies between -0.23 
and -0.17. Both R uv and a^ tend to become zero in the 
freestream. 
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5.3.5 Eddy Viscosity 

Fig. 5.17 shows the distributions of eddy viscosity 
v ^ = -uv OU/3y) . The trends of the experimental 
results are indicated by the dashed lines. It can be 
easily seen that there is a region of singularity exhibited 
at each station. As discussed earlier (section 4.3.6), 
the eddy viscosity becomes either negative or undefined 
in these singular regions. The eddy viscosity distribu- 
tion at x/y c = 100.3 exhibits two regions of singularity 
corresponding to the points of maximum and minimum 
velocity. The eddy viscosity finally drops to zero in 
the freestream. The present eddy viscosity distributi a 
do not show any similarity behavior. 

5.3.6 Mixing Length L m ^ x and Prandtl-Kolmogorov Length 
Scale L^ 

The significance and the definitions of the Prandtl 

mixing length L m ^ x and the Prandtl-Kolmogorov length scale 

Lj. were given in Chapter 4 under section 4.3.7. A 

value of = 0.2 (Kacker and Whitelaw, 1968, 1971) 

war- used in the evaluation of L^. Fig. 5.18 and 5.19 show 

the distributions of L . and L, respectively. The trends 

mix k 

of the experimental data are shown in dashed lines . As 
in the case of eddy viscosity, both the length scale 
distributions show regions of singularity. The length 
scales become either negative or imaginary in these sing- 
ular regions. The length scale distributions at x/y c = 100.3 




EWIS / tUlNFXlELTRU 


187 


03., 


02. 


dI« 

X 


U3 J 


02J 


J.£J 


OD. 




A 

\ 

Ch 


\ 

£D S 


\ 

B\ 


\ 


13 


n\ 


V \ 


Q 


\ 


\ \ / 


B “g" 


x/rn= 1SJ92 

" ~ E 


\ 


Ql 


\ 


\ 

*\ 


\ 


1UL TQ 


/ 

/ 


a * 
a < — £ 




a' 


C 

N xno3 




+ - 




177.7 


•00 


ToT 


E430 3.00 

Y CCM53 


HU30 


5*00 


5430 


Fig. 5.17 Distributions of the Eddy Viscosity 

Trends exhibited by the experimental data 
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Scale 


exhibit two regions of singularity corresponding to the 

points of maximum and minimum velocity. The length 

scales (L m1 - y and L^) beyond the outer point of zero 

shear were about three times higher than the length scales 

in between the inner and outer points of zero shear at 

x/y = 15.92. These length scales gradually become near- 
c 

ly equal at x/y = 100.3. Similarly, the length scales 
(L^x and L k ) in between the inner and outer points of 
zero shear were about three times higher than the length 
scales below the inner point of zero shear at x/y = 100.3 
and 177.7. The present length scale distributions do not 
show any similarity behavior 

5.3.7 Production of Turbulent Kinetic Energy 

The significance of the shear stress production term 
-uv ^ was given in Chapter 4 under section 4.3.8. The 
shear stress production term was evaluated using the pre- 
sent measurements . 

Fig. 5.20 shows the distribution of the shear 
stress production of turbulent kinetic energy given by 

gn 3 

-uv ^ (S/U^) on a logarithmic scale. A small region 
was observed around the points of zero shear where the 
production became negative. This is because of the fact 
that the points of zero shear do not coincide with the 
points of zero velocity gradients due to which and -uv 
assume different signs in the region around the points of 
zero shear. The negative production terms were found 
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Fig. 5.20 Distributions of the Shear Stress Production of 
Turbulent Kinetic Energy 


Trends exhibited by the experimental data 
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to be in the order of -uv (6/U^) = 1.5X10 . The 
negative production points were not shown on the plots . 
It can be seen that the production level in the layer 
between the points of zero shear is about two orders of 
magnitude higher than that of in the layer beyond the 
outer point of zero shear to start with. Gradually, 
this difference in the production levels decreases as 
one goes downstream. 



CHAPTER VI 


PREDICTION METHOD 


The objectives of the theoretical prediction of^ 
the present flows are as follows: 

1. To investigate the applicability of an exist- 
ing prediction method to the present flows 

2. To investigate the ability of the theoretical 
prediction method to show the difference in the flow 
development between a uniform jet velocity profile and 

an asymmetric jet velocity profile 

3. To establish that the asymmetric jet velocity 
profile has distinct advantages over uniform jet velocity 
profile. 

6.1 Selection of the Prediction Method 


A preliminary survey was made of the existing pre- 
diction methods for blown boundary layers in view of the above 
mentioned objectives. The integral methods, for example 
Kind (1971) , tend to be disqualified in view of the second 
and third objectives mentioned above. The asymmetric jet 
velocity profile is a detail connected with the initial 
velocity profile, and integral methods are insensitive 
to such details of the profiles, since they perform inte- 
gration over a region. In comparison, differential methods 




can incorporate the details of the initial velocity pro- 
files like the asymmetric profile into the computational 
scheme . 

Among the available differential methods, the 
method of Irwin (1974) was found to be the best suited 
for the present case of blown boundary layers. In con- 
trast to other methods such as those of Launder and 
Spalding (1972) , Ng and Spalding (1965) , Hanjalic and 
Launder (1972) , and Rodi (1972) , Irwin's method is 
particularly suited for wall jets and blown boundary 
layers and it involves relatively less empiricism. For 
most of the other differential methods, different empiri- 
cal constants are required in going from plane flows to 
axisymmetric flows, from "strong" shear flows to "weak" 
shear flows, and from plane flows to flows' with streamline 
curvature. In comparison, the empirical constants in 
Irwin ' s method were found to be independent of the above 
restrictions . 

Irwin has tested his prediction method for a 
variety of flows. They include wall jets in still air, 
blown boundary layers with pressure gradients , conventional 
boundary layers with and without pressure gradients and 
curved wall jets. Irwin's method gave accurate predictions 
of blown boundary layers with pressure gradients. The 
position of zero shear stress was predicted to be nearer 
to the wall than that of the velocity maximum in accord 
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with the experimental data. Most of the other differen- 
tial methods fail to predict this result since the 

turbulence models used in them imply that -uv = 0 at 
30 

= 0. Irwin's prediction method also gave good results 
in the case of conventional boundary layers, curved wall 
jets, and wall jets in still air. Therefore, Irwin's 
method has been chosen as the best available method for 
the prediction of the present flows. 

6.2 Description of the Theoretical Method 

The prediction method developed by Irwin essen- 
tially uses the finite difference computing method 
developed by Spalding and Patankar (1967, 1969) and the 
turbulence model of Launder, Reece, and Rodi (1973) . 

The turbulence model of Launder , Reece, and Rodi was modi- 
fied to take account of the effect of the wall on 
turbulence and the effect of streamline curvature. 

The starting point of Irwin's method is the 
Reynolds-averaged Navier-Stokes equations for incompress- 
ible fluids without body forces. For steady two-dimensional 
flows, the application of boundary layer approximations 
and the assumption of local isotropy along with the 
relevant mathematical manipulations will lead to the . 
following set of equations for the mean momentum and the 
Reynolds stresses. 


i? 
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DU _ 1 dP 

Dt p dx 
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(uv) 
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D 

Dt 


(uv) 


-7 |2 + 

9y 


E Ai + 9v , 
p 3y 9x ; 


9 (DU) _9_ 

9y P 9y 


(uv 2 ) 


(6.5) 


where 


— = u— + vA 

Dt u 9x 9y 


P and p are the mean and fluctuating components of pressure, 
c is the dissipation rate, and p is the density. 

The viscous term and the normal stresses term 
have been neglected in the mean momentum equation 6 . 1 
because of the following reasons . For points away from 
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the wall, the viscous terms were experimentally found to 
be negligible in fully-turbulent flows. Viscous terms 
are important only at points very close to the wall and 
the law of the wall velocity profile is generally applied 

3 —? ~2 

in that region. The normal stresses term — (u - v ) is 

o X 

generally found to be small in comparison with the other 
terms in the mean momentum equation.. As mentioned in 
Chapter 4 (section 4.2.3) and Chapter 5 (section 5.2.3), 
the contribution of the normal stresses term to the mean 
integral momentum equation was found to be small even in 
the present experiments. 

6.2.1 Turbulence Model 

The solution of the equation 6.1 to 6.5 requires 
that they should form a closed set. Equations 6.1 to 6.5 
do not form a closed set unless the third order correla- 
tions in the equations 6.2 to 6.5 are expressed in terms 
of Reynolds stresses and the mean velocity. Also, the 
closure of the problem requires an equation for the dis- 
sipation rate e . 

The turbulence model of Launder, Reece, and Rodi 
(1973) was used in the closure of the equations 6.1 to 6.5. 
The original turbulence model of Launder et al. (1973) is 
basically applicable to free turbulent flows. Irwin (1974) 
modified the original model to account for the wall tur- 
bulence. This modification affects mainly the pressure- 
velocity-gradient correlations in the equations 6.2 to 6.5. 
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It should be mentioned here that the later work of Launder, 
Reece, and Rodi (1975) does take account of the wall 
turbulence . 

The important features of the turbulence model of 
Launder, Reece, and Rodi (1973) are given below. This 
model is essentially a modified version of an earlier 
model by Hanjalic and Launder (1972) . The major limita- 
tion of the model is that it is applicable only to those 
flow regions where the turbulence Reynolds number is high. 
The turbulence Reynolds number R . is given by R . = 

0 A 0 A 

-~2k k 

(u ) X/v, where X = (15vu /e) . The requirement of high 

turbulence Reynolds number enables the viscous diffusion 
term in the equations for the Reynolds stresses to be 
neglected. Also, under the condition of high turbulencr 
Reynolds number, the very small eddies of turbulence 
responsible for the viscous dissipation are isotropic. 

This local isotropy condition implies that the viscous 
dissipation terms in each of the three equations for 
Reynolds normal stresses are equal to each other, taking 
the value 1/3 e. The above assumptions are already incor- 
porated in the equations 6.2 to 6.4. 

The pressure diffusion terms -,^-(£^ and 7 r-(^jin 

3y p ay p 

equations 6.3 and 6.5 were neglected since they were found 
to be small from experimental data. The triple velocity 
correlations (equations 6.2 to 6.5) were modelled by 
simplifying the exact transport equation using the 
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gradient diffusion principle. The quadruple velocity 
correlations appearing in the exact transport equation 
for the triple correlations were expressed in terms of 
the second order correlations, by assuming that the triple 
correlations are small and their distribution is Gaussian 
For the dissipation rate e, a simplified transport equa- 
tion obtained from the exact transport equation for e 
was used. A number of approximations were made in obtain- 
ing the simplified transport equation for e. 

Using the condition of approximately homogeneous 
turbulence, and Rotta's (1951, 1962) considerations, the 
pressure-velocity-gradient correlations were expressed 
in terms of the mean velocity, Reynolds stresses and e. 
However, Irwin (1974) introduced additional modifica- 
tions in the pressure-velocity-gradient correlations to 
account for the wall turbulence. 

Thus all the terms on the right-hand side of equa- 
tion 6.2 to 6.5 have been expressed in terms of the 
Reynolds stresses and the mean velocity, resulting in a 
closed set of equations 6.1 to 6.5. 

6.2.2 Computational Procedure 

Irwin (1974) has predicted a variety of flows 
including the blown boundary layers using the method 
described above. The finite difference scheme of 
Spalding and Patankar (1967, 1969) was adopted to solve 
the set of equations 6.1 to 6.5. The "source" term 
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is a terminology typical of the Spalding-Patankar scheme, 
which includes all the terms except the advection and 
diffusion for a given equation among the equations 6.1 to 
6 . r . This "source" term has been treated accordingly by 
Irwin for the Reynolds stress equations. However, Irwin 
uses an alternative procedure for the entrainment rate 
which worked well for the kind of flows he predicted. 
Details regarding the treatment of source terms, diffu- 
sion terms, entrainment and the boundary conditions 
as applicable to Spalding-Patankar ' s scheme, can be 
found in Irwin's (1974) work. 

The step size was taken as a fraction of the total 
flow width. Irwin suggested this fraction to be 0.05 
for blown boundary layers. He also suggested the number 
of grid points to be between 50 and 60. 

6.2.5 Starting Procedure 

The solution of the equations 6.1 to 6.5 using the 
computer scheme mentioned above requires the specification 
of the initial starting profiles. The profiles of mean 
velocity, Reynolds stresses and the dissipation rate are 
to be specified at the starting station as a part of the 
starting procedure. 

Irwin has used two types of starting procedures : 

1. Available experimental velocity and turbulence 
profile data as the starting profiles, and 

2. An "automatic starting" procedure. 

Both of these starting procedures were used in the present 
case in order to examine their effects on the computed 
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results. The first method is straightforward. However, 
the second method involves certain assumptions and needs 
co be studied in greater detail. Therefore, the "auto- 
matic starting" procedure is described in detail in the 
following paragraphs. 

The term "automatic starting" stems from the fact 
that it develops the initial profiles automatically using 
the given integral parameters at the jet slot as input. 

In blown boundary layers, the region near the slot really 
requires the solution of an elliptic set of equations 
rather than the parabolic set used in the present calcula- 
tions. However, the "automatic starting" procedure assumes 
that the region over which the governing equations are 
elliptic does not extend very far downstream and that fur- 
ther downstream the flow is mainly affected by the initial 
momentum of the jet, the momentum thickness of the boundary 
layer on the slot lip, the value of its form parameter, 
and the logitudinal pressure gradient. The elliptic 
nature of the flow near the slot was thus ignored and the 
parabolic equations were used right from the slot exit 
plane. A problem arises, however, because at the slot 
there is a region of zero velocity on the downstream face 
of the slot lip and the parabolic method cannot handle 
such a region. This was overcome by replacing the real 
starting profile by an equivalent one which was acceptable 
to the parabolic method. 




The "automatic starting" procedure of Irwin gen- 


erates the starting velocity profile in two following 
steps : 


1. The real velocity profile at the slot was first 
replaced with a top-hat profile in the jet and a power 

law profile in the boundary layer on the slot lip as shown 
in Fig. 6.1 (a) . The velocity Uj in the jet is such that 
the momentum of the jet is the same as the momentum in 
the real jet. The value of i$ BL and the exponent 1/n of 
the boundary layer region were such as to give the same 

value of <S 2 bl and H BL as the rea ^ P r °file. 

2. The equivalent profile shown in Fig. 6.1 (a) 
was replaced by another profile which has the same momentum 
and mass flux as the first one but contains a mixing layer 
joining the jet and the boundary layer regions. The new 
velocity profile is shown in Fig. 6.1 (b) . A cosine 
velocity profile was adopted for the mixing layer. By 
making a mass and momentum balance between the profiles 
shown in Figs. 6.1 (a) and 6.1 (b) , the entire velocity 
profile shown in Fig. 6.1 (b) can be generated. The pro- 
file shown in Fig. 6.1 (b) is the actual starting velocity 
profile used in the computations with automatic starting 
conditions . 


2 j ^ 2 j 

This procedure needs — — , j , H_ T and "a" 

y c 6 2BL BL 


as the input. Here, "a" is the fraction of the slot height 
over which the velocity profile is assumed to be uniform 
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as shown in Fig. 6.1 (b) , 6 ODT is the momentum thickness 
of the boundary layer above the slot lip, H_,_ is the form 

B.L 

parameter of the boundary layer above the slot lip and 


6 2J U J U J 

— — = — (1 - —) . Irwin recommends a value of a = 0.6. 

c 00 00 

6.2.3 (a) Starting Turbulence Profiles 

The starting distributions of the Reynolds stresses 
for the automatic procedure were obtained as follows . The 
starting profiles for the mixing layer and the boundary 
layer are given first followed by the profiles for the jet 
portion . 


The -uv profile was calculated from -uv = v e ff OU/3y) 
with empirical expressions for the eddy viscosity in the mix- 
ing layer and the boundary layer. The three normal stresses 
~~2 2 

(u , v , and w ) were set equal to each other and equal to 
| -uv I/O. 45, with their minimum values not allowed to go be- 
low the free-stream turbulence level. The dissipation rate 
was set equal to the rate of production of turbulent kinetic 
energy, i.e., e = -uv (9U/9y). 

For the jet region, -uv is equal to zero. The 
—% ~2 

three normal stresses (u , v , and w ) were set equal to 

2 

each other and equal to (0.03 U T ) . The dissipation was 

n /o J 


set equal to 


4k 


3/2 


where k is the turbulent kinetic energy. 


Here it should be noted that the starting turbulence 
profiles used in the automatic starting procedure do not 
represent the actual conditions and involve certain arbi- 
trary assumptions. However, as Irwin (1974) has indicated, 
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the wall jet flows were observed to be fairly insensitive 
to the starting turbulence conditions at the slot. 

6 . 3 Comments on the Automatic Starting Procedure 

Even though the automatic starting procedure was 
primarily developed for blown boundary layers, Irwin (1974) 
used it satisfactorily for several other types of flows 
also, such as wall jets in still air, self-preserving wall 
jets, curved wall, jets, and boundary layers. Hence, the 
automatic starting procedure was used with confidence in 
the present case backed by Irwin's satisfactory predic- 
tions using the automatic starting procedure. It may be 
noted here that for any details regarding the theoretical 
method and the computational procedure, the work of 
Irwin (1974) should be referred to. 




CHAPTER VII 


COMPUTED RESULTS 

The purpose of this chapter is to present the re- 
sults of the computations using the prediction method 
described in Chapter 6 and to , theoretically establish the 
effectiveness of an asymmetric jet velocity profile over 
the uniform jet velocity profile. 

One of the major assumptions made in the present 
predictions was that the "automatic starting" velocity 
profile described in Chapter 6 and shown in Fig. 6.1 (b) 
represents the case of a uniform jet velocity pr .file. 

In other words, a uniform jet velocity is assumed to result 
in a velocity profile shown in Fig. 6.1 (b) at a few 
slot heights downstream of the slot. The good quality of 
the predictions made by Irwin (1974) in a number of wall 
jet cases using the automatic starting procedure for a 
uniform jet velocity profile confirms that the above 
assumption is a valid one. 

The present computed results are classified into 
three categories: 

1. Predictions of the zero pressure gradient flow 
were made using the experimental initial velocity and 
turbulence profiles and were compared with the experimental 
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data. The results were then compared with the computa- 
tions using the automatic starting procedure for a uniform 
jet velocity profile with the same jet momentum as the 
experimental asymmetric profile at the slot. 

2. Predictions of the adverse pressure gradient 
flow were made using the experimental initial velocity 
and turbulence profiles and were compared with the ex- 
perimental data. The results were then compared with the 
predictions using the automatic starting procedure for a 
uniform jet velocity profile with the same jet momentum 
as the experimental asymmetric profile at the slot. 

3 . Comparison of the computations for the ad- 
verse pressure gradient flow using the automatic starting 
procedure with two different types of jet velocity pro- 
files; uniform jet velocity profile and a hypothetical 
linear jet velocity profile described later in this chap- 
ter to represent the ideal case of an asymmetric jet ve- 
locity profile. 

7.1 Prediction of the Zero Pressure Gradient Flow 

7.1.1 Computations Usir.'’ the Experimental Initial 
Velocity and Turbulence Profiles 

The computations were done using the prediction 
method described in Chapter 6. The experimental velocity 
and turbulence profiles at x/y c = 16.5 were given as the 
starting profiles. The station x/y c = 16.5 was chosen as 
the starting station because it was the first station 




I 


Fig. 7.1(a) Initial Profile of -uv for the Computations Using 

Experimental Initial Profiles in the case of Zero Pressure 
Gradient Flow 

experimental distribution; assumed distribution 
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1. First normal hot-film data point where 


— C.C 

-uv r 



2. Point of maximum velocity where -uv = 0 

3. Point of maximum |-uv| 


4. First slant-wire data point 


5. Point of minimum velocity 




downstream of the slot and away from the slot where slant- 

wire turbulence measurements were available. 

— S 

7.1.1 (a) Initial Mean Velocity and u Profiles 

The measured normal hot-film data of the mean 
~ 2 

velocity and u at x/y c = 16.5 were used as the initial 
profiles. The grid points coincided with the data points 
in number and position. 

~ 2 “7 

7.1.1 (b) Initial Profiles of v , w , -uv, and s 

T ~2 

The slant-wire data of v , w , and -uv were in- 
terpolated to obtain their values corresponding to the 
position of the normal hot-film. However, the slant-wire 
data did not extend as close to the wall as the normal 

hot-film data. Therefore, the following assumptions were 

~ ~~2 

made for the profiles of v , w and -uv between the first 

normal hot-film data point and the first slant-wire 

~2 ~2 ~ 2 

data point, v and w were set equal to u . A linear 

variation of -uv was assumed with -uv = 0 at the point 
of maximum velocity. Fig. 7.1 (a) shows the assumed -uv 
distribution. The dissipation rate e = 0.3k(3U/3y) was 
assumed as suggested by Irwin (1974) , where k is the tur- 
bulent kinetic energy. The step size was taken as .03 
times the flow width. 

7.1.1 (c) Results 

The computed velocity profiles are compared with the 
experimental data in Figs. 4.4 (c-e) . The computed 
profiles were shown only at selected representative stations 
to prevent overcrowding. The predictions compare well with 
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the data. Figs. 7.1 (b-d) show the predictions of 6^, ($ 2 / 
and Cf compared with the data. The agreement between the 
e. cperimental data and the predictions is generally good. 
However, the predicted values of 6^ and 6 2 are lower than 
the data at far downstream stations and the predicted 
values are lower than the data by about a maximum of 10%. 

Fig. 7.1 (e) shows the predicted development of 

profile length scales, Y max , Y half , and Y min . Fig. 7.1 

(f) shows the predicted variation of U__«, and U .. compared 

max nun 

with experiments. Here also the agreement is good. 

7.1.2 Computations Using the Automatic Starting Procedure 

The predicted results using the automatic starting 
procedure for zero pressure gradient flow are shown in 
Figs. 7.1 (b-f) along with the expert wnr-al data and the 
computations using experimental starting profiles. A 
value of a = 0.6 was used as suggested by Irwin. The 
fractd * the slot height over which the velocity is 
assumed t<- 2 uniform is "a," as shown in Fig. 6.1 (b) . 

The number of grid points used was 60. The step size was 
chosen as .015 times the flow width up to x/y c = 15.0 
and .03 times the flow width after x/y c = 15.0. 

The computations using the automatic starting 
procedure did not show any major difference as compared 
to the computations using the experimental initial profiles. 
On the basis of the arguments given in Chapter 4 , it was 
expected that the growth rate of Y max would be lower, while 
the growth rates of Y^ a ^^, Y m ^ n , 6^ and 6 2 would be higher 
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Fig. 7.1(b) Predicted Development of the Displacement Thickness for 
the Zero Pressure Gradient Flow 

O Experimental Data 

Predictions using the experimental initial profiles 

Predictions using the automatic starting procedure 


Fig. 7.1(c) Predicted Development of the Momentum Thickness for the 
Zero Pressure Gradient Flow 

O Experimental Data 

Predictions using the experimental initial profiles 

Predictions using the automatic starting procedure 
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Fig. 7.1(d) Predicted Variation of the Skin Friction for the Zero 
Pressure Gradient Flow 


O Experimental Data 


Predictions using the experimental initial profiles 


Predictions using the automatic starting procedure 
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in the case of automatic starting procedure. This is due 


to the unique behavior of an experimental asymmetric jet 


velocity profile as compared with a uniform jet velocity 


profile used in the automatic starting procedure. The 


predictions showed that the growth rate of Y max distribu- 


tion was tending towards the universal Y distribution 


discussed earlier in Chapter 4. However, the predicted 


growth rates of and Y were lower instead of being 

higher. The predicted growth rates of 6^ and S 2 were al- 


most the same as those predicted using experimental 


initial conditions 


The present discrepancies between the predictions 


using experimental asymmetric jet velocity profile and 


the predictions using uniform jet velocity profile (auto- 


matic starting procedure) might be due in part to the 


inability of the prediction method to show the difference 


in uniform and asymmetric profiles when no pressure grad- 


ient was imposed on the flow. This is further confirmed by 


the predictions in the case of adverse pressure gradient flow 


The superiority of an asymmetric profile can be 


justified on the following physical grounds even though 


it has not been revealed by the predictions. The asym- 


metric profile has lower velocity gradients at the wall 


than a uniform velocity profile and this feature results 


in reduced momentum losses due to friction at the wall. 


The asymmetric profile has a higher concentration of 


momentum in the upper portion of the slot. This feature 


'iV vy W 
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provides greater momentum to the momentum deficient 
area of the upstream boundary layer and results in the 
even distribution of momentum. As mentioned later in 
section (7.3.3), the mass flow rate in the case of an 
asymmetric profile is less than that of a uniform profile 
for the same total jet momentum. This feature results in 
savings in mass flow rate in the case of an asymmetric 
profile . 

7.1.2 (a) Universal Y__„ Distribution 

££* * 

According to Narayan (1973) , the value of the 
2 2 

quantity S 2s = ^2BL U “ ^ ^ U Ja ve ^c^ l° w to observe 

the universal Y distribution mentioned in Chapter 4. 

max 

Accordingly, a few prediction runs were made using the 

automatic starting procedure with U Javg and y c values 

higher than in the experiments , thereby reducing the value 

of 6 - . The results show that the Y m=v distribution does 
2s max 

not follow the universal Y„,„ distribution when the momen- 

max 

turn deficit of the upstream boundary layer is relatively 
higher. The value of 6 2s was of the order of 0.48 in 
the present experiments . 

7.2 Predictions of the Adverse Pressure Gradient Flow 
7 .2.1 Computations Using Experimental Initial Profiles 




7.2.1 (a) Initial Mean Velocity and u Profiles 


The measured normal hot-film data of the mean 
~ 

velocity and u at x/y = 15.92 were used as the initial 

C 

profiles. The grid points coincided with the data points 
in number and position. , 

~7 ~2 

7.2.1 (b) Initial Profiles of v , w , -uv and e 

~J ~ 2 

The initial profiles of v , w , -uv and e were 
obtained in the same way as given in section 7.1.1 (b) . 
Fig. 7.2 (a) shows the assumed -uv distribution. The step 
size was taken as 0.03 times the flow width. 

7.2.1 (c) Results 


The predicted velocity profiles were compared 

with the experimental data in Figs. 5.4 (b) and 5.4 (c) . 

The agreement is good up to x/y = 72.14. The predic- 

c 

tions deviate from the experimental data from x/y = 

100.3 onwards. However, the predictions were good in the 
outer layer even after x/y = 100.3. The predictions 

G 

were not shown beyond x/y = 177.7 because the outer edge 
of the outer layer was already touching the top wall at 
x/y ~ 177.7. However, the computations proceeded all 
the v/ay up to x/y c = 360.8 without any instabilities in 
the calculations. The predicted velocity profile at 
x/y„ = 360.8 was found to be a normal one, without any 
regions of large momentum deficit. 

Figs. 7.2 (b-h) show the various predicted integral 
and profile parameters compared with experimental data. 
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Fig. 7.2(b) Predicted Development of the Displacement Thickness for 
the Adverse Pressure Gradient Flow 


© Experimental Data 

— : — Predictions using the experimental initial profiles 
Predictions using the automatic starting procedure 
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Fig. 7.2(c) Predicted Development of the Momentum Thickness for the 
Adverse Pressure Gradient Flow 

O Experimental Data 

Predictions using the experimental initial profiles 

1 Predictions using the automatic starting procedure 
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Fig. 7.2(d) Predicted Variation of Skin Friction for the Adverse 
Pressure Gradient Flow 

O Experimental Data 

Predictions using the experimental initial profiles 

Predictions using the automatic starting procedure 




Fig. 7.2(e) Predicted Variation of U for the Adverse Pressure 
Gradient Flow maX 

O Experimental Data 

Predictions using the experimental initial profiles 

Predictions using the automatic starting procedure 





Fig. 7.2(g) Predicted Variation of Y . for the Adverse Pressure 
Gradient Flow 111111 

O Experimental Data 

Predictions using the experimental initial profiles 

Predictions using the automatic starting procedure 
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Fig. 7.2(h) Predicted Development of the Profile Length Scales for 
the Adverse Pressure Gradient Flow 

•, fl, O Experimental Data 

Predictions using the experimental initial 

profiles 

Predictions using the automatic starting 

procedure 
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Here it should be noted again that the experimental data 
beyond x/y = 177.7 must be treated with caution because 

u 

the entire tunnel height was filled with the shear layer 
beyond x/y c = 177.7. 6 ^ and.- ^ were predicted well up to 

x/y c = 177.7. The predictions of are lower than the 
experiments. The predicted U max decay is faster than the 
experiments. However, the decay of U m £ n and the growth 
rate of Y m ^ n are predicted well. The agreement between 
the predicted and experimental growth ratet of Y max and 

Y half 5 00( ^ U P to x/y c = 143.0. 

The disagreement between the predictions and 
the experiemntal data at far downstream stations is 
mainly due to the fact that beyond x/y = 177.7, the entire 
tunnel height was filled with the shear layer and hence 
the full thickness of the shear layer has not been taken 
into account in the experiments . 

7.2.1 (d) Convergence or Divergence Correction 

Irwin (1974) suggested that for flows under severe 
adverse pressure gradients, a correction term involving 
3W/9z has to be applied to the mean momentum equation to 
account for the flow convergence or divergence in the 
lateral direction. Here z is the lateral coordinate and 
W is the lateral component of velocity. The flow con- 
vergence or divergence is caused by the rapid growth of 
the side wall boundary layers and their bleed off by means 
of false side walls. The fact that the flow converges or 
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divei’ges in the lateral direction depends on whether too 
little or too much of the side wall boundary layer flow 
was being bled. 

Sufficient care was taken in the present experi- 
ments to maintain two-dimensionality. The slant-wire 
measurements (Appendix C) show that the magnitude of W is 
very small even at far downstream stations. This implies 
that the magnitudes of 9W/9z may not be significant to 
apply any convergence or divergence correction. Also, 
it was shown in Chapter 5 (section 5.2.3) that the ex- 
perimental data closely satisfy the two-dimensional integral 

momentum equation up to x/y = 177.7. The failure of the 

c 

experimental data to satisfy the two-dimensional integral 
momentum equation beyond x/y c = 177.7 is mainly due to 
the fact that the entire tunnel height is filled by the 
shear layer beyond x/y c = 177.7. Because of the above 
arguments, no convergence or divergence correction has 
been applied to the predictions presented in this thesis. 
Also, one of the main objectives of the predictions is' to 
compare the performance of the asymmetric jet velocity 
profile with that of a uniform jet velocity profile. Since 
it is only a relative comparison, the flow convergence 
or divergence correction should not influence the conclu- 
sions regarding the asymmetric jet velocity profile. 


7.2.2 Computations Usinq the Automatic Startinq Procedure 


MOW 
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in Chapter 6. As mentioned before, the velocity profile 
(Fig. 6.1 (b) ) used in the automatic starting procedure 
was assumed to represent the case of the uniform jet 
velocity profile. The step size was chosen as .015 times 
flow width up to x/y = 15.0 and .03 times the flow width 
beyond x/y = 15.0. The number of grid points was taken 

w 

to be 60. The value of "a" was taken as 0.6. 

The computations using the automatic starting 
procedure predicted reverse flow near the minimum velocity 
point in the wake region at x/y = 130.0. The computations 
broke down after the reverse flow was observed since the 
computing method cannot handle regions of reverse flow. 

In comparison, the computations using the experimental 
starting conditions representing the asymmetric jet 
velocity profile proceeded all the way up to x/y c = 360.8. 
Fig. 5.4 (b) shows the predicted velocity profile at 
x/y c = 100.3 using the automatic starting procedure in 
comparison with the predicted profile using experimental 
starting conditions. A large momentum deficit can be found 
in the minimum velocity region of the predicted velocity 
profile using the automatic starting conditions. Fig. 7.4 
(a) shows a representative predicted velocity profile using 
the automatic starting procedure with uniform jet velocity 
profile just before the occurrence of reverse flow in the 
minimum velocity region. 

Figs. 7.2 (b-h) show the predicted development of 
different integral and profile parameters using the automatic 










starting procedure in comparison with the predictions 


using the experimental starting conditions. The momentum 


and displacement thicknesses shown in Figs. 7.2 (b) and 


(c) increase very rapidly near the point where the 


computations broke down. One can see that the increase 


of 6^ and 6^ starts from the beginning of the computations, 


even though the increment is more rapid near the point of 


breakdown. The skin friction near the point of breakdown 


of the computations shows a trend that the flow might 


encounter separation at the bottom wall, if one proceeds 


further downstream. 


Figs. 7.2 (e-g) show the decay of U max and U m £ n 


and the growth rate of Y m ^ n respectively . U m ^ n increases 


rapidly from the start of the computations and eventually 


goes to zero near the point of breakdown. Y . increases 

min 


rapidly as compared to the predictions using experimental 


starting conditions. Fig. 7.2 (h) shows the development 


of Y , Y . , and Y, Y . _ and Y, - - increase very 

max mm hair mm halt J 


rapidly for the case of automatic starting conditions, 


Y max 9 rowt h r ate is lower than that for the experimental 


starting conditions. 


The above predicted results using the automatic 


starting procedure can be explained in the following way. 


When a pressure gradient is imposed on the flow, the region 


around the junction of the upstream boundary layer and the 


jet needs a considerable amount of momentum to overcome the 


pressure gradient. If sufficient momentum is not supplied 
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to that region by the jet, the result is a highly momentum 
deficient region near the minimum velocity point at the 
downstream stations. In the case cf uniform jet velocity 
profile used in automatic starting procedure, consider- 
able part of the jet momentum is concentrated near the 
wall and not close to the minimum velocity region. Hence 
the outer layer spreads more rapidly into the free-stream 
in search of momentum as indicated by the rapid increase 
of However, because of the severe adverse pressure 

gradient, the entrainment of the free-stream momentum does 
not occur rapidly enough to supply momentum to the mo- 
mentum deficient region near the point of velocity 
minimum. Hence, the velocity in the minimum velocity 
region decreases rapidly as indicated by the rapid decay 

of U . and eventual reverse flow there. Because of the 
mm 

existence of a large momentum deficient region near the 
point of minimum velocity, the momentum and displacement 
thicknesses increase very rapidly. The relative concen- 
tration of the jet momentum near the wall causes increased 
frictional losses, resulting in the rapid decrease of C^. 

The relatively large momentum and displacement thicknesses, 
the lower values of C^, and reverse flow in the minimum 
velocity region are all unwanted effects and they might 
load to eventual flow separation at the bottom wall. 
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In comparison, the flow development in the case 
of the experimental asymmetric jet velocity profile is 
different. Given the same total jet momentum, the asym- 
metric jet velocity profile will have a higher concentra- 
tion of momentum in the upper half of the jet. Hence, the 
momentum requirements of the minimum velocity region is 
more readily met by the higher momentum in the upper half 
of the jet. The outer layer does not grow very rapidly 
as indicated by the slower growth rates of and Y m ^ n 

shown in Fig. 7.2 (h) . Y__„ increases rapidly because of 

the tendency of the inner layer to grow faster to supply 


momentum to the outer layer. The decay of U m ^ n shown in 


Fig. 7.2 (f) is more gradual in comparison to the rapid 


decrease of U m ^ n to zero in the case of predictions using 


automatic starting procedure. The growth of displacement 
and momentum thicknesses is also gradual. In essence, the 
asymmetric jet velocity profile tries to distribute the 
momentum evenly across the layer, thereby preventing the 
possibility of separation of the inner layer near the wall. 
The asymmetric velocity profile also results in lower fric- 
tional losses at the wall by keeping the jet momentum away 
from the wall as indicated by the lower values of in 
the region near to the jet. 

Here it should be mentioned that Irwin (1974) has 
documented one case of flow prediction where the computations 
predicted reverse flow in the minimum velocity region with 
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the experimental starting conditions, even though the ex- 
periments did not show any flow reversal in the minimum 
velocity region. Irwin gave the explanation that the 
predicted flow reversal in the minimum velocity region 
may be due to either experimental inaccuracy or the failure 
of the boundary layer assumption because of the very low 
velocities reached in the region of velocity minimum. 

It was further indicated that the application of flow con- 
vergence-divergence correction, along with the accurate 
specification of the experimental starting conditions has 
carried the computations farther than without them. But, 
the prediction method still predicts the reverse flow. 

In the present case, it is true that the flow velocities 
are considerably smaller, at far downstream stations. How- 
ever, the difference in flow development between the 
uniform velocity profile and asymmetric velocity profile 
can be observed from the beginning of the computations. 
Therefore, the breakdown of computations in the present 
predictions using the automatic starting procedure is not 
a local effect of low velocities, but it is a result of 
a gradual process which starts at the slot because of 
the uniform jet velocity profile. 

7.3 Comparison of the Performance of Linear and 
Uniform Jet Velocity Profiles 

This part of the predictions has been made to 
isolate the effect of the shape of the jet velocity profile 
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on the flow development. Predictions were made with two 
different types of jet velocity profiles in the case of 
an adverse pressure gradient flow by keeping the other 
initial conditions the same. 

The two jet velocity profiles considered are 
given below. A uniform jet velocity profile shown 
in Fig. 7.3 (a), which is the same as the profile shown 
in Fig. 6.1 (b) , and was generated by the automatic 
starting procedure and has the same jet momentum as the 
experimental case. A linear jet velocity profile shown 
in Fig. 7.3 (b) was the second profile. This linear jet 
velocity profile was assumed to represent the ideal case 
of an asymmetric jet velocity profile. It was obtained 
by replacing the profile shown in Fig. 7.3 (a) below 
the velocity minimum point by two linear segments with the 
following constraints : 

1. U • is the same for the profiles shown in 

mm 

Fig. 7.3 (a) and 7.3 (b) 

2. U = 0 at y = 0 

3. U = U. at y/y = a, where a is the fraction 

Jmax J J c 

of the slot height over which the flow is uniform 

4 . The momentum below the minimum velocity 
point is the same for both cases. 

Referring to Fig. 7.3 (b) , U___ v is greater than U • . 

j max j 

However, for the particular case where a =1.0 and there 
is no mixing layer, Uj max =*' ^3 Uj as shown in Fig. 7.3 (c) . 
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Here, U___„ is the maximum jet velocity in the linear 

JluaX 

profile and Uj is the jet velocity in the case of uniform 
profile for the same jet momentum in both cases. The 
above mentioned velocity profiles shown in Figs. 7.3 (a) 
and 7.3 (b) were used as the starting velocity profiles 
in the predictions. 

7.3.1 Selection of the Parameter a 

a is the fraction of the slot height over which 
the flow is uniform with reference to Fig. 7.3 (a) . 

Irwin recommends a value of a = 0.6. When a = 0.6 was 
used, the resulting linear profile was such that the ve- 
locity gradients in the mixing layer were lower than the 
velocity gradients in the jet region below the point of 
maximum velocity . This results in greater momentum 
transfer towards the wall than towards the region of min- 
imum velocity. However, in an actual situation of an 
asymmetric profile, the region near the slot can be 
expected to have higher velocity gradients in the mixing 
layer than in the jet region below the point of maximum 
velocity. Therefore, the value of "a" was taken as 0.95 
in generating the profiles shown in Figs. 7.3 (a) and (b) . 
The higher value of "a" makes the width of the mixing 
layer very thin thereby Increasing the velocity gradients 
in the mixing layer and simulating the actual conditions 
that might be expected in the case of an asymmetric jet 
velocity profile. 
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7.3.2 Starting Turbulence Profiles 


. 2 2 2 

The starting profiles of u , v , w , -uv and t 


used in the predictions for the uniform jet velocity profile 


were kept the same as for the automatic starting procedure. 


2 2 2 

The starting profiles of u , v , w used in the predictions 


for the linear jet velocity profile were also kept the 
same as for the automatic starting procedure. 

The starting profile of -uv for the linear 
velocity profile was chosen in the following way. The 
-uv distribution above the point of minimum velocity 
was kept the same as for the automatic starting procedure. 
The shear stress in the mixing layer was obtained from the 
shear stress distribution in the mixing layer for the 
automatic starting procedure by multiplying with a factor 


equal to <0^ - Wlinear /(Dj - Wunifom • The 


shear stress distribution below the point of maximum 
velocity was taken as linear with -uv = 0 at the point of 
maximum velocity and -uv = v ( 3D/9y) ^i near at tJie first 
grid point. 

The number of grid points was taken to be 80 
to accommodate the details of the linear velocity profile. 
The step size was taken as 0.015 times the flow width 


up to x/y = 15.0 and 0.03 times the flow width after x/y = 

G G 


15.0. 




7.3.3 Prediction Results 


As mentioned earlier, the predictions using linear 
and uniform profiles were made only for the case of ad- 
verse pressure gradient flow. The theoretical method 
predicts flow reversal in the minimum velocity region for 

the case of uniform profile at x/y = 102.0. In compar- 

c 

ison, the computations for the case of linear profile 
proceeded all the way up to x/y = 360.8 without any 
problem. Fig. 7.4 (a) shows a comparison of the predicted 
velocity profiles for the uniform and linear cases at 
x/y c = 100.3. One can see a large momentum deficient 
region near the minimum velocity point for the case of 
uniform profile. 

The predicted development of the different integral 
and profile parameters are shown in Figs. 7.4 (b-g) 
for both the uniform and linear cases . Most of the pre- 
dictions made in the case of uniform profile are in 
agreement with the results obtained in Section 7.2.2. 

The momentum and displacement thicknesses increase very 
rapidly for the uniform profile in comparison with the 
linear profile. The skin friction shown in Fig. 7.4 (d) 
decreases very rapidly indicating higher momentum loss 
due to friction at the wall for the uniform case. In fact, 
the area between the two skin friction curves up to the 
point where the for both cases become equal indicates 
an increased loss of momentum due to friction for the case 
of uniform profile. 
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Fig. 7.4(a) Predicted Velocity Profiles at x/y = 100.: 
for the Adverse Pressure Gradient ^low 


Predictions with linear jet velocity 

profile 

Predictions w itu uniform jet velocity 

profile 






Fig. 7.4(b) Predicted Development of the Displacement Thickness for 
the Adverse Pressure Gradient Flow 


Predictions with linear jet velocity profile 

Predictions with uniform jet velocity profile 

and with the same momentum as in the experiments 

— • — • — Predictions with uniform jet velocity profile and 

with 20% excess momentum than in the experiments 

— x — x — Predictions with uniform jet velocity profile and 

with 30% excess momentum than in the experiments 
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Fig. 7.4 (e) 


Predicted Variation of U for the Adverse Pressure 
Gradient Flow m n 


Predictions with linear jet velocity profile 

Predictions with uniform jet velocity profile 
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The decay of U • and the growth of Y.„ are shown 
mm 3 min 

in Figs. 7.4 (e) and 7.4 (f) . u m £ n decreases very rapidly 

to zero and Y increases considerably in the case of 

uniform profile. In comparison, the decay of U m ^ n and the 

growth of Y m ^ n are gradual for the linear profile. Fig. 

7.4 (g) shows the relative development of the profile 

length scales Y„, = . Y, and Y . . Once again Y. - - and 

3 max' half' mm J half 

Y . „ develop much faster for the uniform case than for 
the linear case. However, the growth rate of Y max is 
almost the same for both cases even though the absolute 
value of Y max is lower in the case of uniform profile. 

The mechanism of flow development in the case of 
uniform jet velocity profile in comparison with that for 
the asymmetric jet velocity profile is given in detail 
under Section 7.2.2 and the same thing is applicable when 
a comparison is made between the uniform and linear jet 
velocity profiles. However, the flow development with a 
linear profile and with the same starting conditions as 
the uniform profile reveals more clearly the distinct 
advantages of having an asymmetric velocity profile in 
the jet. 


At this point, it is reasonable to ask the question 
of how much extra jet momentum is needed in the case of 
uniform profile to achieve a similar flow development 
as in the case of the linear profile indicated by the 
development of different integral and profile parameters. 
To investigate it, predictions were made with different 


‘r''^ •v *’^ f>V-- * '•;“ .W\ '•["■■ . »" : '*. •• ; *- :, >V. / 


values of jet momentum higher than the experimental value. 


The flow development was compared on the basis of develop- 


ment of 6 1 and S Figs. 7.4 (b) and 7.4 (c) show the 
predicted development of 6 and ^ using the automatic 


starting procedure with 20% excess momentum and 30% 


excess momentum. With 20% excess momentum, the predic- 


tions broke down at about x/y_ = 177.70. Reverse flow 

G 


was predicted in the minimum velocity region near x/y = 

c 


177.7. With 30% excess momentum, the flow was predicted 


all the way up to x/y = 360.84 without any reverse flow. 

c 


Therefore, an approximately 30% increase in jet momentum 


is necessary in the case of an uniform profile to achieve 


similar flow development as in the case of a linear 


velocity profile. Also, it should be noted that for a 


given jet momentum, the mass flow rate in the case of 


asymmetric profile is less than the uniform case. This is 


evident from the mass flow calculations for the linear 


and uniform cases shown in Fig. 7.3 (c) . The mass flow 


for the linear case is t/T/2 times the mass flow for the 


uniform case with the same momentum. However, for a sim- 


ilar flow development, the momentum of the uniform jet 


has to be increased by 30%. Therefore the actual savings 


in mass flow in the case of a linear profile for a similar 


/T l 

flow development is (1 - -y- -=pr ) x 100 or 24% along with 


a jet momentum savings of 30%. Fig. 7.4 (d) shows the 


predicted skin friction using the automatic starting 
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procedure with 20% excess momentum and 30% excess momentum 
The increased momentum loss due to friction in these cases 
can be observed as indicated by the high values of skin 
friction. 







CHAPTER VIII 


CONCLUSIONS 

Conclusions 

1. Measurements have been made in a zero pressure 
gradient two-dimensional wall jet with a low jet excess 
velocity and a thick upstream boundary layer. The quan- 
tities measured are the mean velocity, turbulence 
intensities, Reynolds shear stress, spectra and skin 
friction. The derived quantities included turbulent 
kinetic energy, eddy viscosity, mixing length, Prandtl- 
Kolmogorov length scale, correlation coefficients, dis- 
sipation and production rates of turbulent kinetic energy 
and bursting periods. 

2 . Measurements have also been made in an adverse 
pressure gradient two-dimensional wf.l' je with a low jet 
excess velocity and thick upstream boundary layer. The 
quantities measured and derived included ail those listed 
under (1) with the exceptions of spectra, dissipation rate 
and bursting period. 

3. The mean velocity profiles show a defined 
logarithmic region which followed the conventional log- 
arithmic law of the wall with the same constants as 
applicable to normal boundary layers . 
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4. The point of zero shear stress was generally 
found to be closer to the wall than that of zero velocity 
gradient. 

5. The profiles of effective viscosity, Prandtl 
mixing length and the Prandtl-Kolmogorov length scale do 
not exhibit an obvious similarity and consequently pre- 
diction procedures which depend on these concepts are 
unlikely to be wholly satisfactory for blown boundary 
layers . 

6. Tchen's mean vorticity theory with a = 0.8 
estimates the dissipation rates which are in good agree- 
ment with -5/3 law in the logarithmic region of the velocity 
profile. 

7. It is possible to produce an asymmetric jet 
velocity profile and it is a unique feature of the present 
wall jet design. 

8. The asymmetric jet velocity profile affects the 
downstream development of the flow considerably. It causes 
the effect of the jet to be carried to a much longer dis- 
tance as compared to a uniform jet velocity profile. The 
asymmetric profile tries to supply more momentum to the 
momentum deficient upstream boundary layer instead of 
wasting it as friction at the wall, resulting in a more 
even distribution of momentum across the whole layer. 

9. The prediction method of Irwin (1974) gave 
satisfactory predictions of the present class of blown 
boundary layers using experimental starting conditions. 


10 . The prediction method used was able to show 
theoretically the advantages of asymmetric jet velocity 
profile over uniform jet velocity profile more clearly 
in the case of adverse pressure gradient flow. 

11. Blown boundary layers under adverse pressure 
gradients with an asymmetric jet velocity profile result 
in much lower values (or more accurately lower growth 
rates) of momentum and displacement thicknesses and pre- 
vent the development of large momentum deficient region 
downstream of the slot as compared to a uniform jet veloc- 
ity profile. It is well known that large displacement and 
momentum thicknesses result in a situation approaching 
separation and cause lower lift on aerofoils. 

Suacrestions for Future Work 


There is a need for more detailed and accurate 
turbulence measurements taking proper account of high flow 
inclination angles for the present class of flows under 
adverse pressure gradients. A parametric study of the 
effect of asymmetric velocity profile for different slot 
heights and jet velocities can be made. Also, there is 
a need for comparison of measurements with a uniform and 
an asymmetric jet velocity profile under similar conditions 
with the same jet momentum in each case. The skewness of 
the asymmetric jet velocity profile can be enhanced further 
and its effect can be studied on the downstream flow devel- 


opment . 




APPENDIX A 


MEASUREMENTS ON A WALL JET BY BOWLES (1977) 


Bowles (1977) took hot wire measurements of U 

and at the exit of a wall jet with a similar design as 

in the present experiments . The slot height in his case 

was 6.35 mm. Fig. A.l shows the velocity and turbulence 

prof iles at the exit of the jet for the Bowles case. 

The asymmetric velocity profile typical of the present 

wall jet design can be observed in this case also. How= 

ever, the turbulence profile is different as compared 

with the turbulence profile at the jet exit in the present 

experiments shown in Fig. 2.11. The turbulence profile 

shown in Fig. A.l shows a maximum near the wall at about 

0.2y which is well above the outer limit of the thin 
c 

sublayer for this high speed flow (U Jmax = 35 m/sec) . 

The cross sectional area of the nozzle around the corner 
in the case of Bowles' jet was kept constant as compared 
to a reduction in the cross-sectional area in the present 
case. As discussed in Chapter 2, this feature gives rise 
to the possibility of separation around the corner and the 
consequent re-attachment downstream. Therefore, the maxima 
in u t near the wall shown in Fig. A.l may be the remainder 
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APPENDIX B 


ESTIMATION OF THE FRICTIONAL LOSSES 


IN THE WALL JET NOZZLE 


The details on the estimation of frictional 
losses in the present wall jet nozzle are given in this 
appendix. The frictional losses are estimated by making 
an energy balance on the control volume surrounding the 
wall jet nozzle. The measured jet exit velocity profile 
and the measured static pressure inside the nozzle are 
used in the estimation of frictional losses. 

Fig. B.l shows the close-up view of the wall jet 
nozzle. The section (1) corresponds to the position where 
the static pressure connection was made in the wail jet 
nozzle c.s described in Chapter 2. Section (1) is located 
1.27 cn below the surface of the wind tunnel as shown in 
Fig. 2.5. The Section (2) corresponds to the exit of the 
wall jet. Since most of the frictiona'l losses in the noz- 
zle occur between Sections (1) and (2) , an attempt has 
been made to estimate those losses. 

Let be the velocity at any given y 

location at Section (2) 


Uj be the uniform average velocity at 
Section (1) 






b-^ be the width of the nozzle at Section (1) 
be the static pressure at See cion (1) 

above the static pressure at Section (2) 
p be the density of the fluid 
For the unit length of the nozzle: 

Kinetic energy going out at Section (2) = 



Kinetic energy coming in at Section (1) = 


Flow work input at Section (1) = 


p i b i u i 


Let the frictional losses between the Sections 
(1) and (2) = F.L. 

An energy balance on the control volume surround- 
ing the Sections (1) and (2) and the portion of the nozzle 
in between them gives 


F.L. 




dy 


i °i b i + p i b i 


U, 


(B.l) 


In the equation (B.l), can be obtained from the measured 
velocity profile at the jet exit, and the corresponding 
value of P^ can be measured directly. The velocity at 
Section (1) is obtained by making a mass balance between 
the sections (1) and (2) as follows: 

Let M-^ and M ? be the mass flow rates per unit 
length of the nozzle at Sections (1) and (2) respectively. 
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Then 


M i - PD i b i 


“* c 

M,, = p J D. 

t. Jr\ * 


dy 


(B . 2) 


can be calculated from the measured distribution of U 2 
and the mass balance = M 2 . Therefore, the frictional 
losses (F.L) in the equation (B.3) can be calculated from 
the measured velocity distribution U 2 at the jet exit 
and the corresponding measured static pressure (P^) 
at Section (1) . 

In the present nozzle design, the frictional losses 
between the Sections (1) and (2) were found to be about 
15% of the total energy input at Section (1) . These 
frictional losses between Sections (1) and (2) form a 
significant part of the total frictional losses in the 


nozzle. 


APPENDIX C 





MEASUREMENT OF TURBULENCE INTENSITIES 
AND THE REYNOLDS SHEAR STRESS 

This appendix gives the details regarding the 
measurement of turbulence intensities v./U and w. /U 

t 00 L. 00 

and the Reynolds shear stress -uv using the slant-wire. 

A brief description of the probe alignment procedures, 
and the estimation of mean V and W are also given. 

The directional sensitivity of the slant-wire 
makes possible the use of it for the measurement of in- 
dividual components of the velocity vector. Several 
investigators, e.g., Andersen, et al. (1972), have used 
the slant wire successfully to obtain different components 
of the turbulent fluctuations and the Reynolds shear stress 
-uv. A similar method was used in the present experiments 
for the determination of v^., w^, and -uv. A general out- 
line of the method of obtaining these turbulent quantities 
will be given here and for greater details one may refer 

to the work of Andersen et al. (1972) . The method of ob- 
~2 T 

taining v , w and -uv is given below. v t and w t can be 

2 ~~2 2 ~2 
obtained using the relations v t = v and w t = w . 

I 

j 

i 

: 

j 
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C.l Theory of the Directional Sensitivity of the Hot-Wire 
As shown by Jorgensen (1971) , the effective cool- 
ing velocity for a hot-wire may be approximated by 

u t 2 = (u 2 .) 2 + (k 1 ) 2 (v 2 ) 2 + (k 2 ) (w 2 ) 2 (C.l) 

where U 2 is the velocity component perpendicular to the 
wire and lies in the plane of wire supports, V 2 is parallel 
to the wire and W 2 is perpendicular to the wire and the wire 
supports. and K 2 are constants which depend on the hot 

wire probe design and the velocity. These constants are 
experimentally determined for the given probe and the range 
of velocities under consideration. 

Fig. C.l shows the probe geometry and position 
for a slant hot-wire probe. The slant wire makes an angle 
<f> with a plane perpendicular to the probe axis. (x^, y^, z^ 
is the mean flow co-ordinate system in which the mean veloc- 
ity vector has the components (U-^, 0, 0). This co-ordinate 
system is in general not completely coincident with the 

laboratory co-ordinate system (x, y, z) in which the mean 

~2 

velocity vector is (U, V, W) . The measured values of v , 

~2 

w , and -uv will strictly speaking be expressed in the 
(xi, y^, Zl ) co-ordinate system. However, for the flow 
under consideration V and W are small compared with U and 
the measured turbulent quantities can be considered as 
measured in laboratory co-ordinates (x, y, z) with negli- 
ble error. The hot-wire probe axis will be assumed to be 
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aligned with the mean velocity vector, i.e., the probe 
axis is in the direction of x^. 0 is the angle of rotation 

of the probe about its own axis and it is zero when the 
plane of wire supports coincides with the x^, y 1 plane. 

Referring to Fig. C.l, the relationship between 
the velocity components V^, W.^) in the mean flow 

reference system (x^, y^ z 1 ) and the components (U 2 , V 2 , W 2 ) 
in the wire co-ordinates is given by 


1 

K ai 
to 

1 


cos<!> sin<)>cos0 sin<}>sin0~ 


- _ 
U 1 


^2 

— 

-sin<j> cos<}>cos0 cos<|>sin0 


r>u 

V 1 

(C.2) 

U 2 J 


0 -sin0 cos0 


L wj 



Equation (C.3) is obtained by performing the following steps. 

1. Inserting equation C.2 into C.l. 

2. Resolving the velocity components as 


U. = U. + u. 

l l l 


D i - °i + »i 
V x = v x 


W 1 ■ W 1 


3. Recognizing that the mean flow reference 


# ^ ~ ~ 


(CJ^, v r W x ) is the same as the laboratory frame of refer- 
ence (U, V, W) for the flow under consideration where V and 
W are small compared to U. 

4 . Simple mathematical manipulation of the re- 
sulting equations. 


5. Neglecting the higher order terms. 
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4 


(C. 3) 


~2 T T 

Where u , v , w , and -uv are the normal stresses and shear 


T . 


stress in the laboratory frame of reference and u^ is the 
indicated or measured normal stress at any given position 


of the probe. A, D, F are the coefficients given by 


2 2 2 
A = cos <j> + sin <J> 


2 . 


D = (1-K^ ) sin2<J>cos0 
2> 


F = (1-K-^ ) sin2<|)sin0 

For the slant wire used in the present measurements <j> = 45' 


2 2 

C.2 Determination of v , w , and -uv 


Equation (C.3) contains four unknowns, i.e., u , 


T 

v , w , and -uv. Therefore, a minimum of four measurements 


of uf are required with the rotatable slant-wire at four 


~7 ~Z 

different values of 0 for the determination of u , v , w , 


and -uv. However, the nature of the coefficient A does 


not permit the evaluation of u from equation C.3. This 
is because the value of A does not change with different 
values of 0 for a given value of <p , when the probe axis 


is aligned with the mean flow direction. Hence, the normal- 


hot film data of u were used in the equation (C.3) to ob- 


~2 ~2 ~2 

tain v , w and -uv. The u at the slant wire location 


was obtained from the normal hot-film data by interpolating 


~2 


a five point quadratic curve fit. With u being known, a 


minimum of three measurements of u| are required with the 
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rotatable slant wire at three different values of 6 for the 

. J 2 ~ 2 

determination of v , w , and -uv. Measurements of u^. were 

made at seven different values of 6 = (m-l)7r/6 (m = 1,2,3, 

. . .7). Seven measurements were made instead of the min- 
imum required three, to make a redundancy check on the 

~ ~2 

quantities v , w , and -uv. 

C.2.1 Techniques of Obtaining v ^ , w ^. , and -uv 

The basic principle used in the calculation of 
the turbulence quantities was to obtain them in a way in- 
volving a minimum amount of mathematical manipulation of 
the direct measurements in order to reduce the errors to 
a minimum. 

-uv was obtained from equation C.3 in three 
different ways from six different measurements of u^ as 
follows: 


1. 

( “i* e=o° - 

(u i J 0=180° 

2. 

^ U i^ 0=30° 

(u i ) 0=150 

3. 

(u i ) 0=6O° 

(u i } 0=120 


All three values of -uv agreed very well within a maximum 
of 10% variation. The civer age of the three values of -uv 
was finally taken. 

~2 ~2 
v was obtained from the measurements of (u^)q_q 0 

~ 2 

and { u i)0 = i8O°* Th: *- S set of measurements eliminates the 

~2 ~ 2 T 

effect of w on v data. A comparison of v obtained 

~2 

in this manner with the v obtained from measurements at 










other values of 0 gave good agreement between them with 


a maximum variation of about 15%. 


2 2 
w was obtained from the measurement of (u^)g_gQ 0 . 


2 

This measurement eliminates the effect of v and -uv on 


T T 

w data. A comparison of w obtained in this manner with 


the w obtained from measurements at other values of 0 gave 


good agreement between them with a maximum variation of 


about 15%, 


~2 ~7 

The variations in v , w , and -uv obtained from 


different sets of measurements are mainly due to the mean 


velocity gradients present within the finite space occupied 


by the slant wire and they cannot be calibrated. It should 


be noted here that the uncertainties in the measurement 


~2 

of v , w , and -uv given in Chapter 3 are smaller than 


the variations due to velocity gradients mentioned above. 


C.3 Determination of Constants and in the Equation C.l 


Knowledge of is necessary to obtain the tur- 


bulence quantities from equation C.3. However, an attempt 


has been made to determine the value of K 2 also. The method 


of Jorgensen (1971) was used for the determination of and 
K 2 » It essentially involves the use of equation C.l and 


the calibration of the slant wire at three different posi- 


tions in the calibrator. The flow in the calibrator described 


in Chapter 2 has the components of velocity (U,0,0) in the 


coordinate system (x-^, y-^, z^) shown in Fig. C.l. The rela- 


tionship between the non-linearized anemometer output Voltage 


E and the effective velocity U. can be written as 
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E 2 = A + B OjH (C.4) 

Where n is the exponent of the order of 0.45 and A, B are 
constants. Calibration of the slant wire was done at 
three different positions given below: 

1. The wire is perpendicular the flow and faces 
it with the plane of the wire supports in line with the 
flow, with the corresponding non-linearized anemometer 
voltage output being Eq . 

2. The axis of the probe is in line with the flow 
and the wire faces the flow, with the corresponding output 
being E-^ . 

3. The plane of the wire supports, the wire and 
the probe axis are all perpendicular to the flow with the 
corresponding output being E^. 

K-^ and were evaluated from the following equations ob- 
tained from the equation (C.4). 


K. 


K. 



2 . 

COS <P 


(C. 5) 


With <|> = 45‘ 


■ ,J S 


2 . 
cos 0 


(C. 6) 


With 6 = 90° 


The determination of and K 2 was done at different 
flow velocities covering the range of velocities encountered 
in the wind tunnel. The value of K-^ used in the equation 
(C.3) was chosen according to the magnitude of the mean ve- 
locity at the data point under consideration. The value of 


r 

’(1 

Bjft' 


In the case of zero pressure gradient flow, the 
alignment was done in the free- stream. To do the alignment, 
the probe was held in the free-stream and then the pitch 
and yaw of the probe were adjusted so that the mean anemo- 
meter output remained the same at any angular position of 
the probe when rotated about its own axis . Afterwards , 
the probe position was not disturbed during the process of 
taking data for that streamwise station. This alignment 
procedure gave rise to a probe position where the probe 
axis was parallel to the free-stream velocity. 

It was not possible to precisely adjust the probe 
alignment with the mean flow vector in the boundary layer 
and to precisely measure the angular change with respect to 
the free-stream alignment. Hence the mean anemometer out- 
puts were different at different values of 0. The 
misalignment of the probe axis with the mean flow in the 
vertical and horizontal planes can be calculated from 










the mean anemometer outputs at 0 = 0°, 90° and 180° using the 
misalignment analysis of Andersen et al. (1972) „ In the 
region of large velocity gradients, part of the difference 
in anemometer output at different values of 0 should be 
attributed to velocity gradients and finite size of the 
wire. However, in the present case, it was assumed that 
the difference in the anemometer outputs a: 0 = 0°, 

90° and 180° is entirely due to the misalignment of the 
probe. The misalignment angle was never greater than 0.8° 
in the vertical plane, and was never greater than 1° in 
the horizontal plane for the majority of the slant wire 
data in the case of zero pressure gradient flow at x/y c = 

16.5 and beyond. As given by Andersen et al. (1972), the 
error introduced in the measured turbulence quantities 
due to these small misalignment angles is very small. 

The misalignment angles for stations x/y c = 0.292 and -16.4 
were relatively higher because the flow near the ^all 
jet body is inclined at these stations. This is due to 
the special shape of the wall jet body and no attempt has 
been made to correct the data at these stations for the 
flow inclination. 

C.4 (b) Probe Alignment for the Adverse Pressure Gradient Flow 
In the case of the adverse pressure gradient flow, 
the stream-lines were inclined towards the top wall since 
the fluid was escaping through the perforations in the top 
wall. The flow inclination was maximum near the top wall 
and it decreased gradually towards the bottom wall. It 




is not possible to precisely adjust the alignment of the 
probe with the flow at each vertical location. Therefore, 
for reference, the probe was aligned with its axis parallel 
to the bottom wall. 

The misal: gnment of the probe axis with the flow 
in the horizontal and vertical planes was calculated in 
the same way as for the zero pressure gradient flow. The 
misalignment angle in the horizontal plane was found to 
be small and of the order of 1°. However, the misalign- 
ment angle in the vertical plane which is also equal to 
the flow inclination angle was found to be large when 
compared to the zero pressure gradient flow. 

For points near the free-stream, where the shear 
stress was very small, the flow inclination angle reached 
values up to 13°. However, in the region where the tur- 
bulence quantities were significantly large, the flow 
inclination was less than 5°. 

Here, it should be noted that the constant K-^ 

varies with the angle of inclination of the flow vector 

with the probe axis. Littlefield (1978) indicated that 

~ J 

good measurements of v within about 10% error car. be 
made by neglecting higher order terms and using constant 
values of K-^ in the equation (C.3) for flow inclination 
angles up to 5°. The error in the -uv measurements is rela- 
tively higher. In the present measurements the flow inclination 
was less than 5° points where the shear stress is significant. 


An attempt has been made to estimate the mag- 
nitudes of mean V and W with respect to the probe to a 
first order approximation using equations (C.l) and (C.2). 
The principle applied in the calculation of V and W was 
that inside the boundary layer the difference in mean 
anemometer out r uts at different angular positions of the 
probe about its axis gives an indication of the magni- 
tudes of V and W. The same principle was also applied in 
the calculation of misalignment angles. It is easier to 
conceive that the misalignment of the probe in the vertical 
and horizontal planes is a measure of mean V and W respec- 
tively . 

A x'edundancy check was also made on the V and 
W by calculating them from three independent measurements. 

In the case of the zero pressure gradient flow, for stations 
x/y = 16.5 and beyond, the values of W were found to vary 
at random within a maximum of ±.02U and the values of V were 
found to be within a maximum of +.02U. 



In the case of adverse pressure gradient flow, 
the flow was inclined towards the top wall and hence re- 
sulted in a larger V as compared to the zero pressure 
gradient flow. However, the nature and magnitude of W 
variation were comparable with that of the zero pressure 
gradient flow. 



APPENDIX D 

In this appendix the normal hot-film data for the 
zero pressure gradient flow are tabulated. The data con- 
sist of mean velocity U/U^ and the turbulence intensity 
u t /U oo . The skin friction data are also presented at each 


station. 
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TABLE D 

NORMAL HOT-FILM DATA OF U/U. AND u t /D„ FOR 
THE ZERO PRESSURE GRADIENT FLOW 


X/YC - -77.JOO 

UINF - 

CT » .003305 

WXM 

Y (CHS) 

U/UIHF 

.015 

.338 

.018 

.369 

.020 

.403 

.023 

,429 

.025 

. 451 

.028 

.470 

.033 

.200 

• 038 

.522 

.043 

.543 

* .046 

.557 

.053 

,570 

*050 

.579 

.063 

.39: 

.074 

.408 


.102 

.114 

. 17 

• 140 
.143 
.191 
.214 
.341 
.2 47 

• 3'JO 

• 394 
.437 
.321 

• 394 
.711 
.938 
.9*3 

1.092 

1.219 

1.344 

1.473 

1.400 

1.934 

2.100 

2.342 

2.414 

2.070 

3.370 

3.604 

4.394 

4.902 

3.410 

3.919 

4.424 

4.400 


.424 
.434 
.443 
.433 
. 443 
.479 
.492 
.704 
.714 
.724 
.747 
.749 
,707 
.003 
.021 
.030 
.002 
.900 
.933 
.934 
.973 
.904 
.974 
1.000 
1.000 
1.000 
.999 
.999 
.999 
.990 
.998 
. .998 
.997 
.997 
.997 
.9 99 


23.47 H/SEC 
» 1.329 CflS 


.0037 
.0082 
,09 02 
.0924 
.0921 
.0913 
.0007 
.0049 
.0041 
*0014 
.0001 
.0733 
.0771 
.0747 
.0737 
.0720 
.0732 
.0727 
.073a 
.0733 
.0739 
.0738 
.0739 
.0732 
.0723 
.0700 
.0494 
.0479 
.0444 
.0434 
.0394 
.0330 
.0497 
,0404 
.0304 
.0199 
.0119 
.0063 
.0041 
.0030 
.0024 
• 0020 
.0016 
*0014 
>0013 
.0012 
.0012 
.0012 
.0013 
*0014 


X/rc • “16.340 UJNF - 28.21 H/8EC 


DELTA - 1.130 Cftfl 


y<CMB) U/UIHF UT/UIHF 


. 01 ! 

.016 

.020 

.023 

.020 

*033 

.030 

.043 

.040 

.033 

.061 

.074 

.064 

.099 

.112 

.137 

.163 

.201 

.264 

.320 

.391 

.453 

.518 

*643 

.772 

.899 

1.024 

1.260 

1.534 

1.706 

2.042 

2.294 

2.550 

2.604 

3.050 

3.105 

3.312 

3.620 

4.320 

4.034 

5.344 

5.652 


.401 
.579 
.655 
.716 
.730 
.794 
.030 
.663 
.679 
.691 
.697 
.903 
.909 
.914 
.916 
.919 
.921 
• 920 
.922 
.922 
.924 
.929 
.933 
.939 
.944 
. 33 
.963 
.972 
*961 
.995 
1.000 
.995 
.909 
s9 62 
.970 
.974 
.970 
.966 
.966 . 
.961 
.937 
.934 
.952 
.930 


*0634 
.0761 
* 0031 
.0630 
.0622 
.0787 
.0724 
.0676 
*0645 
.0622 
.0602 
.0566 
.0564 
.0345 
.0537 
.0323 
.0322 
.0526 
*0524 
.0528 
.0527 
.0519 
.0505 
.0466 
.0474 
.0463 
.0444 
.0415 
.0377 
.0233 
.0132 
.0069 
.0041 
.0029 
.0023 
*0019 
.0016 
.0016 
.0015 
.0012 
• 0011 
.OO&O 
.0012 
.0016 


Y(CftS) 

.020 

• 023 
.025 
.020 
.030 
. 033 
.030 
.043 
.040 
.053 
.061 

• 066 
♦ 079 
.091 
.104 
.117 
.142 
.160 
.193 
.206 
.210 
.231 
.244 
.249 
.234 
.239 
.264 
.269 
.274 
.204 
.292 
.300 
.307 
.313 
.323 
.320 
.333 
.338 
.343 
.333 
.356 
.366 
.371 
.384 
.396 
.409 
.434 
.460 

. 4B5 
.311 
.349 
.361 
.625 


UJNF - 27.14 H/SEC 
DELIA • 2.174 L*HB 

8 1 .2.61 

U/UIHF UT/UlNF 


.732 
,813 
.942 
1.069 
1.194 
1.323 
1.577 
1.831 
2. 085 
2.342 
2.596 
2.047 
3.101 
3.409 
4.117 
4.623 
3.J33 
3.641 
6.150 


1.203 
1.225 
1.239 
1.233 
1.243' 
1.276 
1.2V3 
1.306 
1.312 
1.370 
1.23: 
1.337 
1.317 
1 .340 
1.36 7 
1.372 
1.303 
1.392 
1.377 
1 . 393 
1 .300 
1.309 
1.262 
1.166 
1.000 
.743 
.393 
.168 
.129 
.136 
. 0*0 
.040 
.043 
,072 
.159 
.236 
.301 
.379 
.427* 
,449 
.445 
« 4B3 
.494 
.521 
.541 
.557 
.507 . 
.414 
.639 
• 6ol 
.491 
.702 
.737 
.762 
.700 
.795 
.821 
,049 
.076 
.396 
.934 
.945 
.985 
.994 
.999 
1.000 
.999 
.999 
.996 
.995 
.995 
.993 
» 993 


-0909 

.0090 

.0071 

.oo/-» 

.0023 
.0 797 
.0754 
.07)4 
.0673 

• U643 

* 0o07 
.0507 
,0540 
■ 0222 
.0500 
.0479 
.0459 
.0454 
.0495 
.6546 
, 0646 
, Ob3L‘ 
.1103 
.1101 
. 1267 

. 1573 
,1263 
.0542 
.0309 
.0360 
.0341 
.0103 
.0152 
.0303 
.0508 
.0725 
.0801 
.0842 
. 0050 
.0854 
• 0858 
.0864 
.0043 
.0859 
.0054 
.0846 
.0034 
.0810 
.0007 
.0790 
.0761 
.0752 
.0712 
.0687 
.0673 
.C *ol 
.0432 
.0416 
.0380 
.0537 
.0452 
.0315 
.0173 
.0081 
.0049 
.0036 
.0027 
.0019 
.0016 
.0013 
.0011 
.0011 
.0014 
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TABLE D (Cont.) 


X'TC - 7.«30 UINF • 27. 99 fl/SEC 

Cf • .013294 DELTA - 2.298 Ctt8 


X/tC - 14.300 UIHF - J4.47 H/BEC 

W - .010934 DELTA - 2.232 CHB 


Y(CHS) 

U/U1NF 

UT/UIi 

.010 

.713 

.1443 

.013 

.749 

. 1328 

.013 

.831 

. 1391 

.018 

,914 

.137? 

.020 

.970 

. 1358 

.023 

1.017 

.1498 

.028 

1.083 

.1410 

.033 

1.129 

.1324 

.038 

1.147 

.1232 

.043 

1.197 

• 1178 

.040 

1.219 

.1128 

.033 

1.737 

,1073 

.044 

1.270 

.0997 

,079 

1.290 

.0930 

.091 

t .301 

.0919 

.104 

1 .304 

.0912 

,117 

1.304 

.091? 

.130 

1.301 

.0930 

.142 

1.293 

.0931 

.133 

1.283 

.0974 

.148 

1.249 

.1002 

• 180 

1.231 

.1021 

.193 

1.231 

.1047 


.244 
.349 
.295 
.320 
.343 
.371 
.394 
.422 
.447 
.472 
.498 
.323 
.349 
.374 
.399 
.423 
.400 
.732 
.SIS 
.079 
1.004 
1.133 
1.240 
1.3Q7 
1.314 
1.441 
1.740 
1.093 
2.022 
2.149 
2.274 
2.403 
2.330 
2.437 
2.7B4 
3.030 
3.292 
3.344 
3,00 0 
4.300 
4.014 
3,324 
3.032 
4.340 


1.127 
1.044 
1.001 
,937 
.070 
.029 
.70S 
.738 
.734 
.723 
.721 
.722 
.724 
.732 
.741 
.740 
.748 
#704 
.793 
.008 
.033 
.839 
.881 
,902 
.922 
.940 
.933 
.947 
.970 
.983 
.989 
. 993 
.993 
.994 
.994 
.997 
.998 
.999 
.999 
1.000 
1.000 
1.000 
1.000 
1.000 


.1112 
.1187 
.1253 
,1206 
.1273 
. 120V 
.1100 


.013] 

.0094 

.0071 

.0057 

.0043 

.0038 

.0029 

.0023 

.0018 

.0014 

.0013 

,0011 

.0010 

.0011 

.0014 


Yf CHS) 

.010 

.013 

• 013 
*016 
.020 
,023 
.028 
.033 
.041 

• 048 
.034 
.069 
.081 
.094 
.107 
.119 
.132 
.143 
.137 
.170 
.183 
.200 
,234 
.239 
.310 
.373 
.437 
.300 
.34/ 
.42 


.0873 

.734 

.0707 

.BIO 

.0733 

.801 

.0712 

1.008 

.0702 

1.133 

.0703 

1.262 

.0702 

1 .522 

.0490 

1.774 

.Q481 

2.030 

.0443 

2.284 

.0433 

2.338 

. 0442 

2.792 

.0637 

3.044 

.0414 

3.343 

.0382 

3.481 

.0354 

3.990 

.0313 

4.232 

.0443 

4,740 

,0410 

5.268 

0340 

3.776 

0267 

4.284 

0195 


Lt/UINF UT/UIAT 


.794 

.1242 

.042 

.1243 

.888 

.1236 

.924 

.1704 

.938 

.1173 

.983 

.1130 

1.024 

.1078 

1.034 

.1034 

1.092 

.0908 

1.119 

.0738 

1.137 

.0931 

1.143 

.0904 

1.181 

.0888 

1.194 

. 0BB1 

1.202 

.0877 

1 .203 

.0874 

1.204 

.0880 

1.204 

.0087 

1.200 

. 0898 

1.193 

.0908 

1.183 

.0919 

1.163 

.0931 

1.140 

.0982 

1.111 

.1010 

1.049 

.1043 

.967 

.1021 

.897 

.0923 

.043 

.0794 

.007 

.0493 

.791 

.0433 

.794 

.0437 

.802 

.0633 

.814 

.0433 

.023 

.043? 

.847 

,0407 

.870 

.0574 

.892 

.0544 

.929 

.0433 

.940 

.0341 

.981 

.0203 

.993 

.0102 

.997 

.0038 

1.000 

.0041 

1.000 

.0031 

1.000 

.0023 

1.000 

.0020 

1.000 

*0010 

1.000 

.0017 

.999 

.0013 

1.002 

.0014 

1.003 

.0014 

1.002 

.0018 


X/TC » 33.030 UINF - 24.64 H/BEC 

CF - .007797 DEITA - 2.1V2C«C 
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APPENDIX E 


In this appendix the slant-wire data for the 
zero pressure gradient flow are tabulated. The data 
consist of turbulence intensities v./U , w./U and the 

X, 00 U 00 

2 u t 

shear stress -uv/U^ . The interpolated data of and 

U oo 

Sd’/U^)/ 3 (y/6) at each slant-wire location were also 
tabulated. They were obtained from the normal hot-film 
data by interpolating a five point quadratic curve fit. 
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APPENDIX F 


In this appendix the normal hot-film data for the 
adverse pressure gradient flow are tabulated. The data 
consist of mean velocity U/U^ and the turbulence intensity 
u t /U ro . The skin friction data are also presented at each 
station . 
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APPENDIX G 


In this appendix the slant-wire data for adverse 
pressure gradient flow are tabulated. The data consist 
of turbulence intensities v t /U oo , w t /U oo and the shear 
stress . The interpolated data of u./U and 

ui 

3 (U/U^) /3 (y/S) at each slant-wire location were also 
tabulated. They were obtained from the normal hot-film 
data by interpolating a five point quadratic curve fit. 


291 



! 

1 

I 

1 

I 


f 

\ 


! 

i: 


i 

i 


o 

w 

fa 


§ 


Eh 


w 


w 


fa 

TtT 


> 

3: 

Q 

O 

*5 

fa 


fa 

fa 


w 

Eh 

Eh 

a 


fa 

fa 

H 

C 

Q 

fa 

g 

C 

e> 

Eh 


C 

fa 

Q 

fa 


D 

0 

W 

fa 

CO 

H 

0 

S 

fa 

1 

fa 

EH 


£ 


< 


fa 


W 





X o 

3 x 


ooooooooooooooo oonrv 
_ ooooooooo il ts « \) n in fv « m 


«Nnnnrtn--ar 


1 «r -<mp cm •* 


OOOOOOOOOOOOOOOOOO^ 
O Mh o o o OOOMnOOn^ ■*■ n (N <N 
— t ^ooofN'Oajcoo’O 


Ncoorfop 


Na»>MHnNriH 4 NoonTN«« 


*Co 
X o 

M O 

*2 

> X 


OOOOOOOOOOOOOOO O *0 O* M 
oooooo.oooors>ooo«4<or>«n>w 
m <o *• >o m rv jrvinnirt'riN-oe'-nc'ne-o 


s o 

3 x 


ooooooooooooooo o o «o is 

ooooooooos]Min»>N>«o<rvr>r>t 
•»nrj < * , -»-*o.ojnor<PNrx^r>. D>on t 

m -i f3 n <* « *■ n rJ oJ ts ck o w w ri m 


«« 0 ' 0 »TnoifintoMionH|fMv 

OMT-O-ONMBSiiaTfiOlflNnhOO 

OM*OOONNCQ«ronmOfSf<rincri 

n nn<n 4ri o *T4n«n ’ t* t 

I MMNNNNNh 
I I I I I I • I 


v m *o o n o ps tm mo>n Mnn xgkMn 
ror< 03 *T-«n--co.f\p>o(a:n<-r) <o<no 

NnrtTnsiMsoNnNOinoinono 


w o 

2 o 

3 x 


ooooooooooeoooooonram 

ooocQTnvioMiHinnMiofinon 

NOO»"t»CDOn'OOD>Vlll«'<IMn MM 

mitv^Mnnniii^innvrtw 


*2 

2 x 


ocooooooooooooooo^^n 

oooooNPinwwN'onnoTHMK 

fso^*fso^rxo**orsf^i>ir»nnn<^rjnp« 

minh TNOTnnnTT<ttMn« 


OOOOOOOOOOOOOOO O T'O N H 
O 0 O 0 OOQ 9 O 0 OOD>O* 0 «rM(NOQ 
Ooooo-«*ro»r*«cr\jMOt>»ra ▼mo a o 

© M «r n n Ifi ri tn *o «g ri ► *o ri m 

▼ ♦ *r n Px M MMMMM 

T 77 7 1 1 

ooooooooo oooaoooo^r o*o 
ooo oona>»or.p>nnMMfs ■oo-oMn 

OJ.ttPIH'OwPlONTftO'CMfflM'PIW 


d o a nn 4 <4 •rm * 


52 

3 X 


M*.nokXrv< 4 CMttv»m<»^fSMa'eMoa 

N»oo<iNPiOTmnDno>' 4 riOMpin 

n«M 9 NDTa<o H nN»Tnn 4 innNi , i 

iioan'o I’tnonffiiw- 
t n m 

i i i i i i 


n'B'40»'onor\*MiftM'«o*(t)'«vn<a 
Mrnfvn'Oo.oivN.on—'oM^ m rr *• »-> <o 
«P»M ▼«HOB>O.OC\*T'*0*rM»rsr>INf'»OI 

" m m m m m ri ri ri in ri *r *> 


292 



1 


*1 


> 


'.t. 


% 


> 

! 



1 






LIST OF REFERENCES 


Andersen, P.S., Keys, W.M. and Moffat, R.J. 1972. 

The turbulent boundary layer on a porous plate : 

An experimental study of the fluid mechanics for 
adverse free-stream pressure gradients. Report No. 
HMT-15, Thermosciences Division, Dept, of Mech. 

Engrg. , Stanford University. 

Bradshaw, P. 1967a. The turbulent structure of equi- 
librium boundary layers. J. Fluid Mech. 29: 

625-645. 

. 1967b. Conditions for the existence of an 

inertial subrange in turbulent flow. Nat. Phys . 

Lab., Aero. Rept. No. 1220. 

, and Gee, M.T. 1962. Turbulent wall jets with 

and without an external stream. Aero. Res. Counc. 

R and M No. 3252. 

Bradshaw, p. and Pankhurst ,k-C. 1964. The design of low 

speed wind tunnels. Progress in aeronautical sciences , 
vol. 5, pp. 1-69. New York: Macmillan. 

Bowles, E.B. 1977. Design, construction and testing of 

wall jet-wall suction boundary layer control systems. 
Internal report WT-5, Dept, of Civil and Mech. 

Engrg. , Southern Methodist University. 

Carriere, P., Eichelbrenner , E. , and Poisson-Quinton . 

1959. Contribution theoretiqur et experimentale 
a 1' etude du controle de la coache limite par 
soufflage. Advances in Aeronautical Sciences , 
vol. 2. New York: Macmillan. 

Coles, D.E. 1962. The turbulent beoundary layer in a 
compressible fluid. Rand Corporation Report 
R-403-PR. 

Corrsin, S. 1964. Further generalizations of Onsager's 
Cascade Model for Turbulent Spectra. Phys. of 
Fluids 7: 1156. 


294 


295 


English, E.C. 1970. Ph.D. dissertation, Cambridge 
University. 

Gartshore, I.S., and Newman, B.G. 1969. The turbulent 

wall jet in an arbitrary pressure gradient. Aero- 
nautical Quarterly 20: 25. 

Goradia, S.H. and Colwell, G.I. 1971. Parametric study 
of a two-dimensional turbulent wall jet in moving 
stream with arbitrary pressure gradient. AIAA 
Journal 9: 2156. 

Guitton, D.E. 1970. Some contributions to the study 

of equilibrium and non-equilibrium turbulent wall 
jets over curved surfaces. Ph.D. dissertation, 

McGill University. 

Hanjalic, K. , and Launder, B.E. 1972. A Reynolds stress 

model of turbulence and its application to thin shear 
flows. J. Fluid Mech. 52: 609. 

Hinze, J.O. 1959. Turbulence , First Edition. New York: 
McGraw-Hill Book Co. 

. 1975. Turbul ence , Second Edition. Ne T -» York: 

McGraw-Hill Book Co. 

Irwin, H.P.A.H. 1973. Measurements in a self-preserving 
plane wall jet in a positive pressure gradient. 

J. Fluid Mech. 61: 33. 

. 1974 Measurements in blown boundary layers 

and their prediction by Reynolds stress modelling. 
Ph.D. dissertation, McGill University. 

Jordinson, R. 1961. Design of wind tunnel contractions. 
Aircraft Engrg. October: 294. 

Jorgensen, F.E. 1971. Directional sensitivity of wire and 
fiber-film probes. DISA Information, No. 11, May. 

Kacker, S.C., and Whitelaw, J.K. 1968. Some properties 
of the two-dimensional turbulent wall jet in a 
moving stream. Trans . ASME , J. Applied Mechanics , 
December: 641. * ~ ~ 

. 1971. The turbulence characteristics of two- 

dimensional wall- jet and wall-wake flows. Trans . 
ASME , J. Applied Mechanics , March: 239. 


296 


Kind, R.J. 1971. A calculation method for boundary 
layer control by tangential blowing. C .A .S . I . 
Trans. 4: 88. 


Klebanoff, P.S. 1954. Characteristics of turbulence 
in a boundary layer with zero pressure gradient. 

NACA Tech. Notes No. 3178. 

. 1955. Characteristics of turbulence in a 

boundary layer with zero pressure gradient. NACA 
Rpp. 1247, 1135-1153. 

Kline, S.J., and McClintock, F.A. 1953. Describing 

uncertainties in single-sample experiments. Mech. 

Eng . 75:3-8. 

Kruka, V., and Eskinazi, S. 1964. The wall jet in a 
moving stream. J. Fluid Mech. 20: 555. 

Launder, B.E. 1969. The Prandtl-Kolmogorov Model of 
Turbulence with the inclusion of second-order 
terms . Trans. AS ME , Journal of Basic Engineering , 
December: 855. 

Launder, B.E., Reece, G.J., and Rodi , W. 1973. Dev- 
elopment and application of a Reynolds stress 
turbulence closure. Imperial Colleae, Mech. Eng. 
Dept., HTS/73/31 . 

. 19 75. Progress in the development of a Reynolds 

stress turbulence closure. J. Fluid Mech. 60: 

537. 


Launder, B.E., and Spalding, D.B. 1972. Lectures in 
Math matical Models of Turbulence . New York : 
Academic Press. 

Littlefield, M.T. 1978. Investigation of hot-wire 

response characteristics in skewed mean flow. M.S. 
thesis. Dept. Mech. Eng., University of Washington. 

McGahan, W.A. 1965. The incompressible turbulent wall 
jet in an adverse pressure gradient. Gas Turbine 
Lab. Rep. No. 82, Mass. Inst, of Technology. 

Narayan, K.Y., and Narasimha, R. 1973. Parametric 
analysis of turbulent wall jets. Aeronautical 
Quarterly 24: 207. 

Newman, B.G., and Irwin, H.P.A.H. 1975. The prevention 
of separation by blowing in two-dimensional flow. 
AGARD Conference on Flow Separation, Conference 
Pre-Print No. 168. 


297 


Ng, K.H., and Spalding, D.B. 1969. Some applications 
of a model of turbulence for boundary layers near 
walls. Imperial College of Science and Technology, 
Dept, of Mech. Eng. Rep. BL/TN/A/14. 

Nicoll, W.B., and Ramaprian, B.R. 1970. Performance of 
conical diffusers with annular injection at inlet. 
Trans. ASME, Ser. D. , Journal of Basic Enqineerincj 
92: 827. 

Papell, S.S., and Trout, A.M. 1959. Experimental in- 
vestigation of air film cooling applied to an 
adiabatic wall by means of an axially discharging 
slot. NASA TN D-9 , August, 

Pa' al , R.P. 1962. Self-preserving, two-dimensional 

turbulent jets and wall jets in a moving stream. 

M.Sc. Thesis; McGill University. 

Patel, V.C. 1965. Calibration of the preston tube and 
limitations on its use in pressure gradients. 

J. Fluid Mech. 23: 185. 

Ramaprian, B.R. 1969. Concial diffusers with annular 
injection at inlet. Ph.D. dissertation, Mech. 

Eng. Dept., Univ. of Waterloo, Canada. 

, 1973. Turbulent wall jets in conical diffusers. 

AIAA Journal 11: 1684. 

. 1975. Turbulence measurements in an equilib- 
rium axisymmetric wall jet. J. Fluid Mech. 71: 317. 

Rotta, J.C. 1951. Statistiche Theorie Nichthomogener 
Turbulenz. 2 ■ Phys . 129: 547. 

. 1962. Turbulent boundary layers in incompress- 

Tble flow. Progress in Aeronautical Sciences , vol. 

2. New York: Pergamon Press. 

Samuel, A.E., and Joubex't, P.N. 1964. Film cooling of 
an adiabatic flat plate in zero pressure gradient 
in the presence of a hot mainstream and cold tangen- 
tial secondary injection. ASME Paper No. 64-WA/HT-48. 

Seban, R.A. 1960. Heat transfer and effectiveness for 
a turbulent boundary layer with tangential fluid 
injection. Trans. ASME , J. of I-Ieat Transfer , 

Series C, 82: 303. 


298 


Simpson, R.L., and Wyatt, W.G. 1972. The behavior of 
hot-film anemometers in gas mixtures. Internal 
report, Dept, of Civil and Mech. Eng., Southern 
Methodist University. 

Spalding. D.B., and Patankar, S.V. (1967 , 1969) Heat and mass 
transfer in boundary layers. London: Morgan- 

Grampain, Parts I and II. Part III: Imperial 

College, Dept. Mech Eng., Report No. BL/TN/G/15. 

Strickland, J.H., and Simpson, R.L. 1973. The separat- 
ing turbulent boundary layer: An experimental 

study of an airfoil type flow. Rept. WT-2, 

Thermal and Fluid Sciences Center, Southern Methodist 
University. Available from NTIS as AD - 771170/8GA. 

. 1975. Bursting frequencies obtained from wall 

shear stress fluctuations in a turbulent boundary 
layer. Phys. Fluids 18: 306. 

Thomas, F. 1962. Untersuchungen iiber die Erhohung 
des Auftriebes von Tragflugeln mittels Grenz- 
schichtbeeinf lussung durch Ausblason. Z. Flugwiss 
10: 46. 


. 1965. Boundary layer control for increased 

lift by blowing. AIAA J. 3: 967. 

Ueda, H., and Hinze, J.O. 1975. Fine structure turbulence 
in the wall region of a turbulent boundary layer. 

J. Fluid Mech. 67: 125. 

Wieghardt, K. 1946. Hot Air Discharge for De-Icing. 

AAF Trslation No. F-TS-919-RE. 

Williams, J. , and Alexander, A.J. 1958. Pressure plot- 
ting measurements on a 8 per cent thick aerofoil 
with trailing edge flap blowing. A.R.C., R and M 
3087. 


